rat 


ret o~_ 
OF MICHIGAN 
FEB 1 21993 


ENGINEERING 
LIBRARY 














»Vouume 20 FEBRUARY, 1953 





CONTENTS 


On the Motion of a Tumbling Body 
Slender-Body Theory—Review and Extension. ... Mac C. Apams AND W. R. SEARS 


Penetration and Deflection of Jets Oblique to a General Stream... Frepric F. Earicn 


An Experimental Study of a Hypersonic Wind-Tunnel Diffuser 
PETER P. WEGENER AND R. KENNETH LOBB 


Analysis of Stiffened Curved Panels Under Shear and Compression 
M. A. MELCON AND A. F. ENsrupD 
Reverse-Flow and Variational Theorems for Lifting Surfaces in Nonstationary Com- 
pressible Flow 
Differential Shear-Flow Analysis for Nonprismatic Semimonocoque Beams 
WruiaM B. KING AND WALTER R. GARRISON 
Readers’ Forum: ‘Some Remarks on ‘A Solution of the Laminar Boundary-Layer 
Equations for a Compressible Fluid with Variable Properties, Including Dissocia- 
tion,’ ’’ by M. Z. Krzywoblocki; ‘‘Some Remarks on the Effect of Flapping Hinge 
Offset on Rotor Blade Stall,” by John R. Meyer, Jr., and G. Falabella, Jr.; “Ap- 
proximate Solutions for the Detached Shock Problem,’’ by S. F. Borg; “Influence 
Functions for Computing Lift and Rolling Moment of a Twisted Narrow Delta 
Wing,” by Beverly Beane; “Energy Removal from a Supersonic Turbulent Bound- 
ary Layer,’ by John Joseph Martin; “Relation of Force, Acceleration, and Mo- 
mentum for a Body Whose Mass Is Varying,” by H. W. Sibert; “Further Com- 
ments on ‘Turbulent Boundary-Layer Characteristics at Supersonic Speeds,’ ” by 
R. J. Monaghan; ‘On the Boundary-Layer Equations in Hypersonic Flow and 
Their Approximate Solutions,’ by Lester Lees; “On the Use of the Neutral Point 
as a Stability Parameter,’’ by Judson R. Baron; “On Asymptotic Solutions for the 
Heat Transfer at Varying Wall Temperatures in a Laminar Boundary Layer with 
Hartree’s Velocity Profiles,”” by H. Schuh; “Graphical Solution for “Turbulent 
Boundary Layer in Compressible Fluids,’ ’’ by John Joseph Martin; “Additional 
Values of C(k),” by W. B. Brower and R. H. Lassen; “Comments on “Two-Dimen- 
sional Airfoils at Moderate Hypersonic Velocities,’’’ by L. E. Flanagan, Jr.; “‘On 
the Transition Through a Weak Shock Front,” by Paul Mandl 136-152 


Published by the 
INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 











JOURNAL OF THE 
AERONAUTICAL SCIENCES 





Editor: Hucu L. Drypen, Director, N.A.C.A. : 
Associate Editor: Rosert R. Dexter Managing Editor: Berneice H, Jarck 





Editorial Committee 


Structures and Materials Aircraft Design 


W. E. Beall 
W. K. Ebel 
Hall L. Hibbard 


Hans Reissner 
Abe Silverstein 
C. Fayette Taylor 
E. S. Taylor 
BS Tho 

. S. Thompson 
Electronics in Aviation H. S. Tsien 
K. C. Black 
E. F. Herzog 
Devi S. Little 

D. K. Martin 

George Rappaport 
Paul Rosenberg 
L. M, Sherer 


Instruments 


Allan G. Binnie 
J. R. Bradburn 
yr ge H. Colvin 


Eskin 
Esval 
Physiologic Problems W. Holzapfel 
L. Howland 
H. G. Armstrong H. Purcell 


Savage 

ugo Schuck 
ard, Jr. 
Statham 
Sylvander 
K, Warner 


Robert J. Benford 

Eugene Du Bois 

W. Randolph Lovelace, IT 
Ross A. MacFarland 


C. §. Dra 
S. G. 
O. E. 
H. 
W. 
Louis H. Bauer rg 
0. 
CL 
L. D. 
x, <; 
D. 


Vibration and Flutter 


Lee Arnold 


M. A. Biot Fuels and Oils 


R. L. Bisplinghoff 
William Bollay 
Chieh-Chien Chang 
I, E. Garrick 


D. P. Barnard 
J. H. Doolittle 
Graham Edgar 

. T. Goodwin 


Meteorology 


C, F. Brooks 
D, M. Little 
E, J. Minser 
F, W. Reichelderfer 


R 
S. D. Heron 
A 


. M. Rothrock H, C. Willett 


Martin Goland 
J. P. Den Hartog 





SUBSCRIPTION RATES 


jo ee or THs Azronauticat Sciences, Published Monthly.—United States and Possessions: 1 Year, $12.00; Single — 
1.50. Foreign Countries Including Canada (American Currency Rates): 1 Year, $13.00; Single Copies, $1.50, f 
AgronauticaL Encinesaine Review, Published Monthly.—United States and Possessions: 1 ot $3.00; Single Copies, $0.50. xt 
Foreign Countries Including Canada (American Currency Rates): 1 Year, $3.50; Single Copies, $0 . ; 
Notices of change of address should be sent to the Institute at least 30 days prior to actual change i abivain: 

Manuscripts for publication, proofs, and all correspondence should be addressed to the Editorial Office of the Jourwat. 


Correspondence regarding subscriptions may be addressed to Publication Offi : “a bev and Northampton Sts., Easton, Pa., or to the. By 


Editorial Office of the Tebenas or THE AgronautTica. Sciences, 2 East 64th Street, New York 21, New York. 


Copyright, 1953, by the Institute of the Aeronautical Sciences, Ine. 


Entered as Second Class Matter at the Post Office, Easton, Pa., May 1, 1937. Acceptance for mailing at a special rate of 
Postage provided for in the Act of August "24, 1912. ‘Authorised April 29, 1987 : 











creat 
coeff 
certa 


vente 





JOURNAL OF THE 
AERONAUTICAL SCIENCES 





—_-— 


VOLUME 20 


FEBRUARY, 1953 


NUMBER 2 





——— 


On the Motion of a Tumbling Body 


A. M. O. SMITH* 
Douglas Aircraft Company, Inc. 


SUMMARY 
© The Douglas D-558 and D-558-2 research airplanes incorpo 
"rated releasable nose sections to facilitate pilot escape in the case 
' of an emergency. Because the free-nose sections did not have 
static stability, there was fear that the noses would fall with a 
steady tumbling type of motion. 

Accordingly, an extensive investigation of the motion of the 
Moses was initiated. In one phase of the work, small models 
dynamically similar to the full-scale noses were released in the 
PNA.C.A. 20-ft. Spin Tunnel to study free-fall characteristics. 
hin these low Reynolds Number tests the noses did tumble. 
» Later, in the high Reynolds Number tests on dimensionally 
similar models that were mounted to spin freely on a shaft, the 
'models could not be made to autorotate. To explain the dis- 
crepancy, a modest investigation of the mechanics of tumbling 
was begun. 

It was learned that the motion can be exceedingly complex. 
Factors that definitely produce a continuous net driving force 
per revolution are: (1) hysteresis effects on the lift coefficient, 

) created by rotation; (2) Reynolds Number effects on the drag 
coefficient; (3) static net rotating moments characteristic of 

| certain body shapes such as the anemometer; and (4) a coupling 

) effect that enables rotational energy to be gained by an otherwise 
neutral body when allowed to tumble freely in space. 

In addition, the following phenomena considerably influence 
the rotation and possibly contribute to the driving moment: 
(1) the effect of the shed vortex field on the moments of a rotating 
body; (2) spinning as with wings in the region, where the C;-vs.-a 

| curve has a negative slope; and (3) the Flettner rotor effect, 

Which may sometimes occur if the body spins about a longitudinal 
axis. 

In the light of the considerations presented herein, it can be 
Seen that tumbling may be greatly influenced by Reynolds 
Number. In particular, it appears that bodies similar to the 
D-558 noses will tumble at a low rotational rate. However, 
Means are pointed out by which tumbling can probably be pre- 
vented if so desired 


(1) NOMENCLATURE 


Cp = drag coefficient 

Co, drag coefficient based on frontal area 

Received April 7, 1952. Revised and received September 5, 
1952. 


* Supervisor, Design Research, El Segundo Plant 


Co, drag coefficient averaged over a revolution, 
Co, = z 
a7 
= lift coefficient 
lift coefficient based on frontal area 
= moment coefficient 


6 average of Cm over one revolution, 


1 wor 
Cm = 9 I Cu de 
aT 0 


moment coefficient based on frontal area and length 


of body 
6 average of Cn, over one revolution, 
1 2n ; 
Gm, = Cm, dé 
a7 0 
= average incremental moment coefficient 
acceleration due to gravity 
= principal moment of inertia of the body with respect 
to axis perpendicular to plane of rotation. 
= length of body 
Mach Number; also moment 
dynamic pressure 
Reynolds Number based on length 
Reynolds Number based on diameter 
reference area 
time 
velocity 
weight 
= displacement in cartesian coordinates 
angle of attack, for airfoils 
angle of attack, for bodies 
mass density of air 
instantaneous angular velocity, rad. per sec 
= angular velocity averaged over a revolution 


(IL) INTRODUCTION 


AS PART OF THE DouGLas D-558 high-speed airplane 
research project, an extensive investigation of the 
motion of a jettisonable nose was performed for the 


Navy Bureau of Aeronautics.t The work was for the 


purpose of determining the loads imposed by the nose on 


+ Submitted April 21, 1950, as Douglas Report ES 15379, 


Parts I through V. The present article is a revision of Part V. 
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Fic. 1. 
D-558-2 airplanes showing the jettisonable nose sections 





The D-558-2 nose in the Southern California Coopera 
tive Wind Tunnel. 


Fic. 2. 


a pilot riding within. Fig. 1, consisting of photographs 
of two balsa models dynamically similar to the D-558 
and D-558-2 airplanes, illustrates the general appear- 
ance of the releasable noses. 

Originally, in the research to determine loads im- 
posed on the pilot by release of these noses, it was 
assumed that the tumble. N.A.C.A. 
Spin Tunnel Tests on a free model! and '/s-scale 
model tests in the Douglas El Segundo Wind Tunnel 
on a shaft-mounted model’ substantiated the belief, 


noses would 


for both types of tests did produce tumbling motion. 
However, the '/s-scale models (Fig. 2) built for testing 
in the Southern California Cooperative Wind Tunnel 
would not autorotate about the fixed shaft even though 
bearing friction was low. Except for inevitable dif- 
ferences in the models due to construction inaccuracies, 
the main difference between the ' ;)-scale models and 
the '/s-scale models was size that is, Reynolds 
Number. However, the instrumentation 
on the large D-558-2 model but not on the small one. 


boom was 


This boom is illustrated in Figs. 2 and 15. 
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Because of the unexpected discovery that the larg 
model would not autorotate, studies were undertaker 
in order to learn whether rotation about a fixed ayjs 
might impose such restriction that the model woul 
not autorotate and whether it would autorotate j 
In the 
process of the investigation several points of interest 
Since the literature on this 


allowed three or more degrees of freedom. 


were brought to light. 
subject is scarce, the findings are described hereiy 
to explain the results of the D-558 and D-558-2 nose 
rotation tests and to assist others in treating the tum 
bling problem. 

In this paper tumbling will be regarded as a free. 
body rotation, predominantly about an axis _perper- 
dicular to the air stream but not exclusively so, fora 
often be discerned 


combination of motions 


Rotation about a fixed axis will not be called tumbling 


may 


because the axis is fixed, but instead it will be des. 
ignated as autorotation. 


(III) GENERAL NATURE OF TUMBLING MOTION 


Consider the free body of Fig. 3, treated as if it wer 
m equilibrium in a vertical wind tunnel, in order ti 
eliminate considerations of translation as in free air 
Because the body undergoes lift and drag fluctuations, 
its venter of gravity obviously will travel in a closed 
path or orbit. Since the velocity of travel along this 
path is small, it can be assumed with only small error 
that the vector wind velocity over the body is not 
affected by the translatory motion and _ direction 
Furthermore, it is assumed that the motion is predom- 
nantly inone plane. This plane of tumbling may gradu- 
ally rotate because there are no forces to hold it fixed 
However, such drift has only a minor effect on the de- 
velopment of the motion. N.A.C.A. motion-picture 
studies in the 20-Ft. Spin Tunnel indicate that fully 
developed tumbling of the D-558 and D-558-2 noses 1s 
substantially a motion in a vertical plane. Often the 
tumbling state is entered by means of a three-dimer- 
sional type of motion to be explained later. However 
since the steady-state tumbling is a coplanar motion, 
explanations based on coplanar analyses should be 
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ON THE MOTION 

Because the axis of rotation is assumed not to tilt 
and because there is negligible 2 motion, the equations 
of motion can be written as follows if it be understood 
that C,, Cp, and C,, represent dynamic coefficients: 


Lift: 
W d*x l pa eae 
g dt? =5eC, SV? (1) 
Drag: 
W d*y I oe : 
2 di? = 9G, 3" = a (2) 
Moment: 
d*6 l een 
dt? = 5 p 4 m S11 # (3) 


Because perturbation velocities for a body will usually 
be small, these three equations can be solved suc- 
cessively instead of simultaneously. 

To solve them, the 6-vs.-time relation must first be 
found. 


Letw = d6/dt. Then 


d*6 ; dw 


dt? dt 


dw d6 2 dw 
do dt 7 dé 


and Eq. (3) can be integrated by means of the above 


substitution. 


w” _ pSlV? of = wo” 
tae ee a 


or 
do lpSlIV? {7% 
=w= ¥ C,, a0 9" (4) 
dt \ I | sles 
where w is the angular velocity when 6 = 0. On 


integrating again, 


i de 
0 pSl V? 


78 
Cm d0 9? 
y ] / dé + w 


Now with the time-@ relation established, x-@ and y-6 


t= aa to (5) 


relations are determinable, since C, and Cp are assumed 
to be known as functions of 8. 

Consider Eq. (4). It shows that the resultant rota- 
tion is a function of the C,-vs.-@ curve. (However, 
in general, this curve will vary with the angular velocity 
because of tumbling 


dynainic effects.) In steady 


/ C0 = 6 
0 


but for any value of @ except 27, 47, etc., the value 
Therefore, a steady- 


or autorotation, 


of the integral is not equal to zero. 
State rotation is one having a reference velocity «w» 
upon which is superimposed a fluctuating one. The 
amplitude of the fluctuations can be considerably 


OF 


“J 


A TURBLING BODY 




















—_ AERODY NAMIC 
Cm 
DRIVING MOMENT 
+ 
w 
——— Vv 
BEARING 
Fic. 4. Schematic breakddwn of rotating forces acting on a 


rotating body 


changed, as is seen from the fact that the variable term 


SIV? f 
L. 10 
iJ 


For instance, note that increases in the moment of 


in Eq. (4) is 


inertia will tend to reduce the fluctuations. 

A matter of great interest in considerations of tum- 
bling motion is the effect of translational speed on the 
steady-state rate of rotation. Eq. (4) can hardly 
answer the problem because all it basically shows is 
that i” C,, d@ must be zero if the rotation is tobesteady 

that is, periodic. Therefore, the prediction of steady- 
state rotational rates reduces to the problem of pre- 
dicting the value of De C,, d6as a function of w/ V, V, p, 
and M. 

Since tumbling normally involves airflow separation, 
predictions of C,, as a function of @ are not likely to be 
accurate. However, if Reynolds Number effects and 
Mach Number effects are negligible, it appears that 
the steady-state value of w is nearly proportional to V, 
regardless of density as will be shown. 

The fundamental moments acting on a body are as 
shown schematically in Fig. 4, where mean values for a 
revolution are used. The case illustrated is that of a 
body mounted on a shaft. 
driving moment as for a windmill or anemometer. 
In addition, greater driving moments may develop when 
Hysteresis in C,-vs.-a@ curves is an 
If the rotation is about a 


There may be a net static 


rotation occurs. 
example of such a force. 
fixed shaft, there is a bearing-friction damping moment 
of relatively constant magnitude. Also, there may be 
damping caused by lift and damping caused by drag. 
The damping caused by lift is the type treated in non- 
steady airfoil theory. Steady-state rotation occurs 
when mean driving moments just equal mean damping 
If the slope of the net moment curve is 
internal 


moments. 
negative, the rotation will be stable. An 
damper within the body is possible because fluctuating 
velocities exist from which energy can be absorbed. 
A piston in oil in a cylinder could absorb such energy. 

For a given body, the net moment coefficient appears 
to be a function of w/V. Reference 3 shows by poten- 
tial-flow analysis that damping of an oscillating wing 
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TYPICAL POSSIBLE END VIEWS 


SIDE VIEW 


Typical bodies having planes of symmetry through the 
axis of rotation. 


lic. 6. 











Cm 
+ 
lies. 
/ v — v o 
Fic. 7. Pitching moment typical of the bodies of Fig. 6 
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Fic. 8. Simple anemometer. 


is a function of w/|’. By inference, a potential-flow 
analysis would reveal the same parameter for a rotating 
airfoil. This type of damping is derived from lift 


effects only. Damping created by drag effects is 


revealed by the following illustration: Consider the 
resisting torque of a circular cylinder of radius ry 
rotating about a transverse axis in a wind, as in Fig. 5. 

According to well-known flow-resolution theorems, 
the effective relative velocity over the element dA,, 
is V.cos@ — wr. The drag of this element is 


dD, = 1/2p Cp dA,(V cos 6 — wr)? 
Its torque is 


dM, = 


1/2p Cp dA, r(V cos 8 — wr)? 


FEBRUARY, 19538 
The effective relative velocity at the same radiys 
on the upstream end is V cos 6 + wr. Then the ne 


torque developed by these two elements is 
dM = dM, — dM, = 2pr (wr) dACpV cos 6 (6 


In coefficient form, 


. dM r dA. 
an = i = 4 Cp cos 6 
l 2p I 2A Yo Yo A ( 


= 
: (7 
V/ 


Again w |’ steps forth as a parameter, determining 


the damping-moment coefficient. In the above illustra- 
tion, Cp has been assumed to be independent of velocity, 

Hysteresis in C,-vs.-a curves constitutes one of the 
most powerful driving forces. When an airfoil, for 
example, is nearly at its steady-state stall, the boundary 
layer is thick as a result of flowing against adverse 
pressure gradients. However, if an airfoil is instantly 
increased from zero lift to a high angle of attack, 
the boundary layer at first remains as it was for zero 
lift. Consequently, no separation occurs immediately, 
even though the usual stall angle is far exceeded. 
Thus it is apparent that for a rotating body the condi- 
tion of the boundary layer at any angle of attack isa 
function of the distance the fluid has traveled while the 
body is reaching that angle. The distance traveled by 
the fluid while a body changes its angular position by 
an amount aisx = Va/w. 

Therefore, it is probable that steady-state tumbling 
rates will usually be proportional to velocity regardless 
of air density. However, such a relation can be ex- 
pected to break down if Reynolds Number effects or 
Mach Number effects Even 
without them, w/ V may not remain constant. Eq. (4 
shows that the fluctuating velocity is affected by p and 
I. With very low values of /, for instance, the fluctua- 
If J is very large, 
In view of 


become appreciable. 


tions in velocity may be very high. 
the rotation rate will be nearly constant. 
the prior discussion of effects of angular velocity on the 
that either 
moment-of-inertia changes or density changes call 


various moments, it becomes obvious 
modify the net C,,-vs.-0 curve of a tumbling body and, 


consequently the average angular velocity. 


(IV) AUTOROTATIONS EXPLAINABLE BY STATIC 
EFFECTS 


Continuous autorotation of a body symmetrical about 
a plane through the axis of rotation, as in Fig. 6, cannot 
be explained by static data alone. This point is met- 
tioned because literature on the tumbling of wings 
frequently implies that tumbling will develop if the wing 
becomes longitudinally unstable. This condition alone 
is insufficient, for instability merely gives the body 
an initial tendency to pitch away from normal flight 
attitudes. 

Consider the bodies of Fig. 6 rotating about a point! 


on the plane of symmetry. With the wind blowing 4s 
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ON THE 
shown, the center of pressure is ahead of the point of 
rotation, and a turning moment will develop with the 
The pitching mo- 


dightest displacement from 6 = 0. 
If this indeed is the 


ment will be somewhat as in Fig. 7. 
moment curve, net work proportional to f{* C,, dé 
issupplied by the rotation to 180°. However, the work 
derived from.a clockwise rotation will be exactly equal 
and opposite to the work derived from a counter- 
clockwise rotation. ‘Therefore, no net work can be 
computed for such a body if static-moment-coefficient 
data alone are considered. Fore-and-aft shifts of the 
point of rotation greatly change the shape of the C,,,-vs.- 
fcurve, but the antisymmetry will remain as long as 
the axis of rotation is on the plane of symmetry of the 


body. 


The Anemometer 


Nevertheless, there are cases of autorotation that can 
be explained by static data alone. One is the anemo 
meter (Fig. 8). It utilizes the principle that a body has 
a lower drag coefficient when heading in one direction 
than when heading in the opposite direction. In Fig. 8, 
cup A has a higher drag coefficient than cup B.  There- 
fore, a net torque is developed which is able to produce 
rotation. 

It is also apparent from cousideration of Fig. 8 that 
cup B could be removed and that rotation would still 
surely occur if the moment of inertia of the system 
were not too low. Furthermore, this type of torque 
differential between positions A and B will continue 
to exist for a single cup as the arm AC is shortened, 
until the centerline of the cup is on the very point of 
rotation C. When, finally, the axis of rotation is on 
the centerline, according to static considerations, the 
cup can produce no work because of symmetry. 

The above discussion of the anemometer should make 
it clear that net torques can indeed be developed by 
symmetrical bodies of the type shown in Fig. 6. Such 
bodies will develop net driving moments if the axis 
of rotation is moved off the plane of symmetry as to 
0. This type of driving torque will be called the ‘‘ane 


mometer-cup effect.”’ 
(V) Dynamic EFFECTS INFLUENCING AUTOROTATION 


The Lanchester Tourbillon 

The Lanchester tourbillon is a fascinating toy wind- 
mill-like device that will be described here for the 
sake of completeness. Its motion is not a tumbling 
type, for its axis of rotation is parallel to the air stream 
like that of an ordinary windmill. 

When given an initial spin in either direction, it will 
continue to rotate in that direction. When pointed 
directly into the wind, no static lift is possible because 
of symmetry (Fig. 10A). However, as soon as rotation 
occurs, the relative wind over a typical element A, 
Fig. 9 becomes as in Fig. 10B or 10C. 
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Fic. 10. Lanchester tourbillon at three rotational speeds 
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Fic. 11. Motion of a strip of paper. 

Lift and drag vectors L4 and D, are shown. If a 
is measured as shown, lift increases as a is reduced 
that is, the slope of the lift curve is negative. When 
the lift and drag vectors are of such magnitude and 
direction that their resultant has a component in the 
direction of travel of the element, a driving force is 
developed. Fig. 10C shows the vector diagram for 
element A at a rotation rate above the autorotation 
limit. The lift and drag vectors are turned so greatly 
by the change in relative wind direction that the 
resultant force is in a direction opposing the motion. 
Consequently, autorotation of this nature can be pre- 
dicted, given the necessary lift and drag data. From 
Fig. 10, it is apparent that R, provides a force in the 
direction of motion if L4/D, > V4/Vw. A flat spin 
is a rotation of this nature. 


A Falling Strip of Paper 


A flat strip of paper, having an aspect ratio of two 
or more will fall in about a 1:1 glide ratio, rotating 
as in Fig. 11. The motion can be explained very 
simply in the light of Fig. 11, which shows successive 


stages of the tumbling motion. 
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Fic. 12. Plan view of rotating sphere test. 


Consider the forces on the strip when it is at position 
2, shown to larger scale at the left in Fig. 11. It is at 
considerable angle of attack, exerting values of lift and 
drag that produce a resultant force R. If the center of 
gravity is behind the line of action of this resultant (25 
per cent M.A.C. at low angles), a stalling couple is pro- 
duced, tending to increase the angle. Rotation is 
rapid. Consequently, the boundary layer is thinner 
than for the steady state and the stall is delayed to 
higher values of C;. At position 3, the strip probably is 
completely stalled. At position 4, it continues to be 
stalled. However, between positions 4 and 5, adhering 
flow over the trailing edge will be established. Re- 
establishment of unsaturated flow is delayed to angles 
below the static stall angle because the thick wake is a 
Hence, the downward lift impulse 
Moreover, the 


retarding influence. 
is less than the upward lift impulse. 
clockwise-rotation impulse is greater than the counter- 
clockwise-rotation impulse, and the inertia of the paper 
carries it past the region of retarding impulse. The 
motion from’5 to 9 is essentially a repetition of that 
from | to 5. It can be concluded that this motion is 
characterized by the following features: 

(1) The primary motive source is a hysteresis in the 
C_-Vs.-a curve. 

(2) Asa result a net upward lift is always developed, 
producing a finite gliding angle. 

(3) Probably the center of gravity must be between 
Otherwise the net lift will 
See the detail diagram for 


25 and 75 per cent M.A.C. 
produce a damping couple. 
position 2, Fig. 11. 
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(4) Damping moments appear to be a mixture of the 
lift and drag types. 

Tumbling of tailless airplanes is a motion of this 
nature. Note that the fluctuating lift and drag forces 
will cause the center of gravity to follow a wavy path at 
a nonuniform velocity. 


A Rotating Dumbbell 


Autorotation of two 5-in. diameter balsa spheres 
(Fig. 12) was demonstrated in the Douglas E1 Segundo 
Wind Tunnel of 30- by 45-in. cross section. 
sembly rotated on miniature ball bearings having ex. 
tremely low friction. The rotation is explained as 
follows: The drag coefficient of a sphere is a function of 
Reynolds Number only, assuming the Mach Number t 
At the sphere critical Reynolds Number, the 


The as- 


be low. 
drag coefficient decreases rapidly. 
upstream will then be at a higher Reynolds Number 
Hence, if the drag 


The sphere moving 


than the one moving downstream. 
coefficient decreases rapidly enough, the sphere moving 
upstream will have less drag than the sphere moving 
downstream. A criterion for such tendencies can be 
derived. 

For one of the spheres, D = (1/2) pSV*Cp, where Cp= 
For autorotation, drag must decrease as |" in- 


creases or 


iD 1 oCp 
om ps|2vr + 52] <0 


dV 2 
_ sl2vc,+ V2 oe | ) P 
ae a ‘sii OR, dV 
Since R, = Vd/v, 
(d/dV)R, = d/v = R/V 
ais sv [2c +R oe 9 
7s hUhE CC 
For neutral stability 
2Cp + Ri(OCp OR,) = (0) 
or for autorotation 
—Ra(OCp/ORa > 2Cp (10 


A typical plot of Ry(OCp/ORa) is shown in Fig. 15 
There is a considerable range of Reynolds Number over 
which autorotation might be expected. The original 
test observations as written by J. G. Brazier are 
quoted verbatim: 

“The tunnel fan was started, and from dynamic pres- 
sures of zero inches of water to three and one-half inches 
of water (Rz = 296,000) the spheres would not autoro- 
tate. Any forced rotation was quickly damped out. 
At the lower dynamic pressures the spheres assumed 4 
tandem attitude and oscillated with an amplitude o/ 
+30°. Either sphere would remain in the upstream 
position. Between dynamic pressures of three and one- 


half and four inches of water the spheres were close t 
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ON THE MOTION 
being neutrally stable. After about one complete revo- 
lution either way, the rotation would damp out. At 
these dynamic pressures the spheres assumed a side-by- 
side attitude and oscillated with amplitudes of + 15°. 
At a dynamic pressure of four inches of water (Ry = 
315,000), the spheres began to autorotate, the autorota- 
tion could be maintained in either direction. The mo- 
tion was self-starting due to the initial oscillatory mo- 
tion of the spheres. The rotational speed was ap- 
proximately three cycles per second. The maxi- 
mum dynamic pressure reached was approximately 
five inches of water (Ry = 353,000). The mounting of 
the spheres to the hub and the hub to the shaft was not 
rigid. After about one-half hour of running time the 
spheres were lost due to the failure of the vertical shaft 
just below the hub. 


“Unfortunately, due to the loss of the model, the 


upper limit of the autorotation was not determined.”’ 


Consequently, a complete check of the theory was 
never made. However, the correlation of the start of the 
autorotation with sphere critical Reynolds Number was 
so good that there appears to be little doubt about the 
explanation of this type of autorotation. Sphere critical 
Reynolds number for this tunnel is about 375,000. 


Shed Vortex Effects on Tumbling Motion 


A tumbling body, especially a tumbling wing, de- 
velops lift of a cyclic nature. Therefore it will shed 
vortices of alternate right-hand and left-hand rotation. 
It is believed worth while to mention here that there 
may be a possibility of deriviiig a driving force from the 
downwash field created, depending on the value of 
w V. However, reference 3 indicates that an airfoil 
oscillating harmonically about its midpoint is affected 
not by the wake but only by the motions and accelera- 
tions with respect to the free stream. Therefore it is 
improbable that the wake will have any effect on an air 
foil rotating at constant speed. Nevertheless, the ques 
tion merits further investigation. 


Effects of Change from Shaft Rotation to Free Tumbling 


If the restraint to rotate about a fixed axis is elim- 
inated, a completely free rotation occurs. On the shaft, 
only one degree of freedom is possible. But when free, 
asin a vertical spin tunnel or falling through the atmos- 
phere, six degrees of freedom are available. In most 
cases observed, free tumbling appears to remain a mo- 
tion in one plane, although the orientation of the plane 
may gradually change. Free motion in a plane then 
differs from rotation about a shaft, in that the center of 
gravity is free to travel up and down and back and 
forth. 


Therefore a major question arises. Is it possible to 
obtain some sort of coupling or phase relationship dur- 
ing free rotation which will generate a driving moment 


that does not exist during pure shaft rotation? 


OF 
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The general problem is to find motions that make 


2) 79 


| C, d6 > | C,, d0 


Free Shaft 
Rotation . Rotation 
In the search for a driving mechanism, detailed knowl 
edge of the free motion is desirable. A good starting 
point is to make calculations of the motion based on a 
C,,-9 curve antisymmetrical about the line 6 = m—that 
is, one having 


/ C,, dd = 0 
Pe Al] 


Such can be supplied by either static or dynamic tests 
of a model supported on a shaft. If 


* 


Jt C.. @= 0 
0 


the model will be rotating at a steady averagerate. Then 
the motion can be examined to learn whether there is a 
net perturbation moment and, if so, what its sign is. 
If no perturbation moment is developed, the initial 
motion, of course, will continue unchanged because 


| C,, dd = 0 
0 


Fig. 14 shows the lift, drag, and moment data used in 
calculating several free-motion orbits. Three c.g. posi 
tions were studied to show the effect of static stability 
on the motion. The curves in this figure represent 
Dynamic values for lift and drag could 
ap- 


static data. 
have been used without trouble; 
peared to be a needless refinement, for differences be 
For 


but such use 
tween the static and dynamic values are small. 

purposes of convenience, moment-coefficient values for 
the center of gravity at 0.676/ were approximated by 


sinusoidal functions as follows: 


0.15 sin 26, —90° < 6 < 90 


270 


Cm, 676 ~ 
= —().075 sin 46, 90° < 6 


The values of C,,_ for the other c.g. positions were 


derived by conventional moment-transfer methods. 
The values can be considered consistent with the test 
results of Runs 55 and 56 of reference 5. 

In making the calculations, the method of Section 


(III) was followed, using Eq. (4) rewritten as 


do 


SIV? i wo” 
= yy f/ c.0+—, (4a) 
dt / Jo pSlV? I 
In addition to the three c.g. positions, the calculations 
covered five rotation rates: 


w 
. r r co) 


= 0.008125, 0.0325, 0.13, 0.52, 2.08 
V pSIV? 1 


Note that the values chosen increase in geometric pro 
gression. The calculations were made in this order: 
The 6-time relation was found; from it, the Cp-time re 
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Fic. 14. Representative static lift, drag, and moment coefficients for three c.g. positions 
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Fic. 15. Model D-558-2 nose. Calculated effect of position of center of gravity and tumbling rate on path of center of gravity 


while tumbling in a vertical wind tunnel. Conditions represent the full-scale nose falling at terminal velocity at 15,000 ft. Com 
puted from data of Fig. 14 corresponding to R = 8.47 X 105. 
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ON THE MOTION OF A TUMBLING BODY 81 
= TABLE 1 
Results of Orbit Studies at 15,000 ft.; D-558-2 Nose; Wt. = 775.5 Lbs. 
Average Average 
Angular Angular Velocity, Angular Falling 
Rotation Drag Coefficient Rad./sec. Velocity Velocity, Average 
C.G. Position Parameter Cox at @ = 0° Ft./Sec wl /V 

0.501 0.138 Rotation impossible 

0.52 1.06 2.26 4.12 254 0.074 

2.08 0.996 17.03 17.27 266 0.533 
0.6761 0.008125 0.586 1.73 0.101 342 0.042 

0.0325 0.746 2.39 0.359 303 0.066 

0.13 0.868 3.07 1.33 281 0.091 

0.32 0.947 5.94 5.10 269 0.184 

2.08 0.961 20.48 20. 24 267 0.6388 
0 880/ 0.18 0.817 3.92 1.26 290 0.118 

0.52 0.924 6.43 4.74 272 0.197 

2.08 0.960 17.64 18.63 267 0.550 

a TABLE 2 
: — 1 2m > 
Values of ACnyz = 5 f, ACn_ dO 
ao 0 

Rotation Parameter 0.008125 0.0325 0.13 0.52 2 08 
C.G. position 
0.5 Will not rotate —().0021 — 0.0007 
0.6761 0.0061 0.0045 0.0033 0.0014 0.0004 
0.80/ 0.0059 0.0031 0.0009 


lation was found; and from them, the steady-state 
falling velocity was calculated. Computations apply to 
the full-scale D-558-2 nose at 15,000 ft. Because varia- 
tions in ¢.g. position and ap Vv pSIV? ‘T both affect the 
(p-time relationship, the steady-state falling velocity is 
not constant for the various cases. In other words 
Co the mean drag coefficient varies. 

Results of these calculations and the orbits of the 
center of gravity are shown in Table | and Fig. 15, re- 
spectively. At the lowest rate of rotation the orbits 
have considerable size; 
rapidly approach shaft rotation in nature. 
the arrows show the orientation of the nose as it moves 
It is in its highest position when 
that is, its drag is 


as the rate increases, they 
In Fig. 15 


around its orbit. 
pointing straight down when 
minimum. 

The arrows are located to intersect the orbits at 0.50/, 
0.676/, and 0.80/ to provide better clues as to the posi- 
tion of the nose during its gyrations. N.A.C.A. test re- 
sults at low Reynolds Number are tabulated in the 
figure. 

The motions shown are based on static data in 
which 


oO) 


/ C2 =9 
0 


[hey have been computed upon the assumption that 
the translational velocities are so low that they produce 
a negligible effect on either the direction or magnitude 
of the true resultant velocity. The data of Table 1 
were used as a first approximation in an iteration proce- 
dure for the purpose of learning the effect of using the 
true velocities and directions. 

In the second approximation, the true velocity rela- 
tive to the nose was computed by vectorial addition of 


the oribital velocity to the free-stream velocity. Then 
for the corrected angle of attack and velocity, an incre- 
mental moment coefficient AC», at each @ position was 
obtained from the moment curves of Fig. 14. If the 
summation of the incremental values were positive, a 
incremental driving moment existed; if negative, a 
damping moment existed. The results of this study for 
each of the ten cases treated are presented in Table 2. 

For the most stable case (c.g. at x// = 0.50), ACn, is 
negative. ACm_ 1s positive for the other two c.g. posi- 
tions, although all the values are small. For example, at 
15,000 ft., at the approximate terminal velocity of 186 
m.p.h. with the center of gravity at 0.676/, the average 
angular acceleration @ is only 0.026 revolutions per sec. 
per sec., using ACn, = 0.0033. The analyses also 
showed that the changes in direction and magnitude of 
the relative-velocity vector were so slight that the 
original assumptions of independence of the x, y, and @ 
motions were exceedingly good. 

The findings recorded in Table 2 are of considerable 
significance. They show that the freedom of translation by 
itself can induce net moments, either driving or damping. 
Therefore. tumbling may at times be possible even 
though static or dynamic data from shaft-supported 
models indicate otherwise. 

An approximate physical picture of the origin of these 
dynamic driving moments can be given. Consider a 
model tumbling in equilibrium in a vertical wind tunnel. 
In the falling part of its orbit, its relative wind velocity 
is higher than the average value. If, during this 
period, the moments are in the direction of rotation, 
they are increased because of the additional velocity. 
In the other half of the orbit, the model is rising and 
the relative wind velocity is below the mean. In this 
phase, moments are resisting the rotation. However, 
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Fic. 17. Transition of a nose from rocking motion to tumbling. 
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Fic. 18. Effect of &@ on net damping moment !/¢-scale model 
results. 


their resistance is reduced by the reduced relative 
wind. Therefore, a net driving moment for an entire 
cycle is created. A careful correlation of Figs. 14 and 
15 will make these relations clear. In effect, the pertur- 
bation velocities distort the C,,_-6 curve. 


Other Possible Motions 


Since free-tumbling motion has six degrees of free- 
dom, there are undoubtedly numerous other possible 
motions of complex nature. For example, an elongated 
body tumbling end over and may also be subject to aero- 
dynamic forces tending to rotate it about its longitu- 
dinal axis as in Fig. 16. 

There are certain combinations of shape and Reynolds 
Number that will provide continuous rotation about the 
longitudinal axis. With such, a Flettner rotor effect 
may be developed which will give a turning moment 
about an axis parallel to the relative wind and through 
the center of gravity. Thus a rotation about all three 
axes may result. The general system of the several 
possible impressed aerodynamic forces can induce com- 


plex motions indeed, particularly when SZYTOscop} 
effects are added. 

Motion pictures taken during the studies of refereng, 
1 showed an interesting motion during the transitio, 
from rocking to steady-state tumbling. A nose mod 
was released broadside to the airflow in the spin tunn¢ 
(Fig. 17). At first it floated around with a rocking 
motion. Gradually the amplitude built up. If ty 
model had been supported on a horizontal shaft, th 
amplitude of rocking might possibly have increased unt 
it went over “top dead center”’ and into continuous p 
tation. Instead, when the nose was within about 20) oj 
pointing straight up, it appeared to fall to the side an 
keep on, instead of swinging back. In short, entry inty 
the tumbling state was not made by motion in a verticy] 
plane. The entry was first a build-up of rocking ampli 
tude in a vertical plane until the first full revolution of 
tumbling was completed by rotation in a plane inclined 
about 45°. After the rotation was well established, the 
plane returned to vertical. One gained the impressio: 
that the model did not have enough kinetic energy to 
swing its weight right over “‘top dead center,”’ so it got 
around the obstacle by tilting its plane of rotation toa 
position where not so much potential energy would be 
required. Asa final steady motion was substantially in 
a yertical plane, the above transitory motion does not 
invalidate the assumptions of three degrees of freedom 
generally used in this report. 


(VI) ANALYSIS OF THE MOTION OF THE D-55S8 AnD 
D-558-2 NOSES 


For elongated bodies roughly similar in shape to the 
D-558 and D-558-2 noses, the following effects appear 
to determine the tumbling motion: (a) the anemometer 
cup effect, (b) the falling-strip-of-paper effect, (c) the 
rotating-dumbbell effect, (d) drag-type damping, and 
(e) the free-rotation effect. 

The anemometer-cup effect exists because the models 
are not symmetrical fore and aft and because the center 
of gravity is not on the axis of symmetry. An 1dea ol 
the magnitude of this effect can be obtained by examin 
ing the static C,,_ values obtained from the tests on the 
D-558 and D-558-2 noses (Figs. 18 and 19). In general 
they are not equal to zero. The values shown on thes 
two figures were obtained by planimeter from the 
moment data of reference + on models mounted 4 
shown in Fig. 2. 

For these bodies the falling-strip-of-paper effec! 
should not be large because they are not primarily lift 
ing bodies. This statement is borne out by examinatiol 
of the C,-vs.-6 data of reference 4, which shows, in 
general, no large increase in dynamic lift coefficient 
However, the effect can definitely be seen. For exatl 
ple, see Fig. 22a of reference 4. 

The rotating-dumbbell effect appears not to be ® 
driving force at large scale. However, it evidently col 
tributes the necessary driving moment for the smal 
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ON THE MOTION 

yscale nose tests of reference 2._ Damping is probably 
nearly all of a drag type, for the reason that the bodies 
are poor lifting shapes. 

The first four effects combine to produce the dynamic 
moment coefficient. Where dynamic test data are 
available, there is no need to consider the factors causing 
the resultant curve. Mean values of Cn, for the two 
noses are Shown as a function of both @ and Mach 
Number in Figs. 18 and 19. Whereas the small '/50- 
sale noses would autorotate, the large noses definitely 
would not. In fact, the damping coefficients (Figs. 18 
and 19) are high. 


With this qualitative analysis of the motion of the two 
noses, we Shall consider whether they are likely to 
tumble in full-scale flight. The known facts concerning 
them are listed below. 


p-558 Nose 

(a) N.A:C.A. !/4-scale model tumbled in a vertical spin 
tunnel. Test Reynolds Number R = V1/»y = 255,000. 

b) Douglas !/;-scale model rotated on fixed shaft below R = 
142,000. Above this Reynolds Number, it would not rotate. 

(c) Cooperative Wind Tunnel '/s-scale model did not rotate 
ona fixed shaft for Reynolds Number between 73,500 and 7,750,- 
00. At the regular Cooperative Wind Tunnel test conditions 
R>2.5 X 10°), damping is high (Figs. 18 and 19). 


D-558-2 Nose 

(a) N.A.C.A. '/,-seale model tumbled in spin tunnel. Test 
Reynolds Number, R = 120,000. Model had no boom or pilot’s 
enclosure. 

(b) Douglas !/;9-scale model rotated on a fixed shaft at all 
tunnel speeds, R < 236,000. Model had no nose boom. 

(c) Cooperative Wind Tunnel '/¢-scale model will not rotate 
on fixed shaft for Reynolds Numbers, R, between 812,000 and 
8,500,000. Damping was extremely high (Figs. 18 and 19). 


Since Reynolds Number effects were a likely explana- 
tion of the nonautorotating characteristics of the large 
models, several attempts were made in the Cooperative 
Wind Tunnel to induce autorotation of the D-558 nose 
by operation at the Reynolds Number of the smaller 
models. the difference in 
model size, the dynamic pressure became extremely low, 
aslow as 0.1 Ib. per sq.ft. Even with Reynolds Number 
the same as for the small models, the nose still refused to 
rotate. The probable reason it did not rotate is that 
bearing friction was then too high. At q = 0.1 Ib per 
sq.ft. and Cn, = 0.1, the average moment for the D-558 


Because of considerable 


nose is only 0.0044 ft.lb. or about the same as 2 Gm. 


acting at 1 ft. Even though the models had excellent 


hearings, their bearing friction was hardly negligible 
compared to 2 Gm.-ft. of moment. 
The above facts seem to indicate that the N.A.C.A. 


and '/59-scale Douglas models rotated because they were 
ina Reynolds Number range favorable to the dumbbell 
type of autorotation. At Ry, > 2.5 X 10%, the dumbbell 
type of driving force does not exist, and the hysteresis 
effect on lift is insufficient to produce rotation about a 
lixed shaft. 
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comparison of static and dynamic data. '/,-Scale model results. 


Typical damping-moment data are presented in Figs. 
1S and 19. The axes of rotation for these data were 


located as follows: 


Model Longitudinal Position Vertical Position 
D-558 0.556/ 0.021/ above centerline 
1D-558-2 0.6761 0.017/ above centerline 


In Fig. 19, where all available data points are plotted, 
the values at / = 0.41 for the D-558-2 nose are notice- 
Checking failed to show any error in 
One other explanation besides erro- 


ably out of line. 
data reduction. 
neous data appears possible. Considera circular cylinder 
rotating about a transverse axis normal to the wind. 
The part moving downwind will be at a lower Mach 
Number than that moving upwind. For speeds either 
well below critical or well above critical, the curve of 
drag coefficient as a function of Mach Number will have 
nearly zero slope, dCp/dM. However, at the critical 
Mach Number the drag coefficient rises rapidly. There- 
fore, the side of the cylinder moving downstream will 
have a much reduced drag coefficient and that moving 
upstream a much increased value. Therefore, resistance 
to rotation in the critical Mach Number range will be un- 
usually high. Because of the long boom on the D-558- 
2 nose, differential Mach Numbers created by rotation 


will be appreciable. Now the critical Mach Number of 
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analysis shows that tumbling can probably be avoided 


(1) Designing a neutral body shape whose net static 


moment, Cn,» is zero—that is, avoid the anemometer. 


cup effect. 

(2) Operating above the sphere or cylinder critica) 
Reynolds Number 
bell effect. 


that is, avoid the rotating-dumb. 
(3) Properly locating the center of gravity fore and 
aft so that the incremental driving moment due to free 


rotation is zero or negative. 
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D-558-2 nose at approximately full-scale conditions. 
a circular cylinder is just above 0.4. The test results 
showing high damping, exactly at the critical Mach 
Number of a cylinder, form a noteworthy coincidence, 
to say the least. This effect is just the reverse of the 
dumbbell effect developed at the critical Reynolds Num- 
ber of a sphere or cylinder. 

The remaining question is whether free rotation can 
develop sufficient driving moments to overcome the 
damping moments due to rotation. Since the data of 
Figs. 18 and 19 are at nearly full-scale conditions, they 
should account for all effects except that of free rotation. 
The effects of free rotation can best be analyzed by means 
of a plot of damping and driving moments versus @// V, 
as illustrated in Fig. 20 for the D-558-2 nose. If the 
incremental driving moment is greater than the damp- 
ing moment, tumbling will occur. The plotted data 
show a distinct possibility of the occurrence of low- 
speed tumbling, even though the static net damping 
moment Cn, is zero and the dumbbell effect on Cn, has 


disappeared. 


(VII) CONCLUSION 


It is hoped that this article has clarified the nature of 
tumbling motions and has pointed out methods of pre- 


dicting the occurrence of tumbling. In particular, this 
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Slender-Body Theory—Review and Extension’ 


MAC C. 


ADAMS? ann W. R. SEARS} 


Cornell University 


ABSTRACT 


The approximate theory of flow about slender bodies and wings 
originated by Munk and Jones is reviewed. It is presented here 
ina form that emphasizes the relation to the source-sink methods 
of von Karman and others. The extension to noncircular bodies 
is made for subsonic flow, paralleling Ward’s extension for super- 
nic flow. The calculation of pressures and forces and the ex- 
tension of the theory to unsteady flows are reviewed, and some 
discrepancies in the published literature are explained. 

Finally, interpreting the Jones slender-wing result as the first 
term of an expansion in powers of a breadth parameter (e.g., 
aspect ratio), it is shown how a more accurate theory can be de- 
veloped by carrying additional terms for both subsonic and 
supersonic speeds. This theory of not-so-slender wings is 
ipplied to some practical wing problems, including direct prob- 
lems of flow past given wings and problems of wing design 
for minimum drag. The accuracy of the new results is assessed 
by comparison with linearized supersonic-airfoil theory for the 


special case of a flat delta wing. 


(1) HtstorICcAL INTRODUCTION 


— APPROXIMATE THEORY of the flow past slender 
bodies, elongated in the direction of flight, which 
isnow generally known as ‘“‘slender-body theory,’’ was 
discovered by Munk! in 1924 and was applied by him 
to the calculation of the forces on airships in low-speed 
fight. Its extension to wings, which gave rise to a new 
branch of airfoil theory, was made in 1946 by Robert 
T. Jones.2> Munk and Jones considered bodies and 
wings, respectively, flying at small angles of attack and, 
in view of the elongation in the flight direction, made 
the approximation that the flow pattern near the body 
in any transverse section is the same as in two-dimen- 
sonal incompressible flow. Jones’s important dis 
covery was the method of applying this approximation 
to cases involving trailing-vortex sheets. 

Tsien, in an important paper of 1938,* pointed 
out the remarkable fact that the Munk airship theory 
applies without modification to the flow about inclined, 
pointed projectiles at supersonic speeds. Subsequently, 
Jones made it clear in his paper that his slender-wing 
theory applied to both subsonic and supersonic flight, 
atleast for pointed plan forms. To be sure, the mean- 
ing of “‘slender’’ is a bit different in the various speed 


Presented at the Aerodynamics Session, Annual Summer Meet- 
ing, I.A.S., July 16-18, 1952. 

* The authors wish to acknowledge the aid given them by their 
dlleagues at the Graduate School of Aeronautical Engineering, 
especially Prof. N. Rott, and also by A. H. Flax and H. R. Lawr- 
ence, of the Cornell Aeronautical Laboratory, Inc. Part of the 
work reported here has received the financial support of the U.S 
Air Force Research & Development Command. 

t Research Assistant in Aeronautical Engineering. 

t Director, Graduate School of Aeronautical Engineering. 


namely, for any supersonic Mach Number, 
means that the body lies well within the 
This, of course, is 


régimes 
“slender”’ 
Mach cone from the nose or apex. 
the reason for the limitation to pointed bodies and 
wings. It also leads to the conclusion that relatively 
blunt bodies and wings may qualify as ‘‘slender’’ at 
low supersonic speeds, while for very high supersonic 
(i.e., hypersonic) speeds the method must fail entirely 
for practical shapes. In the subsonic régime, it is easily 
verified that the word “‘slender’’ becomes less restrictive 
as the Mach Number increases from 0 toward 1, until 
at sonic flight speed all bodies and wings within a wide 
category become “‘slender’’ ones. 

The well-known Prandtl-Glauert differential equation 
for the perturbation velocity potential in either sub- 


sonic or supersonic flight is 


t= M*) rz + Pyy + ?:2 = 0 (1) 


where x, y, and z are rectangular Cartesian coordinates, 
as in Fig. |, the x-axis being aligned with the parallel 
stream far upstream, and where W/ denotes the stream 
Mach Number. The Munk-Jones approximation con- 
sists of neglecting the first term of Eq. (1) in comparison 
with the second and third. The statement at the end 
of the preceding paragraph means that, for a given 
wing or body, this term continually decreases, relative 
to the terms retained, as / increases from 0 toward 1. 
(Since ¢,, increases with increasing / for any given 
body, this conclusion is not so obvious as might at first 
appear, but it can easily be checked by means of the 
Prandtl-Glauert similarity rule.) 

Robinson and Young*? and later Heaslet, Lomax, 
and Spreiter’ observed that the slender-wing theory of 
Jones yields reasonable and self-consistent results when 
applied to wings without thickness at 1J = 1. This 
conclusion, which might seem to follow immediately 











Coordinates 
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from putting MJ = | in Eq. (1), is, in fact, somewhat 
surprising, since it is known that Eq. (1) itself is not 
generally valid in the transonic régime. The treat- 
ment of wings and bodies with thickness, for example, 
requires that Eq. (1) be discarded in favor of a different 
approximation in which a nonlinear term is retained.° 

Thus the theory originated by Munk and Jones is 
seen to be almost unique within the category of small 
perturbation theories, in that it applies, at least to some 
cases, without modification throughout the whole spec- 
trum of flying speeds from J = 0 to high supersonic 
(but not hypersonic) speeds. Furthermore, it is as 
simple and powerful as classical two-dimensional po- 
tential theory. Naturally, it has been widely exploited 
and extended. By means of conformal mapping, 
Spreiter extended the theory to wing-body combina- 
tions.6 Applications to sweptback plan forms; '! 
rolling;* ° yawing,® and pitching* wings; cruciform"? 
and wing-body-tail'' configurations; and other cases 
have also been made. 

A particularly interesting extension has recently 
been made by Allen,'® based on his experiments on ex- 
tremely long inclined projectiles of circular cross sec- 
tion. Allen’s experiments disclosed the presence of 
trailing vortices analogous to the von Karman vortex 
street of two-dimensional viscous flow. This led Allen 
to the concept of utilizing viscous cross-flow fields in 
the Munk-Jones manner; the same idea has been suc- 
cessfully applied to low aspect ratio wings by Flax and 
Lawrence.'* Finally, in unpublished work carried out 
at Langley Field, Clinton Brown and one of the present 
authors (Adams) have considered analytically the ef- 
fects of flow separation from the leading edges of tri- 
angular slender wings. The resulting flow has been 
approximated by considering a pair of vortices near the 
wing’s upper surface. Within the scope of slender- 
wing theory for triangular wings at supersonic speeds, 
it has been possible to determine the location of these 
vortices and their effect on lift, drag, and pressure dis- 
tributions. 

The agreement between experimental results and the 
predictions of the Munk-Jones theory is generally good, 
provided only that the limitations of slender shapes 
and small-perturbation flow are not violated.'*~' 


(II) PRESENTATION OF SLENDER-BoDy THEORY 


The slender-body theory, in its classical presentation, 


begins with the neglection of the term (1 — 1/*)d¢,, in 
Eq. (1), leaving the approximate equation 
dy + 22 = 0 (2) 


The justification for this approximation is, at least in- 
tuitively, the slender elongated form of the wing or 
body under consideration. 
the geometrical properties of the body or wing vary only 
slowly in the x direction, the derivative ¢,,, which is 
the rate of change of the perturbation velocity com- 


It seems clear that, since 
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ponent in the same direction, must also be small, a 
least in the close neighborhood of the body. Thi 
argument can be made more precise by introducing 
dimensionless coordinates £, 7, and ¢ as follows: j 


x=lt, y=sn, 2 = st 3 


where / is a characteristic body dimension in the x direc 
tion, such as its overall length, and s is a characterist; 
width, such as maximum span of a wing or maximyy 
diameter of a body of revolution (Fig. 1). Writing 


o(x, y, 2) = UIF(E, n, £) 


where LU’ is the stream speed, one has, in place of Ey 
(1), . 
(1 — M?) (s?/I?) Fee + Fan + Fy = 0 (4 


In this form one can consider a family of wings or bodies 
differing from one another only in the value of the 
slenderness ratio s//; it is more convenient to consider 
variations of Mach Number as well and to adopt as the 
slenderness parameter the number (1 — AJ?) (s?//*). It 
is for sufficiently small values of this parameter that the 
first term can be neglected in Eq. (4), thus leading to 
Eq. (2). The fact that this neglection is restricted to 
the immediate neighborhood of the body will be made 
clearer later on in this paper. 

Eq. (2) is exactly the equation for the velocity poten- 
tial in two-dimensional incompressible flow; the con 
clusion to be drawn is, therefore, that the perturbation 
potential ¢(x, y, z) has the form 

$(x, ¥, 2) = oly, 2; x) + g(x) 


where ¢(y, 2; *) is a velocity potential for two-dimen- 
sional flow in a yz-plane, in which x appears only asa 
parameter introduced by the geometry of the cross 
section at x and g(x) is any arbitrary function of x." 
All the cases considered by Munk and Jones—and, 
in fact, by all authors cited up to this point fall into 
the category of lifting cases. These are cases in which 
the effects of incidence and camber alone are being cal- 
culated and the effects of thickness ignored, either be- 
cause the wing under consideration has no thickness 
or because thickness effects do not affect the lift and 
are to be evaluated separately by a different method 
They are, moreover, just the cases for which a knowl- 
edge of g(x) in Eq. (5) is not needed, since this term 
makes only a uniform contribution to the pressure 
field at any x (in the small-perturbation approxima- 
tion), which does not affect the lift distribution. In 
fact, in many of the technically interesting cases in this 
category it is clear that g(x) is equal to zero because 


* Since functions of the forms yh(x), sk(x), and ysi/(x) are also 
admissible in Eq. (5), g(x) perhaps should stand for all such terms 
whose second derivatives with respect to y and with respect to* 
vanish. But actually terms of these forms do not appear, 4 
will be shown later; hence, it is safe to consider g(x) to be a func- 


tion of x alone. 
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g(x, )) 2) must be antisymmetric with respect to the 


plane 2 = 0. 


(III) DETERMINATION OF g(x) 


There are important problems, on the other hand, 
where g(x) is different from zero and where a knowledge 
of this function is desirable for the calculation of pres- 
sure distribution and drag. It has not always been 
clear, in the literature of the slender-body approxi- 
mation, why the Munk-Jones approximation repre- 
sented by Eq. (2) seemed to fail in these problems. 
When the solution of Eq. (2) is written in the form of 
Bg. (5), it becomes evident that the failure is not so 
much a fault of the Munk-Jones approximation as it 
is due to one’s inability to determine g(x). One could 
determine this function by bringing in the additional 
boundary conditions ¢ = ¢, = 0 at large lateral distance 
from the wing or body—specifically, at an infinite dis- 
tance if the flight speed is subsonic or at the nose Mach 
cone if the flight speed is supersonic. But solution (5) 
is invalid at such large distances. 


Ward’s Treatment of Pointed Bodies in Supersonic Flow 


Light was shed on this dilemma by Ward" in his 
important and general discussion of the flow past slen- 
der, smooth, pointed bodies at supersonic speeds. Ward 
shows that the additional boundary condition can. only 
be brought into the picture by returning, momentarily, 
to the complete potential equation, Eq. (1) or (4). 
He attacks the solution of this equation by use of La 
place transforms. Let 


£[F(é, n, ©) |] = Fp, 2, ©) > / F(é,n, Oe” dé (6) 
J0 
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Eq. (4) then becomes 
Pe + F.. = mk? p? F (7) 


where m* denotes ’/* — 1, k denotes s//, and use has 
been made of the boundary conditions 


F(0, n. ¢) = 0 = F;(0, n, £) 


It might be pointed out here that these conditions are 
equivalent to putting F = 0 and F; = 0 on the nose 
Mach cone, in view of the well-known nature of the 
supersonic small-perturbation equation—i.e., the ‘‘rule 
of forbidden signals.”’ 

Having thus introduced the important boundary 
condition at the Mach cone, Ward proceeds toward 
the slender-body approximation by expanding a gen- 
eral solution of Eq. (7) in a form appropriate for the 
1.e., for mkr/—E << 1, 


[This technique is not repro- 


region close to a slender body 
where r = V7? + &. 
duced here in detail because it is similar to that used 
immediately below and also in Section (VI.)| His re- 
sult is, in the notation of Eq. (5), 


U m te asl , 
S’ (x) In ion S”(x’) In (xn — x’) dx’ 
; S70 


g(x) = | 5 
(8) 


where S(x) denotes the area of the body cross section at 
x and it is assumed that S’(x) is continuous. The re- 
striction to pointed-nose bodies requires that .S’(0) = 0. 

Thus, Ward has succeeded in applying the slender- 
body theory, in a unified and consistent manner, to 
smooth pointed bodies of arbitrary cross section, cam- 


ber, and incidence in supersonic flight.* 


* Some of Ward’s conclusions were independently obtained by 


Graham.*! 


It is interesting to compare Ward's result with von Karman’s well-known ones” in the particular case of bodies 


of revolution. 


von Karman’s result for these bodies at zero incidence is 


U ies OED dx’ 
d(x, R) = -— os} (9) 
2x JS0 V (x — x’)? — m?R? 
Ufo s m aes : /\ 
2 —< §’(x) In R + S"(x) In— — S”(x’) In (xn — x’) dx (10) 
Il 2 0 { 


where R denotes the radius of the body at x 


i.e., rR? = S(x). 


solution of the two-dimensional flow problem [Eq. (2) | for this particular case. 


But the first term in Eq. (10) is just ¢(y, 2; x), the 
Moreover, Lighthill'* has shown 


that Eq. (10) is fully as accurate as the original differential equation, Eq. (1), so that this form is theoretically 


equivalent to Eq. (9). 
equivalent to combining von Karman’s with Tsien’s. 


Actually, this combination was usually made by practical engineers 


Thus, Ward’s results include von Karman’s and are, for inclined bodies of revolution, 


i.e., the slender-body approximation was 


used for the effect of angle of attack, but von Karman’s theory based on source-sink distribution was resorted to 


lor the symmetrical flow component. 


It is one of the virtues of Ward’s work that it unifies these seemingly differ- 


ent theories, besides extending the whole subject to a wider class of shapes. 


Extension to Slender Bodies in Subsonic Flow 


The same sort of analysis (for slender bodies of arbitrary cross section, camber, and incidence) can be carried out 


lor subsonic speeds of flight. 


Again, one goes back to the complete Prandtl-Glauert equation, Eq. (4), but in this 
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case one transforms its terms in the Fourier, instead of Laplace, manner. Let 


F[F(é, n, ()] =F (a, n, ¢) =f/ F(é, n, f)e'* dé (1] 

Eq. (4) then becomes 
Fn + Fre = B°R?w?F (12 
where now 6° denotes | — A/* and k denotes s// as before. The following boundary conditions have been applied 


in obtaining Eq. (12): 


F(+o,7, 6) = 0 = F(+, 9, ¢) 


From here on, the solution of this problem is practically the same as Ward’s. A solution of Eq. (12) is 


Flw, n ¢) = > K,,(BR|w/r) [A,,(w) cos n@ + B,(w) sin 16] (13 
0 
where the A,, are modified Bessel functions of the second kind and @ = tan~!(¢/n)—.e., r and @ are cylindrical 


coordinates in the £,,f-space. 

Since the correct boundary conditions far from the body have already been satisfied, one can now proceed to 
expand solution (13) in a form suitable for the region close to the body. The expansion formulas for the K,,’'s for 
small arguments are well-known;' if the leading term only is carried in each A,,, the result is 


n n ’ |” 


Ww “if ® r 


- l =. 2\"| A,(w) cos né B,(w) sin 6 
V (a) -) aw as bl slp Se 2 7 : 
F(w, n, £) Ay(w) | In (; Bk wir) + ¥ T5 2 (n (=) a , 


where y is Euler’s constant (y = 0.57... ). The Fourier transforms having served their purpose, one now returns 
to the variables &, n, and ¢ by taking the inverse transform of each term; for the region near the body, then, 


= “1 

* 1 {3 ; , 

F(é, n, ©) + —ady (€) In ( Bkr) + ay’ (&’) In (& — &’) dé” — ay’ (&’) In (£’ — &) dt’ + 
2 2/0 2S% 


| l _ cos né sin 6 A 
ao(O) In € + — ao(1) In (1 — &) + z a, (&) + b,(&) 15 
») ») m 


al , 
] / } 


where d(£) is the inverse transform of Ap(w). The new functions a,(&) and },(&) are related to A,,(w) B,(w) and 
remain arbitrary at this stage, to be determined by the boundary conditions at the body surface. 
The function ao(é) is easily evaluated. From the boundary condition at the body surface, it follows that 


f, (0¢/Ov) dr = US'(x) (16 


where the contour integral is taken around the body contour C in a cross section, v is the normal to this contour 
and S is the cross-sectional area. In terms of the coordinates &, n, and ¢, this condition becomes 

f{ (OF Ov) dr = (1/1°)S'(€) lj 
where now the normal v and the element of arc length 7 are both in the £&,y,f{-space and the contour integral 1s 
taken around the body contour C’ in this space. Also, the cross-sectional area is now considered to be a function 
of the dimensionless variable £, rather than x, and S’(£) denotesdS dé. But F(é, n, ¢) in Eq. (15) has just the form 
of a plane incompressible-flow velocity potential, and the contour integral in Eq. (17) is proportional to the net 
flow of an incompressible fluid out of the contour C’ in a plane-flow problem. Thus, it can be evaluated at any 
convenient contour enclosing C’ and the £-axis, and its value is —27 a(~).. Hence, 


—ag(E) = (1/2?) S’(E) 


The first term and the summation in Eq. (15) are recognized as the two-dimensional potential function ¢() 
| g 
z; x) of Eq. (5). Therefore, introducing the result obtained in Eq. (18), one has for g(x) the following result: 
U B 1 | , yar sia 
g(x) = — | S’(x) In= — — S’(0) Inx —— S’(/) In (J — x) - S"(x’) In (x — x’) dx’ + 
, 24 Z 2 pe 20 


S"(x’) In (x’ — x) dx (19 
*») x 


Again it has been assumed that the body contour is smooth—i.e., that S’(«) is continuous from bow to stern. 


. . . . . . . - mw 
* Apparently this result has also been obtained by Heaslet and Lomax, since they mention it (without proof) in reference - 


where Ward’s results, mentioned above, are obtained by a different procedure. 
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The first term in Eq. (15) represents, in the two-dimensional flow plane, a source at the axis; the summation con- 
sists of the contributions of doublets and higher degree singularities, all located at the axis. It is well known 
that, for any given contour C within a general class, this arrangement of singularities is not unique but is equivalent 
toa distribution of sources around the contour, for example. In practical problems it is likely that ¢(y, 2; x) will 
actually be found by use of conformal mapping. 

Itis important to notice that the flow conditions at any station on a body in subsonic flow depend upon the whole 
body shape; this was to be expected and is confirmed by Eq. (19). In supersonic flow, of course, downstream 
effects are not felt upstream, so that Ward's g(x) is integrated only from the body nose to the station under con- 
sideration. This difference between the two speed regimes makes an essential difference in the treatment of 
bodies that have blunt sterns; these are outside the scope of the subsonic theory given here but are not outside 
the scope of Ward’s supersonic theory. 

Just as Ward's results constituted, for bodies of revolution, a unification of von Karman’s and Tsien’s, the re- 
sult given here for subsonic flow is equivalent to a combination of subsonic source-sink theory with the Munk 








approximation for the cross-flow effect when applied to bodies of revolution. The source-sink theory gives*" *° 
U in dx’ 
o(x, R) = -— a (') (20) 
tr J0 V (x — x’)? + BR? 
ef a le B lL oo l 
o(x, R) — | S’(x) In R + S’(x) In-— — — S’(0) Inx — - S’() In @ — x) - 
2a 2 2 2 


» > » 


I  aineas , , | pre , P 
S”(x’) In (x — x’) dx’ 4 S”(x’) In (x’ — x) dx (21) 
Ze 0 2) x 


Thus, the slender-body theory of Munk and Jones, when presented in the form of Eq. (5), clarifies the relation 
ship between the familiar subsonic and supersonic versions of source-sink distribution theory for slender bodies of 
revolution. That is, a given body of revolution in subsonic axisymmetric flow has the same ¢(y, 2; x) as in super- 
sonic flow [namely, the first term in Eqs. (10) and (21)] but a different g(x). Asa matter of fact, this can also be 
deduced directly from the source-distribution integrals [as in Eq. (10)| and has already been noticed by some 
investigators.** The present study, together with Ward's, shows that this relationship extends to a much wider 
class of bodies. 

Since g(x) in both of the cases investigated above turns out to be a function of Mach Number M, it is not sur- 
prising that one cannot obtain it from Eq. (2), which (together with its boundary conditions) is independent of M. 


(IV) DETERMINATION OF PRESSURE AND FORCES Satisfactory agreement of slender-body results with 


; . . : either experiment or with the results of more ex: 
lhe pressure is calculated from ¢(x, y, z) according ither expe " with the result 1 more exact 


to Bernoulli's equation. In view of the approxima- 
tions of small-perturbation theory already made in Eq. 


> 


theories depends upon the use of Eq. (23) instead of (22) 
in cases where bodies, as distinguished from wings 
|), there is no purpose in using Bernoulli’s equation alone, — involved. Many writers reached this con- 
in its exact form. For wings (“planar problems’’), o— gl sebittiggee by a - ne . — 
lorexample, the consistent approximation is vpammag te iogirenn ey - nagredinn maarontaged 

theory for right circular cones. The same matter is 
p — po = —pUd, (22) involved in the case of subsonic flow past slender ellip- 


. _— soids of revolution presented in Fig. 2. The soli 
On the other hand, as is well known by now, this is not ids of revolut presented in Fig The solid 


+o ; ae ‘ curves in this figure were calculated by Mz ws" by 
a satisfactory approximation for slender bodies (“‘non- pasiaie gu ere calculated by Matthews” by 


use of exact incompressible-flow theory corrected for 
Mach Number by means of the extended Prandtl- 
Glauert rule (the Géthert rule). The dash-dot lines 
p — po = —pl Ud, + (1/2) ($,? + ¢.*)] (23) are the results of slender-body theory using Eq. (23), 
while the dotted curves are Matthews’ results for the 
equivalent source-sink theory but using Eq. (22). It is 
seen that slender-body theory agrees satisfactorily with 


planar problems’’); a more accurate formula for these 
cases is, in wind coordinates, 


Moreover, Eq. (23) is consistent with the approxima- 
tions of Eq. (1) and of the slender-body theory.** * 

In any case, the velocity components are immediately 
available from Eq. (5)—viz., 6, = ¢y ¢: = ¢:—while 


91S given by 


the more elegant method when the proper pressure 
formula is used; this is in contradiction to Matthews’ 
. ' conclusion. 

6. = Oo/e) + EY) — Matthews’ experimental results, also spotted in 
where O¢/Ox denotes the derivative of ¢(y, 2; x) with Fig. 2, are not in perfect agreement with theoretical 
tespect to the parameter x. ones. Presumably this is due to the rather large thick- 
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CFic. 2. Pressure distribution on 1:6 ellipsoid of revolution at 
zero incidence at various Mach Numbers. Solid curves: exact 
incompressible values corrected by extended Prandtl-Glauert 
rule. Dash-dot curves: slender-body theory using Eq. (23). 
Dotted curves: source-sink theory using Eq. (22). Circled points: 
Matthews’ experiments.”*> Bottom: body contour. 
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Fic. 3. Example of wing-body combination treated by Spreiter.® 
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Fic. 4. Examples of configurations terminating in a baseplane. 


ness ratio of the ellipsoid——viz., 1:6. Matthews finds 
better agreement for a 1:10 ellipsoid. 


Lift 

The distribution of the lift along the length of a body 
was calculated by Munk in his original paper! by con. 
sidering the rate of increase of the transverse compo- 
nent of the fluid momentum (impulse) in a coordinate 
system at rest in the fluid. The same principle applies 
to wings and wing-body combinations, provided only 
that proper account is taken of the trailing vortices, 
The conclusion is always that the total transverse force 
(resultant of lift and side force) is vectorially equal to 
(—U) times the impulse per unit x in the two-dimen- 
sional flow field at the rear end of the body: 


= oU $c, pvdr (25 


where C; is the contour of the wing, body, and trailing. 
wake configuration at the stern and v is the unit vector 
normal to this contour, directed out of the body, into 
the fluid. 

Eq. (25) has an accuracy that may be surprising. It 
properly accounts for the effects of ¢,* and ¢,", as in 
Eq. (23), when they are needed. Thus, by neglecting 
¢,” in a wing-body problem, Spreiter® found a discrep- 
ancy ‘between Eq. (25) and the result of integrating 
pressures to determine lift. His error is a relatively 
obscure one, since it appears in the calculation of wing 
lift where ordinarily Eq. (22) is sufficient. But the con- 
figuration he studied (Fig. 3) involves components 
o, due to wing incidence, body incidence, and body 
divergence on the wing. Since the first two of these 
are antisymmetric with respect to the plane z = 0 
while the third is symmetric, the quantity ¢,* contains 
an antisymmetric cross-product term, of the same order 
contemplated in Eq. (23), which contributes appre- 
ciably to the lift. That this error exactly accounts for 
Spreiter’s lift deficiency mentioned above has been 
verified by Harold Mirels, at Cornell University. The 
same criticism of Spreiter’s paper has independently 
been made by Miles.” 

In those special cases where the wing trailing edge 
and body base lie entirely in one transverse plane, such 
as in Figs. 3 and 4, the total lift, according to Eq. (25), 
depends only on the configuration in this plane. More 
generally, chordwise lift loading at any station x on a 
wing, body, or combination thereof, which lies forward 
of every part of the wing’s trailing edge, depends only 
on the geometry at x. 

In many practical problems, however, it is necessary 
to calculate the loading and the lift on parts that are 
affected by trailing vortices. Frequently, the deter- 
mination of the trailing-wake strength is a nontrivial 
problem. Consider first a wing of arbitrary twist 
a(x, y) lying close to the plane z = 0 (Fig. 5). The 
two-dimensional potential problem for a station such as 
B-B involves knowledge of the trailing-vortex distribu- 
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Slender wing with curved trailing edge 


Fic. 5. 


tion from all the wider sections farther forward. The 
boundary value problem for such a section is: To 
determine ¢(y, 2; x) such that 0¢/0z = — Ua(x, y) for 
—¥ < y < y, while ¢ along the strips y, < |y| <s/2 
has the same values as at the trailing-edge points di- 
rectly upstream—e.g., ¢(P) = ¢(P’) in Fig. 5. Thus, 
the solution at B-B involves knowledge of the solutions 
at upstream stations, and, in general, the values of ¢ 
at the trailing edge C-C, required in Eq. (25), depend 
upon the solutions of all these boundary-value prob- 
lems for the stations between A-A and C-C. Fortu- 
nately, boundary-value problems of this type are not 
particularly difficult; however, it is difficult to reduce 
Eq. (25) for such cases to a more explicit form. 


There is another class of problems where the trailing- 
vortex configuration is important—namely, problems 
involving bodies whose after parts contract downstream 
of the wing trailing edge, as in Fig. 6. Here the alter- 
ation of the trailing-vortex positions due to the body’s 
convergence must be accounted for. The vortices must 
follow the stream lines of the flow around the body; 
in the present approximation they are displaced only 
radially, and this displacement may be calculated from 
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the potential ¢(y, z; x) contributed by the body. In 
the problem illustrated in Fig. 6, for example, the radial 
velocity component is 


rey) - Oy 
or Or 


_dR?* | 
l 


dx r 


so that the radial displacement of a vortex originally at 
ya (in Section A-A), when it arrives at B-B, is given 


by 


> 
ax 


. : “sp dR? ; : . 
VB" — W¥a° = ax = R;? = Ra? (26) 


The lift on the portion of the body.aft of the wing 
trailing edge is given by the difference between the 
values of the integral in Eq. (25) at the trailing edge 
and at the rear end of the body. 


Alan = pU fg gv: dr — ¢ gv, dr) (27) 
B B 4 4 


where v, denotes the vertical component of the out- 
wardly directed unit normal vector ». 

Consider, for example, a case in which the body closes 
to a point at B-B. Eq. (27) becomes, upon introduc- 
tion of Eq. (26), 





Fic. 6. 

















Wing-body combination with noncyclindrical afterbody. 
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ALan = —pU | / gv. dr + (vue 
body at wing at 


A-A A-A 


td 7 
Pratiom)( _ Vy? - =) dy | (28) 


Obviously, this formula is based on the assumption that 
the vortex wake does not roll up and is not otherwise 
distorted or shifted before it reaches B-B. If the body 
is extremely long between A-A and B-B, this assump- 
tion may be unjustified. The motion of trailing vor- 
tices near cylindrical bodies of revolution has been 
studied by Lagerstrom and Graham.” 

It is clear from Eq. (27) that, in the absence of rolling- 
up or other secondary effects, there is no lift on a cylin- 
drical afterbody. 


Drag 

Ward has given a rather comprehensive discussion of 
the drag of slender bodies in supersonic flow in the paper 
mentioned above.'® It may be useful to summarize his 
results here and to try to relate them to the familiar 
engineering concepts of wave drag and induced drag. 
Ward’s result may be rearranged in the form 


D=D,+D.+ Dz (29) 


where 


| 1 si ; 
D = - Zt f / S’(x)S"(x’) In |x — 
tor - 0 0 


7; 


x’| dx dx’ — 2s) S"(x’) In (1 — x’)dx’ + 
Jv 


n 


[.S’(2) ]? In (30) 


= 
_——_ 


») 


— =— —— (3 
. 2 m ( Ov “a oF 


Dp = S(l) (bo — pp) (32) 


As before, C; denotes the contour of the body section 
at the stern, which, in the cases of wings and wing-body 
combinations, must include the trailing-vortex wake. 
Here pz is the base pressure, so that Dz, is the base drag, 
which presumably must be evaluated from experiment 
or viscous-fluid theory. 

The first part of D,; is recognized as von Karman’s 
famous result!” for supersonic pointed bodies of revolu- 
tion at zero incidence. Thus Ward has generalized this 
result, both for bodies not of revolution and for bodies 
where S’(/) is not zero. Moreover, in Ds he has in- 
cluded the contributions due to incidence and lift. It 
is interesting that the term [.S’(/)|? In (m/2) in D, is 
the only term affected by Mach Number within the 
supersonic régime. 

In engineering terminology, D, is the wave drag, and 
Dz might, in most cases, correctly be called the induced 


drag. The exceptions are, again, the cases where S’(/) 
is not zero, for which ¢(y, z; /) includes a source term; 


for other cases, D, is just the kinetic energy produced 
in a unit cross-flow plane at the stern per second. Fo, 
wings, D, includes the vortex drag and accounts cor 
rectly for the leading-edge suction. When the leading. 
edge drag discovered by Jones,” which occurs with 
rounded leading edges, is involved, it is correctly ae. 
counted for in Ds. 

Ward has also derived an especially interesting form 
for drag at incidence: the drag is equal to its value at 
zero lift plus one-half the product of the lift and the 
angle measured from zero lift: 


D = Dizon + (1/2)aL (33 


The simplicity of this general result —and especially the 
absence of cross-product terms —seems remarkable. 

In subsonic flow, the drag situation is even simpler, 
of course, since bodies are required to be closed ut both 
ends and since there is no wave drag. The drag formula 
for this speed régime is simply 


D = DD. (34 


where D, is defined in Eq. (31). Since the body closes 
at its stern, D, arises from the trailing wake and is truly 
induced drag. Again, the formula correctly accounts 
for leading-edge singularities. 


(V) APPLICATION TO UNSTEADY FLow 


The original slender-body idea of Munk was actually 
to replace a steady three-dimensional flow: by a non- 
steady two-dimensional one. Namely, fixing his atten- 
tion upon a transverse lamina of fluid through which the 
body was flying, Munk noticed how the impulse of flow 
in this lamina changed with time because of the pass 
age of the body through it. In this sense, it seems clear 
that unsteady three-dimensional problems do not re- 
quire any essentially different treatment but only a 
simple extension of Munk’s idea. Specifically, one may 
still treat the flow in any transverse plane as two-dimen- 
sional and incompressible, now applying boundary 
conditions that are not only local in terms of x but also 
in terms of ¢. Thus, the calculation of the potential 
distribution ¢(x, y, 2, ¢) is the same as in steady flow, 
and only the formulas for the pressure, and therefore 
the forces, are different. 

The approximation involved becomes clearer if the 
differential equation for ¢@ is written down. Two new 
terms must be added to Eq. (1)—viz., 


[1 — (U?/a?) oir + by, + 622 = (1/a?) (Gu + 
2U¢d.1) (35 


where a denotes the speed of sound and UL’ is the undis 
turbed stream speed, both constants. The present ap- 
proximation is to discard all of the terms in Eq. (35 
except the last two on the left-hand side, whence, using 
a notation analogous to that of Eq. (5), 


Puy + o> -_ 0 t 
g(x, y, 2, t) = oly, 3; x, t) + g(x, dS 


(36) 
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SLENDER-BODY 


To evaluate this approximation, consider, for ex- 
ample, a problem involving sinusoidal! motion for which 
it may be assumed that ¢(x, y, 2, f) is proportional to 
e"; Eq. (35) becomes 


) 


eg + 25” Moe + (M*— l)bee = by + x2 (37 
a’ a 


where ./ is the stream Mach Number. In terms of the 


y/s, and ¢ = 2/s introduced 


variables £ = x//], n = 
above [Eq. (4)] and writing «x for the usual ‘‘reduced 
frequency’’ n// LU’, this is 

_ > 

J = F,, + F.. 


ft 
ss 


ST (-eMeF + QicM2F, + (M2? — 1) Fe 
(3S) 


where 


g(x, vy, 2, t) = UIF(§, n, oe" 


Now, assuming s << / as in slender-wing theory, one 
can see that the approximation represented by Eq. (36) 
is correct, provided only that neither J/ nor « is large 
compared to 1. This is the correct approximation for 
many practical cases, such as flutter of low aspect ratio 
wings in any speed range below the hypersonic. Miles*’ 
has extended Ward’s slender-body theory on this basis. 

Recently, however, Lomax, Heaslet, and Fuller*! 
have presented a different sort of slender-wing theory 
viz., they approximate to Eq. 


a \Drr 


for unsteady flow 
by neglecting only the term ¢,,, but not (—l 
or the right-hand terms—1.e., they write 


(35) 


_ 


bu + 2Ug, + U*brr = a*(hyy + $22) (39) 


This theory was subsequently criticized by Muiules* 
as being inconsistent. Miles stated that the approxi- 
mation of Eq. (36) should have ben used. While 
Lomax, et al., do not explain why they have retained 
the left-hand terms in Eq. (39), a plausible reason may 
be advanced: They are interested in “‘indicial’’ prob- 
lems—i.e., transients--rather than oscillations; it is 
well known that these require, in effect, adequate treat- 
ment of all values of x in Eq. (38), including « >> 1. 
Consequently, the approximation of Eq. (36) is unsatis- 
factory, and the first two left-hand terms in Eq. (38) 
must be retained. To be sure, the term (s?//*) (7? — 
|)F might well be dropped for any moderate 1/. But 
no harm is done and the order of approximation is not 
changed if a part of this small term—namely, (s* 
P) MF 
simplification results, for then, upon transformation 
to a coordinate system fixed with the stream, the ve- 
locity potential is found to satisfy the wave equation 


is retained; a considerable mathematical 


for two-dimensional motion—i.e., 


bi = A7(dyy + bez) (40) 


The result is an approximation quite analogous to 
Munk’s and Jones’s but employing the acoustic equa- 
tion, Eq. (40), in place of the Laplace equation in 
transverse planes in a coordinate system at rest in the 
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fluid. It appears from the analysis above that this will 
be more accurate at high frequencies and in transient 
problems. 

Thus, it may be concluded that Lomax, et al., are 
justified in carrying certain time derivatives of ¢ and 
that they carried another ‘“‘inconsistent’’ term to 
achieve a mathematical simplification. This term is 
small and has no significance to the order of approxi- 
mation involved so long as the wing is slender.* It is 
interesting to note that Eq. (40) is the same as for 
small-perturbation hypersonic flow** past bodies and 
wings, the latter not necessarily in the slender-wing 
category. 

(In a private communication, Dr. Miles has also 
pointed out an error in the application of this new theory 
to the transient motion of a delta wing in reference 31 
viz., the increase of span with time in the spaced-fixed 
coordinate system was neglected in setting up the 
boundary-value problem for Eq. (40). This error does 
not, of course, affect our conclusions regarding the 
validity of the theory itself.) 


(VI) EXTENSION TO ‘‘Not-So-SLENDER"’ WINGS 


It has already been noted in Section (III) how the 
slender-wing theory may be obtained by expanding 
the velocity potential of linearized theory in powers 
of a slenderness parameter |such as Bk in Eqs. (12) 
(15)] and retaining only the largest terms. An essential 
extension of the power of the slender-body theory seems 
to result if additional terms of this expansion are re 
tained. The methods of Ward and the analogous meth- 
ods for subsonic flow, both presented in Section (IIT), 





are appropriate for this purpose. 

One further generalization of these methods is re- 
quired, however, in order to make this idea useful. This 
is to distribute line singularities at various spanwise 
distances rather than only at the axis. 


Subsonic Flow 
This extension, which might be called the theory of 
wings, will now be outlined briefly 
In what follows, only pianar prob 
i.e., wings having small incidence, 


‘“‘not-so-slender”’ 
for subsonic flow. 
lems are considered 
camber, and thickness and lying nearly in the plane 
z= 0. 

Once again, for subsonic flow, Fourier transforms are 
used, and the Prandtl-Glauert Eq. (4) is reduced to 
Eq. (12). A solution analogous to the first term of 
Eq. (13) but permitting the placement of singularities 
at arbitrary spanwise locations n = 7’, say, is 


F(w, 0, ©) = fw, 7’) Ko(Bk|w\r’) (41) 


where again Ky is the Bessel function and now r’ de 

9 9 1'/e in ——. 
notes [f? + (n — n’)?] In the &,,¢-space this is 
a line source of variable strength lying along the line 


* It would not be correct, for example, to make the same ap- 


proximation for a nonslender wing near M = 1. 
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n = n’,& = 0. The strength of the line source is pro- 
portional to the inverse Fourier transform of flo, n’); 
say /(&, 7’), and a distribution of such line sources in the 
£, n plane can be used to satisfy the boundary conditions 
of a symmetrical wing of small thickness at zero inci- 
dence. The velocity potential ¢(x, y, z) is obtained by 
inversion of the Fourier transforms. In fact, this in- 
version leads immediately to the familiar surface in- 
tegral over a source distribution of linearized subsonic 
wing theory. 

To obtain the desired expansion in the width pa- 
rameter 6s//, one must restrict the solution to the neigh- 
borhood of the wing and expand Ko in Eq. (41) for small 
values of its argument—viz., 

e < e 
Ko(e) = -(1 + <)in —-ytiQl “Poe 


? 
O(e* In e) 


(42) 


The perturbation potential $(x, y, 2) is therefore ob- 
tained, to terms of the desired order in the width param- 


sU , | r| , . S’(é) 
$(é, 1,0) = — fon) 0 |= | dy +4 ad 


is sUT, Oo? 
S”(é’) In (&’ — &) ae’ | 2p? | 
| g + 6 te | Oe 


l ha ve , , 4 , ”» / , , 
| in (& — 8) de f(y ~ Pale, a! — 
Z Of /0 
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eter, by substituting Eq. (42) into (41), integrating 
. . - 5 

over 7’, and inverting the transforms term-by-term 

The relation between the arbitrary distribution function 


/(&, 9) and the surface boundary value is determined by 


investigating the behavior of the elementary solution 
(41) at the plane z = 0; the result is 


f(é n) = (1 1 )b(&, n, + 0) (43 


i.e., f is proportional to the vertical velocity component 
at the wing surface. Thus, an expression for $(x, y, 2), 
correct to the second order, appropriate for thin sym. 
metrical wings without incidence, is obtained. The 
surface boundary condition for this case is 


¢.(x, y, + 0) = +Uoa(x, y) 
or 

$;-(&, 7, + O) = +sUa(é, n) 
where o denotes the local inclination of the wing surface 
in the x direction. The desired result, evaluated at 


the plane z = 0 in particular, is 


we = : lf, S”(é") In (£ — &’)dt’ — 
2 4rl 0 


=m +) a + 


4 "a . 
I oP nn Sie ( — ‘\2 t/ , 
2 de" J In (é E)dé / n n )*a(&', n in | (44 


where all integrals with respect to n’ are carried across the plan form. 
The first four terms in this formula constitute the slender-wing approximation for the case in question and have 


their counterparts in the results of Section (IIT). 


The first term is ¢(y, 2; x); the next three constitute g(x). It 


has been assumed in Eq. (44) that S’(0) = 0 = S’(/), as is appropriate for most practical cases, and that S’(x) is 


continuous. 


The new theory provides the additional terms in Eq. (44). 


The analogous treatment of lifting wings without thickness is easily accomplished by using a distribution of 


doublets in place of sources. 


tial, Eq. (41), with respect to z—1.e., 


F(w, n, £) = 


h (w, n’) being arbitrary. 


The corresponding potential is found by expanding A, for small argument, integrating over n’, and inverting the 
In this case, the doublet strength h(£, 7) is proportional to the surface potential 


transforms. 


The surface boundary condition, however, ordinarily involves ¢, rather than @ 
o:(x, y, 0) = 


consequently, the integral expression*for ¢(x, y, 2) must be differentiated, and an integral equation for h(é, 7) re- 


Be / 1 Bk 
Qs (- .*™ ) 


sults: 


l , ay’ 
fiw é, n ) = — 
§ 7 ~ 9 


Ualé, n) = 


l d* o / , / , , l d? ; ; 7 
sah In (& — &’) dé free yn’ )dn +3 ef In (&’ — é)dé fie, n’) in’ | (47) 


Here, the integrals with respect to 7’ must be carried across the wing plan form and also across the trailing-vortex 


— Bk| | + h(w, n’)Ki(Bk|w|r’) 
r 


F (1/m)(&, n, +0) 


The doublet solution is found in the usual way by differentiating the source poten- 


ViZ., 
(46) 
ie: 


— Ua(x, y) 
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wake at stations farther aft than the position of maximum span. 
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In this case the Jones theory is represented by 


the first integral alone, all the subsequent terms being provided by the new technique. 
The procedure for solution of this integral equation for h(g, 7), which may seem rather formidable at this point, 
will be discussed immediately below after stating the analogous formulas for supersonic flow. 


Supersonic Flow 


The ‘‘not-so-slender’’ wing theory for the supersonic case is carried out exactly as described above, except that 


Laplace, instead of Fourier, transforms are used. 


are. 
(a) For symmetrical wings without incidence, 
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This treatment has been outlined in reference 34. 


The results 
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where it has been assumed that .S’(0) = 0 and that .S’(x) is continuous. . 


(b) For lifting wings without thickness, 


[ 1 f ‘ dn’ m’*k? 
Ua(t, n) = hy (& 9°) = + = 
S yn —- 9 25 


[(-3 +m 


Solution of the Integral Equations 


As has already been mentioned, the customary ‘‘di- 
rect’’ problems of wing theory require that Eqs. (47) 
and (49) be solved as integral equations for /(£, 7). 
The difficulty of this mathematical problem disappears 
immediately when it is remembered that solutions to 
order B°k? (or m*k*), only, are being sought. To clarify 
this, imagine the doublet distribution / to be written 
explicitly as 


h=h® + h@ (50) 


where h‘” is the doublet strength in Jones's slender- 
wing approximation, which is independent of 6k, and 
h® includes the higher order terms, presumably of 
order 6°k* but the exact order need not be assumed 
here. Upon substitution of Eq. (50) into (47) and re- 


calling that 
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(by definition of h"), one has, in place of Eq. (47), 
3 dn’ ore = ? —. 
hy (, 9’) — = Bk? X (functions of h)| 


. 7 —7F ( 
= Bk? X (functions of h“) 


(52) 


where terms of order higher than 8°k? have been neg- 
lected, as originally proposed. 

Thus, in place of the formidable integral equations, 
Eqs. (47) and (49) for h, one has to solve only equations 
of the form of Eq. (52) forh®. Fortunately, this form 
is one for which explicit solutions can immediately be 


written down. It is, in fact the same equation as oc- 
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curs in Jones’s theory and elsewhere in thin-airfoil 
theory.*® 

To date, this second-order theory has been applied 
to only a few practical cases. The results of these ap- 
plications will now be presented briefly, and some fea- 
tures of the method will be pointed out. 


Example 1: Flat Delta Wing 


The flat delta wing at incidence, in the supersonic 
régime, presents an ideal problem on which to test the 
new theory, since for it, by virtue of its conical sym- 
metry, results of the complete Prandtl-Glauert theory 
are available for comparison.* The calculation pro- 
ceeds without difficulty, treating Eq. (49) in the man- 
ner sketched above. The significant results are sum- 
marized in Fig. 7, where the lift-curve slope divided by 
k, C./ak is plotted against the breadth parameter 
mk/2. When mk/2 = 1, the leading edge of the delta 
coincides with the tip Mach cone, and the wing surely 
cannot be considered slender. Jones’s theory gives the 
value a for C,/ak for all k. 

Also spotted in Fig. 7 are experimental values taken 
from reference 14 and (one point) from reference 2. 
It may be concluded, on the basis of this test that in- 
clusion of second-order terms in the breadth parameter 
provides an important improvement of Jones’s theory 
over a considerable range of this parameter. 

The present example is not well suited for a test of the 
improvement of pressure distribution afforded by the 
second-order theory, since all three theories give the 
same spanwise and chordwise distribution functions. 
In other words, a comparison of pressure distributions 
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L'(x) de- 
Dash line: 


\/2 (upper) and M = 0 (lower). 
Solid line: present theory. 


sketched, at M = 
notes the lift loading. 
Jones’s approximation. 


at any point of the plan form leads only to the results 
plotted in Fig. 7; the pressures bear the same ratios to 
one another as C,/ak for any given mk/2. 


Example 2 (Fig. 8) 


It is interesting to apply the new method to the case 
of a wing such as the flat-plate wing shown in Fig. 8, 
* . 9 . . . 
for which Jones's theory predicts no lift loading at sta- 


96 


tions aft of the maximum span for either subsonic o; 


supersonic speed. 
the plan form of the wing selected for this example has 


For reasons to be explained later 


continuous slope, curvature, and third derivative at the 
station of maximum width. The calculated chordwige 
lift loading is plotted in Fig. 8 for both subsonic anq 


supersonic speeds. 


Example 3 (Fig. 9) 


A certain difficulty, apparently an essential one jn 
some cases, arises in the application of the not-so- 
slender wing theory to wings whose chordwise lift 
loading according to Jones’s theory is discontinuous 
or has discontinuous slope or curvature. This trouble 
arises in the last integrals in Eq. (47) or (49). The re. 
sults in such cases exhibit infinite values of h near 
such stations of discontinuity.* Naturally, for super- 
sonic wings such discontinuities do not affect any parts 
of the plan form further forward, so that there is no 
difficulty in treating a wing (e.g., Example 1, above 
whose Jones loading is discontinuous at the trailing 
edge. 

To handle the after portions of the flat-plate wing of 
Fig. 9, however, requires a special technique, even in 
supersonic flow. One considers first a fictitious wing 
having the same plan form and the same incidence for 
ward of maximum span but having its incidence a(x 
adjusted along the x direction aft of station AB so as 
to render the Jones loading suitably continuous. For 
the flat wing of Fig. 9, for example, a suitable a(x) is 
quadratic in x aft of AB. Having solved this problem 
by the present second-order methods, one must proceed 
to correct the loading aft of AB for the presence of the 
fictitious incidence. 

* Actually, this difficulty was introduced in the term-by-term 


inversion of the Fourier (or Laplace) transforms, above. In 
this step, smooth behavior of the functions was implied 
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Fic. 9. Chordwise distribution of lift on the flat-plate wing 
sketched, at M = 7/2. (Designations of curves same as in Fig 
8). Dotted line: estimated fairing to asymptotic curve. 
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Let a, a constant, be the true incidence, and let 
a(x) denote a — a(x). If ¢’ is the potential 
calculated by the second-order theory for the fictitious 
wing, the correct potential for the original wing is 


~ 


o(X, y, 3) = eo”) +9" (5: 
where ‘? is the potential for the flow past the rec- 
tangular wing ABCD with incidence a(x). 

Now the calculation of ¢"” is outside the scope of 
the second-order theory, but it lies within the category 
of cases treated easily by the complete Prandtl-Glauert 
theory. In fact, the incidence a‘! (x) being a function 
of x only, a technique due to Rott® is especially suited 
jor this calculation in the region ABEF, and the chord- 
wise lift loading is obtained very simply. Fig. 9 shows 
the result of applying this method to a specific problem. 

It may be remarked that the second-order results 
obtained for such ‘‘discontinuous’’ wings by direct use 
of Eq. (47) or (49), without benefit of the correction 
technique just outlived, are not without some signifi 
cance in spite of their infinite behaviors at the discon 
tinuities. Namely, they are asymptotically correct 
far downstream of the discontinuity (and upstream as 
well, for subsonic cases). In this statement, moreover, 
the word ‘‘far’’ must be interpreted with relation to the 
width, and for supersonic flow, as Stewartson has 
shown,** ‘‘far downstream’’ probably means the dis 
tance of only three or four Mach-wave reflections (Fig. 
9). Thus the technique used here is analogous to Light 
hill’s technique* for slender bodies of revolution which 
are not smooth—1.e., whose meridian profiles exhibit 


discontinuous curvatures. In this spirit the corrected 


D= 


we 


second-order curve in Fig. 9 has not been calculated 
aft of EF but has been faired roughly (dotted curve) 
toward the asymptotic second-order curve. 

It is important to recognize that, although the second- 
order wing theory could not be used directly for the 
“discontinuous” wing, it nevertheless plays an essential 
part in reducing the problem from one that was in- 
tractable in Prandtl-Glauert theory to one that could 
be handled easily. Unfortunately, a comparable tech- 
nique has not yet been worked out for subsonic cases, 
so that to date the theory of not-so-slender wings in sub 
sonic flow is severely restricted in application. 


Supersonic Wings of Minimum Drag 


The present extended slender-wing theory permits 
the determination of wing configurations for minimum 
drag. In subsonic flow, as is well known, all wings with 
elliptic spanwise lift loading are minimum-drag wings 
for given lift and span. In supersonic flow, on the other 
hand, the phenomenon of wave drag appears, and the 
minimum-drag problem is more complex. 

For the classical problem of the distribution of lift of 
a wing without thickness, to produce minimum drag, 
slender-wing theory gives the same result for super- 
sonic as for subsonic flow, because wave drag is absent 
in the approximation. The second-order approxima 
tion introduced here brings in wave-drag terms, so that 
the minimum-drag problem yields information about 
the chordwise, as well as spanwise, loadings. 

The drag for a supersonic wing without thickness 
can be obtained to second order by use of Eq. (49), 
taking account of the leading-edge suction. The re- 
sult is 
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where /(£) denotes f°h;(é, »)dn, so that the chordwise 
lift loading is —2mrpUk/(é). Here 


tions with respect to » are carried across the wing and 


again, the integra 


wake. 

It is observed that there is an almost exact resem- 
blance between the last three terms of this equation and 
Eq. (30) for the drag D, of slender bodies or slender 
wings with thickness. In this comparison the chord- 
wise loading in one case and the slope of the cross 
section area in the other case play the same role. 

A few results that arise from a study of Eq. (54) can 
be mentioned briefly : 

(a) For the case of a wing with the parameter mk < 
0.44, Eq. (54) indicates that minimum drag is achieved 
for given lift, span, and chord when the spanwise load- 
ing is elliptic and the first-order (Jones) chordwise 


loading is also elliptic. This makes both induced and 


wave-drag minima, to second order. The drag in this 


case is given by 
wA) + (m?*k?/2) (C,?/2A) (55) 


where A denotes the aspect ratio, s*/(area). Jones* 
seems to have arrived at this same result, although his 
methods were not disclosed. 

For wings of higher aspect ratio, or mk > 0.44, the 
theory indicates that the chordwise loading should not 
be elliptic and that the wave drag can be reduced by 
maintaining a finite loading at the trailing edge. 


example, a flat delta wing produces less drag than a 


For 


wing with elliptic chordwise loading for mk > 1.1. 

(b) For the special case of conical wings the mini 
mum-drag wing is a flat plate. This is the same result 
as is well-known in slender-wing theory; 


of second-order terms does not alter it. 


the inclusion 
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Fic. 10. A wing satisfying the conditions for minimum drag 
for given lift, span, and chord when mk < 0.44; viz.,a flat plate 
having the plan form show. 


(c) A simple example of a wing with elliptical distri- 
butions as discussed in (a) is a flat-plate wing having 
the same plan form as the meridian contour of von 
Karman’s famous ogive," as sketched in Fig. 10. 
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Penetration and Deflection of Jets Oblique to 
a General Stream 


FREDRIC F. EHRICH* 


Westinghouse Electric Corporation 


SUMMARY 


The dynamics of the entry of jets of fluid into a general stream 
itan oblique angle from slots and orifices is studied by means of a 
simplified two-dimensional potential flow analysis. The wake is 
ipproximated by an extensive area of dead fluid separated from 
the jet by a vortex sheet at uniform pressure. The jet is sepa 
rated from the main stream by a stream line, and the jet velocity 
is restricted to one particular value so that there is no discon 
tinuity across this stream line. Analysis using the conventional 
Helmholtz-Kirchhoff method leads to computations of the various 
relationships between the geometrical and velocity parameters 
Calculations are also made of the shapes of jets in some common 
geometries 

An equivalent analysis is made of the situation in which the 
trailing edge of the jet does not separate from the downstream 
wall, and the results are compared to the analysis with separation. 

The solution for the flow from a two-dimensional orifice is com 
pared to experimental data for the roughly comparable case of 
flow from a circular orifice, and a qualitative agreement is shown. 

The results have a possible utility in serving as a skeletal struc 
ture of a more complete analysis of the mixing process in combus 


tion chambers and cooling films 


NOMENCLATURE 


A’, B’,C’, D’, E’, F’ = points in the jet 
B.C = points in the main stream 
a = width of the jet source 
= offset of the jet source 
( = asymptotic width of the jet in the 
stream 
F = a function of ¢ 
h = a distance in the ¢ plane 


= absolute value of the velocity at a 


point 

Y= In U/q + 0 = logarithm of the conjugate hodograph 
plane 

= distance along the vortex sheet 

l = velocity of the main stream 

l = velocity of the jet in the slot 

‘ = coordinate in the direction of the main 
stream 

y = coordinate in the direction perpendicu 


lar to the main stream 
mx + ty = the physical plane 
angular direction of the forward side 
of the jet source 
8 = angular direction of the rear side of 
the jet source 
= a small distance in the ¢ plane 
= the mapping half plane 
6 = angular direction of the velocity at a 


point 


Received June 13, 1952 
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potential function 


o = 
¥ = stream function 
wo=o+1y = complex potential plane 


INTRODUCTION 


QO" SOME BASIC AND GENERAL INTEREST is an 
analytic description of the penetration and de- 
flection of a jet of fluid entering a large body of fluid 
from an orifice set oblique to the direction of the main 
stream. 

The complete problem is one in compressible, viscous, 
unsteady fluid flow in a complicated three-dimensional 
geometry. In order to make the problem amenable to 
analysis, it is necessary to restrict attention to some 
fluid dynamic aspects in a simplified geometry. 


ANALYSIS 


Making the usual assumptions of steady, incompres 
sible two-dimensional potential flow reduces the prob 
lem to the schematic flow picture shown in Fig. la. It 
is assumed that the densities and total pressures of the 
jet stream and the main stream are equal. This last 
stipulation restricts our attention to the peculiar case 
where the boundary between the leading edge of the 
jet and the main stream is not a line of discontinuity but 
a stream line common to both flows with a continuous 
variation of velocity existing across that boundary. 
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A stream of fluid from an orifice formed by two ob- 
lique walls meets the main stream at a stagnation point 
BB’ and forms a common stream line BC. The jet 
separates from the wall at point /’, and a vortex sheet 
or free stream line E’D’ marks the boundary between 
the jet and an extensive wake assumed to be composed 
of still fluid. 
in the wake, the pressure along the vortex sheet must be 


Since it is assumed that there is no motion 


a constant, and, hence, the velocity must also be con- 
stant. A second possibility is the situation in which the 
downstream boundary of the jet does not separate from 
the wall as depicted in Fig. 1b. This, of course, in- 
volves an infinite velocity of the fluid in turning the 
sharp corner at E’ but is probably approximately 
realizable in certain real flows. 

The flows under consideration are bounded only by 
straight walls and constant-pressure vortex sheets and, 
hence, fit into the class of problems which can be 
treated by the Helmholtz-Kirchhoff method.! The 
analytic procedure consists of a series of conformal 
transformations which links the complex physical plane 
z to the complex potential plane w. The degree of 
difficulty of the various transformations depends in 
general on the number of discontinuities in the flow 
sources, sinks, doublets, stagnation points, points of 
infinite velocity, and points of separation. In order to 
reduce the number of points of discontinuity under con- 
sideration, two specializations of the generalized geome- 


try will be considered. Firstly, a slot-type jet source 
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with a equal to @ will be studied to eliminate the zero 
velocity at the source of the jet A’F’. Also, an orifice- 
type jet source with a equal to zero can be considered, 
since the zero velocity at the joining point BB’ is thus 
eliminated. These geometries are shown in the sum- 
mary of transformation planes in Fig. 2. 

Although the shape of the vortex-sheet boundary is 
not known in the physical plane, it can be drawn in the 
hodograph plane by virtue of the stipulation of constant 
velocity along it. With all the vortex sheets becoming 
ares of circles and all straight boundaries transforming 
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Fic. 3. The mapping functions and analytic relationships 


into radial lines, the boundaries of the flow are readily 


drawn in the hodograph plane as in Fig. 2. 


The conju- 


gate of this plane is functionally related to the physical 


and potential planes, since 


10 


—dw/dz = ge 


Taking the logarithm of the circle sector transforms the 


bounded area to a degenerate rectilinear polygon, 
since 

Q = —In [(q/U)e™™ In U/qg +70 (2) 
The boundaries of constant velocity become vertical 


lines, 
horizontal lines as in Fig. 2. 

The Schwarz-Christoffel transforn 
polygon to be transformed into the 


new plane, ¢, giving the functional 


and boundaries of constant inclination become 


1ation® permits this 
upper half of some 
relationship 


stream lines are horizontal lines by virtue of their con 


stant stream function (as in Fig. 2 
also be transformed directly into the ¢ 
utilizing the Schwarz-Christoffel transformation. 


This plane may 
¢ upper-half plane 
Be- 


cause of the basic similarity of the configurations under 


consideration, they all give the same transforming func- 


tion: 


dw/dé = 


The set of transformations given by Eqs. (1), 


+ h)?Uc/r(e + h(E — 1)? (4) 


(2), (3), 


and (4) completely determines the problem functionally 
and permits the calculation of various parameters. 


CALCULATIONS 


The shape of the vortex sheet may be found by virtue 
of the fact that the velocity is constant so that 





QO = QO(¢) (3) 
which is tabulated in Fig. 3. This function must be j 
° . - . . . as = 
left in terms of a parameter / that indicates the relative 
location of the vertices of the polygon in the ¢ plane and ng 


is a function of the relative orientation of the edges of 
the jet source in the physical plane. 

The boundaries of the flow system may be drawn in 
all 


the complex potential plane w immediately, since 


dz 


dy 


dx 


dz/dwidw/déidt\ = (1 


= ds cos 6 


ds sin 6 = 


U’)\dw dt\dé 


N47 
= [real part (e°) |ds 


° . ) 
imaginary part(e°) |ds 


(5) 


(6a) 


(6b 
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Eqs. (7a) and (7b) give a pair of parametric equations 
in ¢ which determines the shape of the vortex sheet in 
terms of the constant / and are presented in Fig. 3 for 
the geometries in question. 

A similar procedure, with some modification, is effec- 
tive in evaluating the relative orientation of the edges 
of the orifice in terms of the parameter h. An integra- 
tion in the physical plane from edge B’ to EF’ is equiva- 
lent to an integration from B’ to KE’ inthe ¢ plane. But 
an integration in the ¢ plane along the real axis must 
avoid the pole at the point A’F’. This is most readily 
accomplished by the path shown in Fig. 4 where the 
integration is carried out along the real axis to a dis- 
tance (—h — e¢). Then the integration is calculated 
along a semicircle of radius e from (—h — e€) to (—A + 
e), from which point the integration can be carried out 


along the real axis. Thus, 


om —1 
/ ds/¢ = f (dz/dw)(dw/dé)deé = 
J B' ie 


e2(1 + h)%d¢ 
J. E+ ME — 1)? 


(Sa) 


CE 


4 —-h—e 
/ dz/c = lim | f F(¢)dé + 
Jy e—>0 =. . 


0 a1 a 
F(ee” — h)d(ce”)+ j Fide | (Sb) 
™ J —h+e 
The complex quantity given by this equation may be 
split into two geometry ratios, a/c and b/c ,which indi- 
cate the width of the orifice and the offset of the rear 
edge from the horizontal plane of the forward edge of the 


orifice. For the case of the slot-type jet source, part of 
this calculation is unnecessary since continuity requires 
that 


Uc = Usa (9) 


so that a/c is more easily calculated from Eqs. (2) and 
(3) such that 
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In U/U, + ta = QO(-A) (10) 


Although the problems have been set up for general- 
ized values of a and 8, the pertinent integrals cannot be 
They can be solved 
But for 
practical purposes, the integrations can be carried out 
only for a = 2, a = w/2, and 6 = f, utilizing partial 


integrated for the general case. 
only if a and @ are rational fractions of 7. 


fraction reductions. (The cases a = 0 and 6 = 0 are 
trivial.) Fig. 3 summarizes the pertinent analytic rela- 
tionships, and Figs. 5, 6, 7, and 8 show plots of the 
geometry and velocity parameters for the various con- 
figurations under consideration. It must be realized 
that only one particular velocity ratio corresponds to 
any set of specific values of the geometry ratios, and 
this velocity ratio is governed by the criterion of 
smooth entry—.e., it is the only velocity such that the 
total pressure of the jet is equal to the total pressure of 
the main stream. 

Figs. 9, 10, and 11 show typical plots of the shape of 
the jet from various types of orifices. The leading edge 
of the jet or the stream line common to both the jet 
flow and the main flow is not readily calculable, but its 
shape is easily and accurately sketched from the knowl 
edge of its asymptotic behavior. The stream line must 
bisect the corner at BB’ because of the symmetric na- 
ture of a stagnation point. Also, the jet quickly 
reaches an asymptotic width of value c, which quantity 
has entered the calculations of the other geometric 


parameters. 
CONCLUSIONS 


Although the analysis is severely limited in its applic- 
ability because of the extreme nature of the simplifying 
assumptions, certain general conclusions may be ascer 
tained. 

The concept of smooth entry of a jet with no discon 
tinuity in velocity across the stream line common to the 
jet and the main stream is of extreme significance in 
certain applications. A film introduced at such a 
velocity would have little tendency to mix with the 
main stream and, hence, would be extremely effective as 
an insulating layer. Similarly, if mixing is desired, as 
in the case of the introduction of secondary air into a 
combustor, such velocities must be avoided. 

Fig. 11 includes a curve from an empirical formula 
deduced by Callaghan and Ruggeri’ from extensive 
measurements on heated air jets from circular orifices 
perpendicular to the main stream. Their formula, re 
duced to the nomenclature of this paper and specialized 
for the assumptions under consideration, gives 


y/a = 1.91[(x + a/2)/a)}?-* (11 


Of course, this is only a qualitative comparison, since 
the circular jet is by no means directly equivalent to the 
two-dimensional jet under consideration, especially 
since it in no way approximates the extensive wake that 


has been assumed. 





OFFSET RATIO 


% 


DEPTH OF PENETRATION 


i) 









































103 
















































































































































































































































































JETS OBLIQUE TO A GENERAL STREAM 
- _— 
ee, i ene 
u “ Pi 
ral / 
S | —-— 
0 b 
06 WITH ho 
uf SEPARATION a POSITIVE @ POSITIVE 
b NEGATIVE g b POSITIVE 
U a 
04 a a 4 
ae 3 
3 = 
, = g 
NO Ps FE 
z 
- ut | SEPARATION “4 S 2 Pa 
° — 
- 2 
« oO 
b/o OFFSET RATIO 
bo —Z 
= 02 04 06 08 10 12 -6 -4 -2 2 4 5 
Uy, = Sf ET VELOCITY 
x RATIO NO SOLUTION 
-02 + "4 -2 
i) 
RAT 
Ph 
ate ® | 
-04 / -4 
re, 
-06 & _ 
$ oa 
° — 
-08 
/ © NEGATIVE © NEGATIVE 
b POSITIVE b NEGATIVE 
Fic. 5. Smooth entry velocity for a lateral jet. Fic, 7. Contraction ratio from an orifice (with separation). 
——— sx 
/ b ¢ 
08 a — 
- to 
© POSITIVE © POSITIVE 
b NEGATIVE b = POSITIVE 
06 Pa ° 
- 6 
_ 
= Wa 
z 
Pe © 
04 7 S 4@ Z 
s 
a 
<= 
z 
°o 
oO 
02 
2 
é x 
* b 
Yo OFFSET RATI 
$ (*) SE ATIO 
g o 2 4 6 
. a 08 1.0 12 - - = 
JET VELOCITY RATIO 
= 
a 
& -02 -2 
ic} 
> 
-04 ° 
Le} a 4 
aes 
WITH SEPARATION 
-0.6 L -6 
u eens ee (== 
-08 ae =<. 
NO SEPARATION 4 
Oo NEGATIVE G NEGATIVE 
bs POSITIVE b = - NEGATIVE 
Fic. 6. Depth of penetration for a counterflow jet Fic. 8. Contraction ratio from an orifice (without separation 





104 JOURNAL OF 


THRE AERONAUTICAL SCIENCES 


FEBRUARY, 1953 























Fic. 9. Shape of a lateral jet (6/a = 0; c/a = 0.691 

















\ 
\ 
MAIN 
STREAM a JET 
oe ~~ +> —_— 
U hr. ur 
a i a ee cll 
Fic. 16. Shape of a counterflow jet (b/a = 0.358; c/a = 0.800) 











i 
R onriee ~ 
ae 
\ ee 
a _ 
Wi yet at ae 
penn a 
al — 
4 — 
a 
MAIN STREAM 7 Pt 
u a 
Ps 
/ 
JA WAKE 
& 
/ s/ 
/ \ 
Fic. 11 Shape of a jet from an orifice (b/a = 0; c/a = 0.504 


Probably no other solutions of the type discussed in 
this paper are obtainable in terms of elementary func- 
tions. Solutions for general values of a@ and 86 are 
probably available by numerical integrations but are of 
a difficult nature because of the presence of the unde 
termined constant /. As yet, the author is unaware of 
an effective procedure for the inclusion of the parame 
ters of variable viscosity, unequal densities, and jet 
velocities other than the particular values studied here. 


REFERENCES 


' Milne-Thomson, L. M., Theoretical Hydrodynamics, 2nd Ed., 
pp. 269 ff.; The Macmillan Company, London, 1934 

> Churchill, R. V., Complex Variables and Applications, \st 
Ed., pp. 171 ff.; McGraw-Hill Book Company, Inc., New York 
1948 

Callaghan, E. E., and Ruggeri, R. S., /nvestigation of the Pen 
tration of an Air Jet Directed Perpendicularly to an Air Stream, 
N.A.C.A. T.N. No. 1615, p. 6, June, 1948 











Re 
ous 
Tun 
viou: 
expe 
air ¢ 
only 
press 
num 
Nun 
with 
genc 
sure 
Mac 
(or | 
Dat 
nozz 
and 
thes 
vith 


A 


the 
rev: 
not 
cha 
Po 
ciel 
Will 
sur 
for 
pre 
and 
test 
) 
sigi 
resi 
pre 
sur 
nor 


R 


Bur 
Lat 
Arn 
rep: 
the 











O4 








An Experimental Study of a Hypersonic 


Wind-Tunnel Diffuser 


PETER P. WEGENER? ano R. KENNETH LOBBt 


LS. Naval Ordnance 


ABSTRACT 


Results of a variable-area diffuser investigation in the continu 
Naval Ordnance Laboratory Hypersonic 
A brief introduction discusses pre 

The diffuser investigated and the 
The results show that 


ous 12- by 12-cm. 
Tunnel No. 4 are presented. 

vious supersonic diffuser work 
experimental technique are then described 
air condensation in the test section at high Mach Numbers has 
only a minor effect on diffuser performance. Data on overall 
pressure ratios for starting and maintaining hypersonic flow for a 
The test Mach 
a diffuser 


number of diffuser configurations are presented. 
Numbers range from 5.9 to 9.6. From these data, 
with a single-peaked throat and a three-degree plane wall diver 
gence aft of the throat was selected as most practical. The pres 
sure recovered by this optimum diffuser varies, depending on 
Mach Number, from 1.8 to 2.3 times that recovered by a pitot 
(or impact) tube operated at the test-section Mach Number 
Data on pressure-distribution measurements throughout the 
nozzle and diffuser, spark schlieren photographs of diffuser flow, 
and data on Reynolds Number effects are also given. Finally, 
these data are compared with those of other investigators and 


with one-dimensional theory. 


INTRODUCTION 


A MAJOR COMPONENT OF supersonic wind tunnel 
working sections is the diffuser that decelerates 
the flow from supersonic to low subsonic speeds. A 
review of this problem is given by Ferri.' Diffusers do 
not operate free of losses, and their performance is 
characterized by indicating their overall pressure ratios 
(po/p., see Fig. 1) or pressure recovery. This is suffi- 
cient because the air leaving the diffuser has practically 
wind-tunnel supply temperature. A low overall pres- 
sure ratio or high pressure recovery is then desirable 
for wind-tunnel power-plant operation. This operating 
pressure ratio is different for different types of diffuser 
and test-section configurations. It is also a function of 
test-section Mach and Reynolds Numbers. 

Many supersonic wind tunnels have diffusers de 
signed primarily on the basis of subsonic experience, 
resulting in high overall pressure ratios.” Diffuser end 
pressure in these cases is below test-section pitot pres- 
sure. (A pitot tube decelerates the flow through a 


normal shock and a subsequent isentropic compression. ) 
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This lower recovery is due to the fact that, aside from a 
normal shock, viscosity effects in the diffuser introduce 
additional losses, thus lowering the end pressure. Since 
a diffuser employing a system of oblique shocks should 
have a better pressure recovery than one with a single 
normal shock, efforts were made to improve supersonic 
This method had al- 
Variable 


wind tunnels along these lines. 
ready proved successful in ram-jet diffusers. 
area diffusers whose throats can be closed after flow 
has been established were of interest here because of 
their higher pressure recovery. 

To the authors’ knowledge, the first supersonic 
tunnel diffusers giving end pressures higher than pitot 
pressure were discussed by Kurzweg* and Neumann and 
Lustwerk.* Diggins (N.O.L., unpublished) extended 
Kurzweg’s work up to a Mach Number of 4.9. The 
first hypersonic diffuser was investigated by Bertram.° 
His tests were carried out in the 11l-in. Hypersonic 
Wind Tunnel of the N.A.C.A.’s Langley Laboratory at 
a Mach Number of 6.9. End pressures up to 2.1 times 
pitot pressure were achieved by him. The present 
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Fic. 2. Working section of the N.O.L. Hypersonic Tunnel No. 4. 
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Fic. 4. Overall pressure ratios for three diffuser configurations 
investigation with a variable-area diffuser covers a 
Mach Number range from 5.9 to 9.6. 

The overall pressure ratio needed for starting super- 
sonic flow is important because it alone determines the 
maximum performance requirements for the power 
plant. At the present time little is known on this sub- 
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ject.°-* Starting requirements were also investigated 


and the results are presented. 

Finally, it is interesting to note that, although the 
open jet diffuser type is amenable to theoretical treat- 
ment,® the closed jet diffuser of interest here cannot yet 
be treated by a unified theory. This is because the 
variation of pressure and shearing stress along the 
solid walls must be included in the analysis. Such a 
calculation would require presently unavailable knowl- 
edge of turbulent boundary-layer characteristics in 
coverging and diverging channels and an understand 
ing of shock-wave boundary-layer interaction. 


WIND-TUNNEL AND EXPERIMENTAL METHOD 


The continuous N.O.L. 12- by 12-cm. hypersonic 
tunnel" operates from a 3,000 Ibs. per sq.in. supply of 
dried air. Pressure-regulating valves control the 
supply pressure of the tunnel in the range from | to 30 
atmospheres. Electric heaters control the supply tem 
perature of the tunnel in the range from room tempera- 
ture to above 500°C. 
air condensation in the test section is to be avoided. 


Heating of the air is needed if 


A steel wedge-type nozzle with a built-in cooling system 
to maintain nozzle wall temperature at a low value 
during high-temperature operation expands the air to 
a desired Mach Number in the range from 5 to 10 (see 
Fig. 2). The Mach Number distribution in a traverse 
across the exit plane of this wedge nozzle, exclusive of 
the boundary layer, is uniform to within +0.011/ at 
M = 7.2, the Mach Number for which most of the data 
are shown. 

Each diffuser wall has three plates, the first of which 
is linked directly to the nozzle end. Because of the 
high Mach Number in the test section, the shock angle 
caused by the flow deflection at this point is small, and 
long models can be used without having the conven- 
tional parallel wall test section. Three pairs of jacks 
control the motion of the three plates with respect to 
the working section wall. Sliding gaskets of Silicone 
rubber stripping fastened to the diffuser plates seal the 
diffuser in all positions. A flexible compartment seal 
between the first plates and the outside tunnel walls 
keeps the air behind these plates at near test-section 
pressure. The diffuser can be operated electrically 
during the run from fully opened to fully closed position. 
The diffuser throat opening is recorded to within 0.005 
in. Steel sidewalls with pressure taps along the center- 
line or walls with 1-in. thick circular commercial plate- 
glass windows enclose the diffuser. Photographs are 
taken with a spark light source of '/s-microsec. dura- 
tion. 

The diffuser leads into a 12-in. pipe connecting the 
tunnel to vacuum pumps. Since the high-pressure air 
supply lasts for hours and the vacuum pump capacity 
exceeds tunnel requirements, tunnel operation is con- 
tinuous. Arbitrarily chosen overall pressure ratios can 
be established to the order of 1 per cent accuracy. 
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HYPERSONIC WIND 

To determine operating conditions, supersonic flow 
was established with a high pressure ratio. (Pressure 
ratios of the order 10* are available.) The diffuser 
throat was then closed to a point where supersonic 
flow broke down in the test section. After this mini- 
mum throat area had been established, minimum over- 
all pressure ratios were measured for all diffuser throat 
areas larger than the minimum area. The overall pres- 
sure ratio (po/p,) at the moment of ‘“‘breakdown’’ of 
supersonic flow in the test section was noted. Flow 
breakdown was determined by any of three methods, 
which yield identical results: visual observation of the 
flow in the test section with the schlieren system, change 
of noise, and change of test section static and pitot 
pressure. 

Minimum starting diffuser throat areas and mini- 
mum overall starting pressure ratios were determined 


in a similar manner. 


OPERATING REQUIREMENTS 


If a supersonic tunnel is operated from a reservoir of 
dry air at room temperature and pressure at M/ > 4.8, 
a fraction of the air will condense at, or shortly after 
reaching, its saturation point in the nozzle. '! Such 
air condensation affects the commonly measured flow 
parameters differently. In particular, static pressure 
is extremely sensitive and pitot pressure nearly insensi- 
tive to the presence of condensed air in the flow. Fig. 3 
shows overall pressure ratios for one diffuser configura- 
In one test the air was not 
In the other 


tion and nozzle-area ratio. 
heated and air condensation occurred. 
test the air was preheated above the supply tempera- 
ture needed to keep the thermodynamic state of the air 
outside the condensation region during the entire ex- 
pansion process. (To minimize Reynolds Number ef- 
fects in this comparison, the supply density was kept 
the same for both the high and low 7.) It can be seen 
that overall pressure ratios in the case with air conden- 
sation are somewhat lower than in the heated case 
free of such condensation. This is to be expected since 
the actual flow Mach Number is higher in the ‘‘hot’’ 
case. However, the difference is small enough for 
engineering purposes to state that diffuser performance 
is nearly insensitive to air condensation as previously 
found for the performance of the pitot tube.” The 
following tests (except the sensitive static pressure 
measurements and checks on starting) were therefore 
made without preheating the air, and the Mach Num- 
bers indicated as M’ are to be taken as those derived 
from the insensitive pitot measurements with the aid 
of a flow table. This procedure (of running the tunnel 
with cold air) makes it possible to photograph and ob- 
serve the flow without the danger of cracking glass 
windows because of the heat. It also speeds up testing. 

The pressure recovery of the diffuser in terms of 
pitot pressure over the range of Mach Numbers tested 


was generally of the same order. A detailed investiga- 
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tion was therefore carried out at the optimum Mach 
Number 7.2 only. 

Overall pressure ratios as functions of diffuser-area 
ratio for three different diffuser configurations are 
given in Fig. 4. The Reynolds Number, based on 
tunnel width, for the comparable condition of no air 
condensation (Fig. 3) is 3,500,000. It can be noted that 
the case of parallel second diffuser plates® requires the 
lowest pressure ratio for large throat openings. No 
improvement is obtained when the angle between the 
second plates is increased from parallel to 1° to allow 
for boundary-layer growth. Diffuser performance is 
also practically unaffected if the angle between the 
third plates is increased from 3° to 14° and presumably 
flow separation takes place between the last plates. 
Fig. 5 gives a comparison of the pressure ratio for two 
different throat locations corresponding to the first and 
second jack position. It can be seen that the first 
throat position gives a better pressure recovery. It is 
assumed that in the second case the advantage of an 
increased number of shock reflections is more than off- 
set by the increased friction along the longer duct. 

A reduction of supply pressure increases the viscous 
effects at the diffuser throat. If the supply pressure is 
lowered, the minimum throat opening must be in 
creased to accommodate the thicker boundary layer as 
shown in Fig. 6. It was also found that for a given 
diffuser throat an increase in supply pressure resulted 
in a decreased overall pressure ratio. 

Fig. 7 gives a set of spark schlieren photographs taken 
with the configuration of a single peaked throat and 
subsequent 3° wall divergence. Pressure-ratio data 
corresponding to the conditions photographed in Fig. 7 


are shown in Figs. 3 and 4. 





108 JOURNAL OF THE AERONAUTICAL SCIENCES—FEBRUARY, 1953 

















30 T 

25r a 
= 
= 
<20-F _ 
° 
a 
WwW 
s 
oO 15r 4 
o 
uw 
er 
a 
> 
zx 10 be 4 
a 
= 
” 

M'=7.2, T9=!15°C 
5 4 7“ 
DIFFUSER THROAT MEASUREMENT (MM) 
8 10 12 14 16 18 
Oo 1 rt 4 1 | l 1 3 
0.06 0.08 0.10 0.12 0.14 0.16 


DIFFUSER AREA RATIO A/a. 


Fic. 6. Effect of supply pressure on minimum running area ratio 





(b) (c) 


Fic. 7. Spark schlieren photographs of diffuser flow; 17’ 
7.2, po = 30 Atm., To = 15°C. (a) Test section and converging 


diffuser plates, A;/A2 = 0.20 (optimum); (b) converging diffuser 
plates, A;/A2 = 0.60; (c) converging diffuser plates, A;/4. = 
0.40 


In the photographs it may be noted that the bound- 
ary layers along the test section and diffuser walls are 
turbulent. 
and temperature measurements in the layers with 7» 


This was confirmed by detailed pressure 


high enough to avoid air condensation. 
Finally, Fig. 8 gives pressure distribution measure- 
ments taken along the sidewall centerline. These 


measurements were carried out with the supply air 


heated to 7) = 320°, 30°C. higher than is needed to 
avoid air condensation throughout. 

Diffuser data for the same configuration taken at 
other Mach Numbers as check points are presented in 
the collection of all available results in Figs. 12 and 13. 

Using the data given above, we may postulate the 
following diffusion process: The oblique shocks emanat- 
ing from the junction of first diffuser plates with the 
nozzle exit and their single reflection slightly ahead of 
the throat cause the significant increase in pressure. 
The Mach Number at the throat outside the large area 
taken up by boundary layers is still rather high (in the 
case of J = 7.2 in the test section it is of the order 4 to 
5). Further diffusion from the throat to the end of the 
diffuser is relatively inefficient; in fact, it is about 
equal to that obtained in any diverging pipe with an 
entrance Mach Number of about that at the diffuser 
throat. Aft of the throat (see Fig. 8), transition to 
subsonic flow takes place through an oscillating shock 
system. Noise measurements at this point show a 
predominant frequency of the order of 5,000 cycles per 
sec. for a diffuser setting giving optimum pressure 
recovery. This frequency is a function of diffuser 
throat area, decreasing to about 1,000 cycles per sec. 
if the diffuser throat is wide open. Parallel second 
diffuser plates apparently contribute little to the stabi- 
lization of this shock system.!” 

Aside from structural reasons, it is desirable to 
eliminate the parallel section, because at MJ = | 
the product of density and velocity shows a maximum 
and, consequently, the heat transfer from the flow to the 
tunnel walls is maximized. This is important because 
of the high stagnation temperatures involved, result 
ing in the possible need for diffuser cooling systems in 
addition to nozzle cooling. This heat transfer was 
sufficient to cause glass cracks near the diffuser throat 
during heated runs. 

In a further series of tests with cone cylinder, sphere 
and missile models, and supports in the test section, it 
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was found that, aside from an increased minimum dif 
fuser throat area, overall pressure ratios were little 
changed.’ Apparently, the favorable additional ob- 
lique shock system introduced by the model counter- 
acted the detrimental viscous effects due to surface 


friction and wake. 


STARTING REQUIREMENTS 


Minimum pressure ratios and minimum area ratios 
to start the tunnel are shown in Figs. 9 and 10. It can 
be seen that, although the minimum starting area ds 
appreciably larger than the minimum operating area, 
it is about one-third smaller than that predicted by one- 
dimensional theory.'! If the diffuser throat is larger 
than the minimum needed to establish supersonic flow, 
the starting pressure ratios are nearly equal to those 
given previously for operation. The minimum starting 
pressure ratio is somewhat lower than the pitot pres- 
sure ratio for the test-section Mach Number, and it is 
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considerably below the pressure ratios previously antic- 
ipated by hypersonic wind-tunnel designers. This 
result its important because it may eliminate costly 
additions to hypersonic tunnel power plants for the 
momentary attainment of extremely high pressure 
An effect of the rapidity with which flow was 
A distinction was made 


ratios. 
established could not be found. 


“slow” and a “‘fast’’ start. In the case of the 


between a 
“slow’’ start, the overall pressure ratio was slowly 
increased, while the tunnel operated subsonically. 
The pressure ratio at which steady supersonic flow oc 
curred in the test section was noted. For the “‘fast’’ 
start, the starting pressure ratio was established within 
1/500 sec. by means of the fast-acting valve. No change 
in this ratio could be found with respect to the previous 
test. During the starting process, a supersonic jet 
detached itself from one nozzle wall and passed into the 
diffuser as shown in Fig. 11. Only for sufficient pres 
sure ratio did this jet spread to the full width of the 
nozzle. Neither the diffuser configuration nor the 
presence of models had an appreciable effect on start 


ing requirements. 


COMPARISONS AND CONCLUSIONS 


The present results are compared with data from 
other diffusers in Fig. 12. A general increase of pres- 
sure recovery with Mach Number in terms of pitot 
pressure can be noticed. The drop at the highest Mach 


(b) 





Fic. 11. Schlieren photographs of jet in test section before 
complete starting of tunnel; p) = 10 atm., JT) = 25°C. (a 
po/Psp = 60.8 and (b) po/pPs» = 62.8. (Starting pressure ratio 
po/Psp = 64.4; M’ = 7.2.) 
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tunnels. 


Numbers might be avoided by optimizing the length 
of the first diffuser plates. In fact, a special configura- 
tion would be needed for every Mach and Reynolds 
Number to attain the best performance. However, 
for practical reasons, a diffuser with a fixed first plate 
length is desirable. Its shape would then be selected 
to be reasonably efficient in a whole range of Mach 
Numbers and have maximum efficiency at either the 
highest Mach Number or at some other critical point 
dictated by the power-plant requirements. Although 
recovery appears high in terms of pitot pressure, it 
must be remembered that the latter drops to low values 
at high Mach Numbers (at W/ = 10, po’/po = 3/1,000). 
This is illustrated in Fig. 13 by a plot of diffuser effi- 
ciency as defined on the figure. Again the available 
results of other tunnels are given, and all data appear to 
lie on a single curve. 

No tunnel starting data suitable for direct comparison 
could be found in the literature. 

Finally, the following may be stated: 
single peaked throat, variable-area diffuser with plane 
walls was found to work as well or better than more 
complicated configurations in the range of Mach Num- 
bers 5.9 to 9.6. In this range of Mach Numbers the 
pressure recovery of this diffuser was better than that 
of a pitot tube operated at the same Mach Number. 
It was found that the essential part of the pressure 
recovery takes place up to the throat and that viscous 
effects play an important role in this process. Since 
neither the parallel plates nor the angle of the diverging 
section had any significant effect on the pressure re- 


A simple, 
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covery, it isconcluded that the diffuser section after the 
throat may be kept short. This minimizes the total] heat 
transfer to the wall for high-temperature operation oj 
hypersonic tunnels. It was also found that hypersonic 
tunnels may be started at pressure ratios about equal to 
the pitot pressure ratio, thus making additional start. 
ing devices unnecessary. The minimum diffuser throat 
area at which the shock system could be ‘‘swallowed’ 
and supersonic flow could be established was only two- 
thirds of that predicted by one-dimensional theory in 
spite of the viscous effects present. Apparently, an 
oblique shock system is immediately built up during 
the starting process. The presence of models and sup. 
ports increased the minimum diffuser throat somewhat 
but did not materially affect overall pressure ratios. 
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Analysis of Stiffened Curved Panels Under 


Shear and Compression 


M. A. MELCON* anp A. F. ENSRUDt 
Lockheed Aircraft Corporation 


SUMMARY 


This paper presents a method of analysis for curved panels rein 
forced by longitudinal stiffeners subjected to shear and compres 
sion. New formulas have been derived for the critical buckling 
stress of curved panels in shear and for the various effects of the 
diagonal tension field in the postbuckling state of the sheet 
These latter effects are expressed by fictitious compressive stresses 
in the stiffeners which are then combined by an interaction equa- 
tion with the compression stresses due to the external loading. 
In addition, a formula is given for the ultimate shear strength of 
the sheet. The method proposed in this paper has been com- 
pared with available experimental data, and satisfactory agree- 
ment is found. The method has been put in a form that requires 
the solution of certain mathematical relations and a minimum of 
chart reading and, hence, is readily adaptable to IBM or other 
high-speed computing techniques. A sample interaction diagram 
is included showing how the results of this method may be pre- 
sented for practical application. 


INTRODUCTION 


I SPITE OF GREAT ADVANCEMENTS in the design of air- 
craft structures, the analysis of stiffened shells still 
presents a challenging problem. The difficulty arises 
mainly from the fact that no practical mathematical 
treatment of the stress pattern in the postbuckling 
state has been developed as yet. Therefore, for the 
time being, the designer has to be satisfied with semi- 
empirical solutions that show sufficient agreement with 
test results. 

Based on a study of the available literature giving 
theoretical and experimental data pertaining to this sub- 
ject, a method has been developed for the prediction of 
the ultimate strength of sheet and longitudinal stiffeners 
in a curved panel. 


NOTATION 


Ay = stiffener area 

a = ring spacing 

b = stiffener spacing 

C,; = rivet factor, ratio of net to gross area of web 

F.' = allowable compression stress for stiffener alone; F,’ is 
the lower of either the column allowable (use a 
fixity of 2 for stiffeners continuous across rings) or 
the crippling cutoff of the stiffener 

F. = allowable compression stress for stiffener plus effective 


skin 


Meeting, 
Revised and 


Presented at the Structures Session, Twentieth 
LA.S., January 28-February 1, 1952, New York. 
received August 27, 1952. 

*Group Engineer, Structural Methods and Dynamies 


+ Structures Engineer. 


effective compression stress in stiffener due to diagonal 
tension 

a measure of stiffener bending stiffness required to 
break up shell into panels 

applied compression stress based on stiffener area plus 
effective width of skin 

applied compression stress based on stiffener area plus 
total area of skin 

Allowable gross area shear stress for web failure 

basic allowable gross area shear stress (for homogene- 
ous diagonal tension field at 45°) 

critical buckling shear stress for shell 

critical buckling shear stress for flat panel 

F,' + Fs," = critical buckling shear stress for curved 
panel 

ultimate shear stress of the material 

applied gross area shear stress 

part of applied shear stress carried in diagonal tension 

ultimate tensile stress of the material 

moment of inertia of stiffener about centroidal axis 
parallel to the tangent of the skin contour 

torsional stiffness factor: for open sections, J = 
A gtst?/3; for closed sections J = 4A %ts:/p, where A 
is the enclosed area and p is the perimeter 

shell parameter 

radius of curvature of panel 

thickness of web 

thickness of stiffener 

a factor reflecting the various effects of the diagonal 
tension field on the stiffener 

a factor indicating the intensity of diagonal tension in 
web 

a factor reducing the critical shear stress under com 
bined loading 
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FIGURE |. 
TYPICAL SHELL STRUCTURE 
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METHOD OF ANALYSIS 


A semimonocoque shell reinforced by rings and longi- 
tudinal stiffeners (Fig. 1) may fail in three ways: (1 
failure of stiffeners, (2) failure of sheet, and (3) failure o 
rings. If several components fail simultaneously, the 
type of failure is usually referred to as general insta- 
bility. It is assumed that the attachments between 
skin and supporting members will not shear or pop and 
thus reduce the ultimate load carrying capacity of the 
structure. 

The following deliberations are limited to an investi- 
gation of the first two types of failure. The results of 
this investigation are then checked against the test data 
of specimens that failed by buckling of stiffeners or 
tearing of the skin. Formulas derived in this paper 
apply specifically only to thin-walled shells of circular 
cross section with constant values for shear flow, stiff- 
ener area, stiffener spacing, and sheet thickness. Adap- 
tations for variation in these items are easily made. 

The proposed method may be outlined as follows: 

The critical buckling stress of a curved panel is given 
as the sum of the strength of the shell and that of the 
panel. If the applied shear stress is below the critical 
shell-buckling stress, no longitudinal stiffeners are re- 
quired. On the other hand, it is obvious that the stiff- 
eners must have a certain amount of bending stiffness 
in order to subdivide the shell into panels. As soon as 
the applied shear exceeds the critical buckling stress for 
a curved panel, a rather complex interaction between 
sheet and supporting members is started. The mathe- 
matical treatment becomes extremely involved, if not 


) 
f 


impossible. 

At the start, the shell between the rings bows outward 
like a barrel. With increasing shear flow, the radial 
components of the diagonal tension field will pull the 
stringers inward. In addition, a nonlinear relationship 
between the applied torque and the compressive stress 
in the stiffeners can be observed. Failure usually 
occurs from an individual buckle in the shell forcing the 
stiffener out of its original location, thus bringing about 
its collapse by an intricate beam-column action. 


Critical Shear Buckling Stress 

A simplified expression has been derived for the criti- 
cal buckling stress of a curved panel. For a flat panel 
the critical buckling stress may be expressed by 


F,” = KE(t/b)? (1) 


8 id 


where the magnitude of the coefficient A depends on 
the material, the degree of edge restraint, and the aspect 
ratio of the panel. For aluminum alloys and standard 
construction, this value is taken equal to 5.25 irrespec- 
tive of aspect ratio. The variation in panel length is 
taken care of in the term that reflects the effect of curva- 
ture. The critical buckling stress of the curved panel 
may be expressed by 


Fy... = KE(t/b)? 
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where in this case A = 5.25 + (m/4)(b/a)*S”* and 
S = a*/Rt denotes the shell parameter. 
for the critical buckling stress of a curved panel may be 


This formula 


transformed to the form 


te _ Fr’ cs F, ” (2 


where 


F,! = (4/AWS )(Et/R) 
FPF.” = 5.25E(t/b)? 


Postbuckling Behavior of a Shear Panel 


With the introduction of thin sheet construction for 
aircraft structure, engineers began to accept the new 
idea that the buckling of structural components did not 
necessarily indicate failure. Since that time, a great 
deal has been written on the postbuckling behavior of 
structures. Of foremost importance in the study of 
this subject, is the theory of the incomplete diagonal 
tension field. 

The stress pattern in a flat sheet subjected to shear 
forces beyond its buckling strength is known as di- 
agonal tension. A rigorous formulation of the transi- 
tion from the unbuckled state to Wagner's’ ideal di 
agonal tension field has not yet been accomplisiied. 
Semiempirical formulas developed by Kuhn® are most 
widely used. 

When the stress in a plane web starts to exceed its 
initial buckling strength, the applied shear forces are 
gradually taken by a combined truss action of the web 
and stiffeners. The sheet acts more and more like a 
diagonal, while the stiffeners take the place of uprights. 
There is a tendency for a buckle to form from corner to 
corner of the panel, provided this pattern is compatible 
with the deformation of the stiffeners supporting the 
edges of the panel. 

The various methods proposed for the analysis of flat 
panels in the postbuckling state differ in the assumption 
of magnitude of compression stress the sheet is able to 
Additional complications 

The diagonal tension in 


sustain in its buckled shape. 
arise when the panel is curved. 
the sheet tends to reduce its curvature in the direction 
of the wrinkles. This action induces nonuniform radial 
loads on the longitudinal stiffeners. 


Analysis of Longitudinal Stiffeners 

The longitudinal stiffeners in a reinforced shell per- 
form three functions: (1) They subdivide the shell 
into panels; (2) they sustain axial and radial loads in- 
duced by the tension field; and (3) they sustain directly 
applied axial loads. As long as the shell stiffened by 
rings only is able to sustain the applied shear flow with- 
out buckling, no stiffeners are required for the purpose 
of reducing the panel size. With increasing shear flow, 
it usually becomes more desirable to raise the initial 
buckling stress by the addition of longitudinal stiffeners 
than by thickening of the shell. The required bending 
stiffness of the stiffeners which will raise the initial 
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FIGURE 2. 














FIGURE 3. 











buckling stress to the full value of shell plus panel may 
be determined according to Seydel’s” formula for the 
critical shear stress of flat orthotropic plates, 


F,... = (32/ta’) V D,D»' 
where 
D, = D 
D, = EI,,/b 
D = Et*/12(1 — pn?) 


” and eliminating on 


Setting this expression equal to F, 
both sides, there follows the expression for the portion of 
the moment of inertia, A/,,, of the stiffener required to 


divide the shell into panels, 


Al,, = (bt?/5)(a/b)’ 


As mentioned previously, this requirement will be 
transformed into a fictitious compression stress indi- 
cating the portion of the stiffener strength needed for 
effective subdivision into panels. It is assumed that 
the ratio of the stiffener moment of inertia required for 
this purpose to its total moment of inertia may be taken 
equal to the ratio of an additional area to the total area 
of the stiffener. The following relation for the fic- 
titious compression stress /, 


f"A,, = F,'\A, 


” 


is established : 


st 


Then, 


nEI,,/a°A,, reflects the column strength of 
Combining and solving these equations 


where F,’ = 
the stiffener. 
gives 


f." = F."(at/A.,) V0.053b/a (3) 


This equation indicates that heavy skin and wide ring 
spacing require strong stiffeners in order to avoid exces- 
sive values of f.”._ However, this portion of the total 
effective compression stress in the stiffener in most 
structures is usually small. 

The effects of the diagonal tension field on the stif- 
eners in a buckled shell are rather complex. The axial 
load build-up in the longitudinal stiffeners caused by 


the diagonal tension (Fig. 2) is given by 
P= (f, — F,.. )bt cot a 


In addition, there are radial components of the diagonal 
tension field which produce bending moments in the 


stiffeners. The pull per running inch (Fig. 3) exerted 
by the tension field along this chord line is 


faitii@-f 


ji tan a 


and the radial load per running inch is 


P, = Tb/R = (f, — Fi.,.)bt(tan a/R) 


The unknown angle a may be found by trial and error. 
However, this standard approach to the analysis of 
longitudinal stiffeners does not agree with tests for 
many reasons, and, hence, a different approach to the 
problem was chosen. 

Fig. 4 pictures a stiffened shell subjected to torsion. 
Assuming a diagonal tension field of 45° and no bending 
in the stiffeners, the induced compression stress in the 


stiffeners is given by 


fe’ = (fy — Peer.) (0t/A st) (4) 
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FIGURE 4. 
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The actual angle of the tension field is usually smaller 
than 45°; in addition, there is a pull P along the direc- 
tion of the wrinkles which forces the stringers toward 
the center of the cylinder. It can be seen from Fig. 4 
that the pull P has a greater effect on stringers of smaller 
flexibility, because the shell support in this case starts 
at a later state. Furthermore, it is obvious that this 
effect is increased with sharper curvature. These, in 
addition to other considerations, led to the determina- 
tion of the empirical factor v, which was introduced into 
Eq. (4) to give a fictitious compression stress that is a 
measure of required stiffener strength to sustain the 
effects of diagonal tension. Eq. (4) then transforms 
into 

tf.) = W(f, — Face.) (0t/A,,) (5) 
where 


vy = 1+ (a/R)V (I,,/Jq,)(t/b) 


This formula gives satisfactory results for shells with 
R < 100in. The total effective or fictitious compres- 
sion stress for which the stiffener must be designed be- 
comes, then, 


f.’ +f." = total effective compression stress in stiff- 
ener 


The allowable stress to be used for determining the 
strength of the stiffener is F.’, which is the lower of the 
column allowable (using a fixity of 2 for stiffeners that 
are continuous across rings) or the crippling cutoff for 
the stiffener alone. Allowables determined from tests 
on the stiffener by itself may often be unsuitable for this 
analysis, since the stiffener may fail in a mode not possi- 
ble for a stiffener that is attached to the shell. Table A 
summarizes the effective compression stress to be used 
at different values of f, when the panel is subjected to 
shear alone. 


TABLE A 
f.’ and f-” When Panel Is Subjected to Shear Only 


hag y," 


f-’ = 0, when fs < Fer fe” = 0, when f, < F, 


fa — Fag, ot mn (!: - re) ( F, is 
af > Ae pie le = eg As 


when f, 2 F, 3 me 
; = 4 , when F,’ < fo < Fi 


a 


ail as F,"at 3/0.053b 
. Ax ie 


when f, 2 F, 


, 


If, in addition to shear, the panel is also subjected to 
direct compression, the effects of shear on the stiffener 
are combined by an interaction equation with the ef- 
fects of the direct compression. The only exception to 
the method as previously outlined for determining the 
effects of shear on the stiffener is that the critical shear 
buckling stresses, in this case, must be reduced because 
of the effects of compression. There are well-known 
interaction equations that give the initial buckling of 
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panels under combined loadings. However, for the 
purpose at hand, these equations are not satisfactory, 
since they will not properly account for the amount of 
the shear being carried in diagonal tension nor will the 
angle of diagonal tension be the same as when the sheet 
buckled under shear alone. 
conservative to assume that, after buckling occurs in 


In other words, it would be 


combined shear and compression, any additional shear 
is carried in the same manner as if the buckling were due 
to shear alone. If linear interaction were assumed for 
initial buckling in compression and shear, the following 
equation may be written: 


fF) + / Fe) = | 


Cer 

Then the shear stress at which the panel buckles be 
comes 

l , 

Ty = . . . F,.. 

i+ G/U]? a) 

Based on this reasoning, a factor \ was arbitrarily 
selected to reflect the reduction in the critical shear 
buckling stress due to compression, 


d= V1/[1 + (fe/f,)] (6) 


where /,, 


stiffener area plus total area of skin. 


is the applied compression stress based on 


The prediction of stiffener failure due to the com 
bined action of shear and direct compression is based 
upon the following interaction equation: 


de PWN 1.125 7 1.125 
(‘ tie) + (*:) =1 (7) 
c c 

where f, is the direct compression stress based on 
stiffener area plus effective width of skin and F, is the 
allowable compression stress based on the same area. 
In other words, the ratio f,/ F, is the ratio of the applied 
compression load to the allowable compression load of 
the stiffener plus skin. Table B summarizes the effec- 
tive compression stress to be used at different values of 
f, when the panel is subjected to shear and compression. 
For \ = 1, Table B is identical with Table A. 


TABLE B 
f-’ and f-” When Panel Is Subjected to Shear and Compression 
Se’ fe" 
fe’ = 0, when f, < AF, fe” = 0, when f, < AF,’ 


ee Wfs — AF; )bt (‘: = ak at) 
re = cr : _ - — , 
. Ase F, Ast 


when fs 7 AF, 310.053 


a 


,when AF,’ S f.< AF 


er 


oa AF,"at_ 3/0. OF3b 
adi 7 é f ot a : 


when f, 2 AFs,, 


Little test data could be found for shells reinforced by 
Whether the factor 
still applies in this case could therefore not be sub- 
The only 


stiffeners with closed section. 


stantiated by comparison with test data. 
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FIGURE 5 - ULTIMATE ALLOWABLE .GROSS AREA SHEAR STRESS. 


two specimens that had closed stiffeners failed by tear- 


ing of the web. 


Analysis of Shear Web 


The formulas for the allowable gross area shear stress 
presented here are based on the assumption that the 
net area of the sheet is equal to, or greater than, 75 per 
cent of its gross area. In the unbuckled condition of 
the web, the maximum tension stress in the sheet is 
equal to the shear stress. For aluminum alloys, the 
ultimate tensile strength is approximately 1.67 times 
the ultimate shear strength of the material. This in- 
dicates a marginal tension capacity of 67 per cent. In 
the ideal tension field, where the compression strength 
of the sheet is assumed equal to zero, the maximum ten- 
sion of the sheet is equal to twice the shear stress. This 
indicates that in the homogeneous tension field the sheet 
cannot develop the ultimate shear strength of the ma- 
The curve of Fig. 5, giving the basic allowable 
However, a 


terial. 
Shear stress F,,, takes account of this. 
correction is also needed to provide for an angle a dif- 
ferent from 45° and for wrinkles forming from corner to 
corner of the panel. This correction is reflected by the 
empirical factor 


7 = W11 + (a/b)?}/2 


ifa/b < 1, use 1; ifa/b>3,use3. The ultimate allow- 


able gross area web stress is, therefore, 


F, = F.. + ((Fe — Fee.)/a] (9) 


8 


where F,, is obtained from Fig. 5. 

The curves of Fig. 5 are similar to the lower curves of 
reference 11, Fig. 9. The curves of this reference were 
modified so that for f,/F;,. = 1 the basic allowable 


er. 


shear stress becomes F,, = 0.75F,,. This correction is 
made so that net area shear stress will never exceed 
| ee 


In exceptional cases where the net area of the sheet is 
less than 75 per cent of the gross area, an additional 
check is recommended. The circular interaction for 
combined tension and shear is given by the formula 


| fy | k fi, — Fre. | 
a sa san = j (10) 
CF, CrP in 


where C, denotes the rivet factor. 

If compression is superimposed on the shear stress in 
the panel, the initial buckling stress is again multiplied 
by the factor \ in both Eqs. (9) and (10). 


Analysis of Attachments 


At either end of the sheet or at splices, the load in the 
sheet must be transferred by rivets or fasteners to other 
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FIGURE 6 - TYPES OF SPECIMENS TESTED 


parts of the structure. The shear load for which each 


fastener should be checked is given by 


| (Uj, — F,...) }? 
P = qgp¥i + ke a | (11) 
\ i 


where p = pitch of fasteners and g = shear flow in 


sheet. Again, reduce F,,, by the factor \ if compres- 


r 


sion is present. 


COMPARISON WITH TEsT DATA 


The substantiation for the methods proposed herein 
are based upon the test work performed by the N.A.C.- 
A.*8, the Aluminum Company of America,'? Douglas 
Aircraft Company, Inc.,* and the National Bureau of 
Standards.’ A sum total of 54 test specimens has been 
Fig. 6 shows the general 
As can 


analyzed in this investigation. 
configurations of the types of specimens tested. 
be observed, some of the specimens are complete cylin- 
ders stiffened with longitudinal stiffeners and rings, 


while others are portions of cylinders. The radius of 


specimens varied from 15 to 60 in.; skin gage, 
from 0.020 to 0.051; and the factor v had a range 
on the specimens examined from approximately 1.10 
to 3. 


Critical Buckling Shear Stress 

Table 1 gives a comparison of the critical buckling 
shear stress computed by the method of this paper, by 
the method of reference 9, and observed test values. 
Twenty-six panels were checked. Test values averaged 
2.0 per cent higher than predicted by the proposed 
method and 5.8 per cent higher than by the method of 
Figure 7 shows graphically the relation- 
There is no 


reference 9. 
ship between test and predicted values. 
significant difference between the predictions of the two 
methods. The advantage of the proposed method lies 
in the fact that the buckling stress of the curved 
panel is expressed explicitly in terms of parameters 
of the panel and that reference to charts is not re- 


quired. 


Panels Subjected to Shear Only (Stiffener Failures) 


A total of 24 specimens that were subjected to shear 
only failed in the stiffener. Fig. 8 shows the relation- 
ship between the shear strength predicted by the 
method of this paper and test shear strength. Including 
all specimens, the average conservatism of the predic- 
tion is 7.8 per cent. The results from specimen 2 of 
reference 2 and specimen 10 of reference 4 were overly 


conservative. In the first case the stiffener apparently 
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. TABLE | 
Comparison of Initial Buckling Stress Computed by Fyer. = Fe’ + Fy” by Method of Reference 9 and Observed Test Values 


Panels Subjected to Shc 


Proposed method, 


Reference Specimen rau oe FF” + 
l 14 2,869 
20 2,844 
21 3,455 
2 2 3,018 
5 4,128 
Pe 4,854 
11 5,442 
3 4,245 
6 5,580 
9 5,755 
12 6,710 
3 l 3,318 
2 4,897 
3 3,643 
4 4,908 
5 6,317 
6 8,38 
7 «fe 
8 8,é 
9 3,§ 
10 2.4 
11 iP 
12 j 
4 8 
28 
27 
32 
16 
3 
10 
5 20a 
30a 
6 l 





Initial Buckling Stress 


Method of 


ar Only 


Ratio of Test to 
Predicted 


Observed Proposed Method of 


ref. 9 by test method ref. 9 
2,840 2,950 1.028 1.039 
2,750 2,870 1.009 1.044 
3,230 4,120 1.192 1.276 
3,160 3,460 1.146 1.095 
3,730 4,130 1.000 1.107 
+, 480 +, 660 0.960 1.040 
5,400 5,160 0.948 0.956 
+,790 4,720 1.12 0.985 
5,700 6,040 1.082 1.060 
5,480 5,840 1.015 1.066 
6,120 7,150 1.066 1.168 
2,950 3,200 0.964 1.085 
+, 580 5,000 1.021 1.092 
3,480 3,700 1.016 1.063 
1, 600 4,800 0.978 1.043 
6,700 6,400 1.013 0.955 
8,360 8,500 1.014 1.017 
6,530 6,700 0.942 1.026 
8,640 8,800 1.034 1.019 
3,940 3,040 0.774 0.773 
2,350 2,380 1.02: i.014 
7,110 7,900 1.046 1.111 
6,170 7,270 1.123 1.179 
1,470 

1 ,230 

1,460 

2,200 

1,230 

1,460 

2,200 

12,780 12,200* 0.846 0.955 
10,850 13,150 1.073 1.212 
3,280 3,700 1.083 1.128 


* This observed test value is probably low because of some initial condition, since identical panels with larger radii buckled at a 


higher stress. 


developed a fixity of 4 instead of 2 as determined from 
the compression stresses read in the strain gages, and in 
the second case the test is questioned, since this speci 
men was one of a family of specimens and the test re 
sults were incompatible with the rest. If these two 
specimens are excluded, the average conservatism be 


comes 5.8 per cent. 


Panels Subjected to Shear Only (No Stiffener Failures) 


As a negative check of the method for predicting 
stiffener failures, 15 specimens that did not fail in the 
stiffener were analyzed for stiffener failure. The pur- 
pose of this was to determine to what extent stiffener 
failures would be predicted by the proposed method 
when they did not actually occur in the test panel. 
Fig. 9 is a graph of predicted shear stress vs. test shear 
stress and indicates that for most of the specimens the 
test stress is lower than the predicted stress. This is as 
it should be, since these panels did not fail in the stiff- 
ener. 

Those panels that had a test stress close to the 
predicted stress based on stiffener failure failed by pop- 
ping of the rivets, and it is questionable whether the 
panels would have carried much more stress had the 
rivets not failed. 


Panels Subjected to Combined Shear and Compression 

(Stiffener Failures) 

A total of 15 specimens subjected to shear and com 
pression were analyzed. Fig. 10 shows the interaction 
between effective compression in stiffener due to shear 
and the direct compression for the test specimens ana 
lyzed. 

The 
the method of prediction of 5.0 per cent. 
from one specimen from reference 4 with a high conser 


results indicate an average conservatism for 
The test result 


vatism was considered to be questionable, since it was not 
compatible with other testsin the same seriesof specimens. 
If this specimen is neglected, the average conservatism 
drops to 3.8 per cent. The actual conservatism in an 
interaction diagram between shear flow in skin and 
compression load on stiffener would be less than the 
above value, since the conservatism enters into only 
that portion of the shear flow supported by the stiff- 


eners. 


Panels Subjected to Shear Only (Web Failures) 


Nine specimens failed in the web. Fig. 11 shows a 
graph of test stress vs. predicted stress for web failure. 
The results indicate that the method is conservative by 


an average of 0.3 per cent. There is reason to question 
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INITIAL BUCKLING BY PROPOSED METHOD - ksi 


FIGURE 7.- INITIAL SHEAR BUCKLING STRESS IN CURVED 
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PREDICTED SHEAR STRESS BASED ON STIFFENER FAILURE - ksi 


FIG 8 - SHEAR STRESS AT FAILURE vs PREDICTED SHEAR STRESS FOR 
PANELS WITH STIFFENER FAILURES. PANELS SUBJECTED TO 
SHEAR ONLY. 


the results from three of the specimens. Specimens 20a 
and 30a of reference 5 failed by tearing the sheet in the 
end panels and was precipitated by shearing of rivets 
connecting stringers to jig. In testing specimen 1 of 
reference 3, the loading jig bound, and the applied 
torque was actually about 10 per cent greater than the 
test gage indicated. If these three specimens were 
omitted, the average conservatism becomes 4.0 per 


cent. 


-FEBRUARY, 


1953 
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TYPICAL APPLICATION 


A typical application to which this method of analysis 
may be applied is the determination of allowable 
strengths for a fuselage shell structure. The following 
is an outline of the procedure for constructing an inter 
action curve of combined shear and compression on a 
curved panel. An example of such curves is shown in 
Fig. 12. 
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PREDICTED SHEAR STRESS BASED ON STIFFENER FAILURE ~- ks 


FIG 9 - SHEAR STRESS AT FAILURE vs PREDICTED SHEAR STRESS FOR 
PANELS WITHOUT STIFFENER FAILURES. PANELS SUBJECTED 
TO SHEAR ONLY 














FIGURE IO-INTERACTION OF COMPRESSION DUE TO SHEAR 
AND DIRECT COMPRESSION FOR PANELS WITH 
STIFFENER FAILURES. PANELS SUBJECTED TO 
SHEAR AND COMPRESSION. 
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ANALYSIS OF 


Interaction for Stiffener Failure 

(1) The ring spacing, stiffener spacing, stiffener sec- 
tion, skin gage, and the radius of the panel are known; 
hence, the factor vy may be calculated. 

(2) Determine F,’ and F,. These may be calculated 
by standard procedures. The allowable compression 
strength of a panel is often determined as a load per 
stiffener and skin rather than as a stress on stiffener plus 
effective skin. For the rest of this example, it will be 
assumed that P is the allowable compression load and p 
is the applied compression load. 

(3) Assume a value for p. 

(4) For the assumed value p, assume several values of 
f, and obtain corresponding values of \ from Eq. (6). 
\ involves f,,, Which is determined by p/(A,, + 0t).] 

(5) From Table B, determine f,’ and f,” for assumed 
values of f, and corresponding values X. 

(6) Determine value of f, which makes 


‘es +s re -(2) 
F,' : P 


(7) This value of f, times ¢ is then the shear flow that 
may be applied with the compression load p. Repeat 
this procedure for other values of p to complete the in- 
teraction for stiffener failure. 


Interaction of Web Failure 


The effect of compression on the ultimate shear 
strength of the web is small, since its only influence is to 
reduce the critical buckling shear stress, F;,., by the 
factor. The web failure curve with the effect of com- 
pression included may be obtained as follows: 

(1) Compute 7. 

(2) Perform steps 3 and 4 under section above. 

(3) Determine the value of f, such that f, = F,, 
where F, is given by Eq. (9). It should be noted that 
in these determinations F,,, must be multiplied by the 
factor \. If the rivet factor C, < 0.75, then an additional 
check must be made in accordance with Eq. (10). The 
value f, times ¢ is then the shear flow that will produce 
web failure. 

The envelope of the stiffener failure curve and web 
failure curve forms the complete interaction curve for 
the panel. In plotting the curves, it is recommended 
that p = O be the first point investigated, since in this 
case \ is always equal to one and the trial-and-error 
solution involved in steps 4, 5, and 6 of the section above 
iseliminated. Once the value of f, at p = 0 has been 
determined, it is possible to make a close first approxima- 
tion as to what f, will be at other values of p, and again 
labor of trial and error will be greatly reduced. 


CONCLUSIONS 


(1) Simplified formulas for the calculation of critical 
shear buckling stress of curved panels have been 
derived. Test values averaged 2.0 per cent higher than 


predicted values. 
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PREDICTEC ALLOWABLE SHEAR STRESS — ksi 


FIGURE |i. - SHEAR STRESS AT FAILURE vs PREDICTED ALLOWABLE 
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FIGURE 12.- EXAMPLE OF INTERACTION CURVES OF COMBINED 
SHEAR AND COMPRESSION ON CURVED PANELS. 


(2) A method has been determined for predicting the 
ultimate strength of a stiffened curved panel subjected 
to shear. The average prediction is approximately 6 
per cent conservative for longitudinal stiffener failure 
and 4 per cent conservative for web failure. 

(3) A method has been determined for predicting the 
ultimate strength of a stiffened curved panel subjected 
to combined shear and compression. The average pre 
diction is slightly less than 3.8 per cent conservative for 
longitudinal stiffener failure. 


(Concluded on page 126) 
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SUMMARY 


A reverse-flow theorem for compressible nonsteady flow, valid 
within the limits of linearized theory, is derived. This theorem 
gives a general class of relations between linearized solutions for 
lifting surfaces in direct and reverse flow. Based on the same con 
siderations used to establish the theorem, an adjoint variational 
principle, which may be useful in approximate solutions of non- 
steady lifting surface problems, is obtained. To illustrate the 
uses of the reverse-flow theorem, it is applied to the determination 
of relations between aerodynamic coefficients in direct and re- 
verse flow and to the obtaining of influence functions for total 
lift, pitching moment, and rolling moment for a wing oscillating 
with arbitrary motion and surface deformation, in terms of the 


pressure distributions for simpler cases in reverse flow. 


INTRODUCTION 


aww )W THEOREMS for various cases of steady 
flow, valid within the limits of linearized theory, 
have been presented by von Karman,' Hayes,’ the 
author,*? Munk,* Harmon,’ and Brown.§ The most 
general forms of these theorems have been given by 
Ursell and Ward’ and by the author® using different 
methods of attack. It has been found that the general 
form of the theorem may be extended to nonstationary 
flows as well. The proof of this result for incompres- 
sible flow has been given previously, but the more 
general proof for compressible flow, both subsonic and 
supersonic, was left for the present paper. The deriva- 
tion is completely analogous to that given for nonsteady 
incompressible flow, but due account must be taken of 
the fact that the velocity potential for compressible 
nonsteady flow satisfies a form of the wave equation 
rather than Laplace’s equation. For convenience, the 
case of harmonic time variation will be considered. 
Other cases of variation in time follow from the result 
for harmonic motion by a direct application of the 
Fourier integral theorem.® Since the results obtained 
will be valid for all frequencies for which the linearized- 
flow equations themselves are valid, it follows that the 
results will also be valid for any time dependent motion 
of the airfoil which possesses a Fourier integral trans- 
form in the time variable. The unit step and unit im- 
pulse motions of the airfoil are, of course, among such 
motions, as are most others of practical interest. 

Received February, 1952. Revised and received October 23, 
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The results obtained will be seen to bear a certain 
formal relation to Helmholtz’s reciprocal theorems jn 
the theory of acoustics.!* In the case where the airfoil 
is at rest with respect to the fluid, the classical Helm 
holtz theorems are, of course, directly applicable, but 


this case is of little interest in aerodynamics. 


MATHEMATICAL FORMULATION 


A nonstationary compressible inviscid flow to the 
usual linear approximation is characterized by the well- 
known wave equation in a moving coordinate system,'! 


ea (5 +u a) = () l 
ao” \Ot — ” Ox an 


where ¢ is the usual velocity potential, such that the 
perturbation velocity components along the x, y, and z 


axis (u, v, and w, respectively) are given by 
Op Op Od 
= —, t= —., w= 
Ox Oy Oz 


dy is the velocity of sound in the free stream, and Lis 
the free-stream velocity along the x axis. The axes are 
shown in Fig. 1. The flow is irrotational and, in par- 
ticular, Ou/Oz — Ow/Ox = (0. The pressure at any 
point (to within an arbitrary additive function of time 
is, according to the linearized Bernoulli equation, 


x = P/po = Uou + (0¢/0t) (2 


where pp is the free-stream density and P is the pressure 
at any point. It may readily be seen that u = (0@ Ox 
and 0¢/0f satisfy Eq. (1) by simply differentiating Eq, 
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Wake of direct flow 


Wake of reverse flow 


1) with respect to either x or ¢. Thus, z also satisfies 


Eq. (1). 
The specified motion of the boundary of the lifting 
surface lying approximately in the plane z = O is given 


for harmonic oscillation by 


g = 2(x, y, t) = 20(x, ye (3) 


The boundary condition for the flow is 


Ody O20 P 
. = Wy = Us + iw (4) 
Oz Ox 


where @ go(x, y, se’ and w = w(x, y, Z)iwt. It is, of 
course, necessary to admit complex values of ¢@» and wo 
to take care of the phase differences that may exist. 
Over that part of the plane z = 0 occupied by the lifting 
surface, there is a pressure discontinuity, Az, between 
the upper and lower surfaces, but everywhere else in 
the fluid Aw must be zero. At the trailing edge, the 
Kutta condition requires that Aw = QO, in order that x 
go smoothly into its value in the wake. 

Now if there be considered a flow over the same plan 
form with the direction of the free-stream velocity Ly 
reversed, oscillating at the same frequency, there will 
be associated with this flow functions, @ and #, which 
have properties and imposed conditions similar to those 


Consider the volume integral 


- a * o” w Oo w? 
I, = / / / (ve - M? = See fe 
2 a ao Ox ay Ox ay 


where 7, the volume of integration, is taken to be all space outside the wing. 
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of @ and zw. For the reverse flow case, however, the 
wake trails from what was in the first case the leading 
edge, the Kutta conditions being applied there. It 
should be noted in each case that, across the plane z = 0, 
A¢ and Au differ from zero only on the wing and on the 
respective wakes. The wakes in the two cases are il 
lustrated in Fig. 2. In each case, there is a singularity 
in uw and x of the type 1 /+/d at the leading edge, d being 
the distance from the edge.*~* 

If the directions of positive x, y, and s and positive 
perturbation velocity 4, 3, @ are held the same in the 
reverse flow case for convenience, the equation satisfied 
by the velocity potential ¢ is 


4 _— —_ — 0 = ») 
? Ao” Ot Ox ¢ 
The pressure is given by 
t = P py = —Ui + (06/2f) .. (6 


® also satisfies Eq. (5) 
The expressions for and mw may be specialized to 
harmonic motion to give 
r= Lu + two (7 
t= —U ti + iwd (8) 
Similarly, Eqs. (1) and (5) may be written for the case 


of harmonic motion 


> 


0’ _@ Od w” 
Vv" — AP — M +4 =@Q (9 
? Ox? ay Ox Ay" 
0°o _w .. 0o 3 
V% — M*—+2—M + p= 0 (10 
- Ox ay Ox Ao” ¥ 
where J = l/do is the free-stream Mach Number. It 


is of interest to note that the two partial differential 
equations are adjoint in the usual sense, thus further 
illustrating the fact, already brought out for the integral 
equations of lifting surfaces in steady flow,* that re 


verse flows satisfy adjoint equations. 


It should be noted that Eqs. (1) and 


5) are satisfied by ¢ and ¢ everywhere outside the lifting surface except in the respective wakes and that, as v ap 


proaches + © in the wake of the direct flow, ¢ and # go to zero, while as x approaches — ~ in the wake of the re 


verse flow, @ and z go to zero. 


cause of the Kutta condition in the reverse flow and, similarly, for 7 and ¢. 
theorem may be applied to the volume integral, using as the fluid boundary, the plane z = 0, from x 


x= 4t+o., 


Further, at the leading edge, where the derivatives of ¢ are singular, 7 is zero be 


Thus, a simple extension of Green's 


© to 
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It will be convenient to consider the terms of the in- 
tegral J) in pairs of corresponding terms from the first 


and second parentheses. The first pair of terms gives 


n= ff f we — pV*x)dv 
» 


Green's theorem may be applied directly to /; to give 


oo SS or 
I, = 1S — 1S (13) 
SS ro ‘ 


where the area of integration A is the entire plane z = 0 
—o tox = +o and from y= —o to y= 


(12) 


from x = 
+o. The subscript + denotes the positive side of the 
plane, while the subscript —denotes the negative side. 
The n direction is positive for the outward normal 
of the bounding plane. Because of the antisymmetry of 
the lift problem, 


(%*) _ Op - (s*) “ op 
On/ + oz’ \On/— — oz (14) 
= =F). % = —F. 
Thus, 
25 
n= f fz? as f fis = dS (15) 


0 need 
now be considered, but P and 7 may be taken, in surface 
integrals over the plane z = 


where only values on the + side of the plane z = 
0, to represent the difference 


in pressure between upper and lower surface. 
The second pair of terms from the two parentheses in 


- wf ff (* es oo) dv (16) 
, Ox? 


This may be written 


- Oo Or 
a M a = = PT 
SSS Ox (* Ox , 4 . 


Carrying out the integration with respect to x gives a 


I) gives 


te 


(17) 


surface integral, 


M? -w ff (*° ot) ) dS 18 
I, a = 1 COS “X) ds 8) 
on ri) mys COS (N°X) ¢ ( 


The direction cosine of the angle between the m direc- 
tion (in this case, the z 
I = 0. 

The third such pair of terms in Jo gives 


— 24— uf ff (e+ 027) ae (19) 
x Ae 


which may be written 


direction) and x is zero, so that 


I; -_ 


i = 


.w ‘ }O 
— 2 . M JSS (7%) dv (20) 
ao V Ox 


Performing the integration on x gives 
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.w J 
21 M / J (1) cos (n-x) dS (2] 
ao A 


= 0 as for Jo. 
The remaining pair of terms in Jp gives 


ay~ V 


which is identically zero. 
Ip is equal to the sum of J,, Js, 73, and J4. 
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I; = 


which is zero, since cos (n-x) 


It has been 
shown that J», 73, and J, are zero, so that J>5 = J,... J, 
however, is clearly also zero, since ¢ and 7 satisfy wave 
equations of the forms Eqs. (1) and (5), respectively. 


Thus, /; 
re) 
® 1S _ ff "dS = (23 


= OQor 


SLE 


To put this result in a more symmetrical form, it is ex- 


pedient to write, from Eq. (8), 


OT — Ou goa Od 04 
Oz Oz Oz 
From the irrotationality condition, 
Ou/Os = OW/Ox 
Also, by definition, 
w = 06/02; w = 0¢d/O0z 
Then, 
O7/0z = —U((O®/Ox) + tw (25 


Substituting this result into Eq. (23) gives 


JS rwdS — SS o(-u 0 a a iw) dS =0 


(26 
The first term of the second integral may be integrated 


by ere on x, and since at x = —o, @ = OQ), while at 


= w@ = 0, there results 


uf fe Y ae ws —_ l fi / J 
x A 
2 wdS = — ur f | umdS (27 
Ox A 


Substituting Eq. (27) into Eq. (26), 


SS mw dS — J / (Usu + twh)@dS (28 
A or ae 


But according to Eq. (7), the term in parentheses in the 


second integral is inate 7, so that 
SS Fwds — SS zadS = 0 
or 


SS twdS = Sf Sf rods (29 


where the integrations are over the wing plan form only, 
since m and 7 are zero over the remainder of the plane 


z= 0. 
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This is the desired reverse-flow relation for com- 
pressible nonsteady flow. In the case of steady flow 
w = 0), Eq. (29) reduces to the theorems of Ursell and 


Ward’ and the author,* 


SS uwe=SS iw (30) 


THE VARIATIONAL PRINCIPLE 


As in the case of steady flow, the foregoing considera- 
tions may also be applied to the derivation of an adjoint 
variational principle, which provides a convenient basis 
for the approximate solution of lifting-surface problems. 
The type of variational problem considered is one in 
which solutions satisfying the partial differential equa- 
tion exactly but not the boundary conditions are 
employed. This is the common case in lifting-surface 
problems and is typified by the methods that involved 
building up a solution from distributions of sources, 
doublets, and elementary horseshoe vortices. Thus, 


consider 


my 7 og 24) 
bH = 6 w in) — 7 S 
I ff (war + ww) 5 (3 an +x d 


(31) 


where 6 is the usual symbol for variation according to 
the calculus of variations and w and @ are the specified 
arbitrary boundary conditions for 0¢/On and 0¢@/0n, re- 
spectively, on the wing plan form. The two problems, 
one in direct flow and the other in reverse flow, are in- 
dependent, and, accordingly, variations with respect to 
io and 6¢ or 67 and 67 are independent. The variations 
are arbitrary over the class of flows that satisfy the 
Kutta condition for direct and reverse flows, give finite 
downwash over the airfoil, and satisfy the differential 
equations, Eqs. (1) and (5). The variation of H7 may be 


written 
bH = JS (wit + Wor) — 
A 
l oo oo  _ 5g 258) 
j 5 - S (32 
5 (or 3° + ‘at oe . — 


Returning to the derivation of Eq. (29), it will be ob- 
served that the only requirements on ¢, 7, ¢, and 7 re 
lated to their satisfying the appropriate wave equation 
lor moving coordinates; no requirements as to bound- 
ary condition were imposed. Thus, Eq. (29) may be 
applied to the last two terms of Eq. (32) to give 


* { _ 06 “[ _ 6 
/ J T ? aS = / / on ° as 
A on JA On 


(3:3 ) 
“( dé¢ a 
If T ° aS = SS OT ® is 
‘ A On A on 
Substituting Eq. (33) into (32), 
, s dg 0 
oH -[f (w = *) om + (2 — *) br |dS (34) 
A On On 
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In the vicinity of the correct solution, ¢ and @ corre- 
sponding to w and @, respectively, 


w = 0¢/0n; @® = 0¢6/On (35) 
so that 
mas ff 
A 
Op oe 
(we -+- De) — = (= — 7 = ds =0_ (36) 
2 on On 
which is the variational theorem for this case. The 


variational method for approximate solution of non- 
steady lifting-surface problems may be applied in the 
same manner as for steady-flow problems, by applica- 
tion of a Rayleigh-Ritz type of procedure. This has 
been more completely discussed in a previous paper on 
the steady-flow case.* Further extensions of these ideas 
have been developed; it is possible to choose the varia- 
tion functions # corresponding approximately to such 
@’s as will make the total lift, total pitching moment, 
total rolling moment, ete., approximately stationary. 
The details of such procedures will be discussed in a 
later paper. 

It is also possible to derive a variational principle in 
which 5¢, 62, 66, and 6% do not satisfy the partial differ- 
ential equations but rather, the boundary conditions. 
At present, these appear to be of less interest in lifting- 


surface problems. 


APPLICATIONS OF THE REVERSE-FLOW THEOREM 


Applications of the reverse-flow theorem for non 
stationary flow are similar to those for stationary flow 
which have been given in detail in several papers. First, 
they may be used to establish relations between the 
aerodynamic coefficients for a wing in direct and reverse 
flow. These are of interest in the case where one of the 
solutions is already known or more easily obtained than 
the other. Second, and perhaps more important, in- 
fluence functions for total lift, pitching moment, and 
rolling moment are easily obtained in terms of the pres- 
sure distributions over the wing in reverse flow with 
certain simple prescribed boundary conditions. These 
are analogous to the generalized Munk formulas already 
given for steady flow;* they are so called because, in the 
limiting case of two-dimensional incompressible flow, 
they reduce to Munk’s integral formulas for the total 
lift and moment on an airfoil with an arbitrary distribu- 
tion of local angle of attack. The cases of two-dimen- 
sional flow, lifting-line approximations, or wings having 
a spanwise axis of symmetry are of particular interest, 
since, in these cases, a reversal of flow direction gives a 
flow around essentially the same wing. Thus, it is not 
necessary to consider a separate reverse-flow problem at 
all. A few specific applications will be illustrated below. 
These will serve to illustrate the manner in which the 


general theorem may be applied. 
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(1) Relations Between Lifts Due to Vertical Motion 


Let 


ai/c = (h/c)e™) 


‘ (37) 

Zi/c = (h/c)e™S 
corresponding to vertical motion of the surface, where c 
is any representative dimension of the plan form. By 
Eq. (4), 


w, = iwe(h/c)e™| 


@, = iwc(h/c)e"™4 


Applying the general theorem, Eq. (29), it may be 
noted that since w and @ are equal and constant, they 
may be taken out from under the integrals and can- 
celed, leaving the integrals of the pressure over the sur- 
Thus, for vertical motion, the lifts in direct 
In all that follows, lifts up- 


face equal. 
and reverse flow are equal. 
ward are defined as positive. 


2) Relation Between Lifts and Pitching Moments for 
Pitching Rotation 
For the case when the wing is rotating about the axis 
v = 0 (actually, the choice of origin is arbitrary) in both 
the direct and reverse flows, a more complex relation is 


obtained. Let 


Ze = xe'™ \ 20 
. ‘ P (ed! 
w. = (Uy + iwx)e'f 
and let 
Z = —xe | 10 
2 eee ; oP (40) 
wW= [—l 0(OZo Ox) -+- tw30 Je’ “ 
OT 
@, = (Uo — iwx)e™ (41) 


Applying the general theorem, 


f [2ls-s2s- 
Uoc 
c e ~weox 
J fe. (: +i = *) dS 
Uo Cc 


The integrals of P, and P; over the surface are simply 


the lifts L. and /»2, respectively, and, also, 


S SPA -x/oOdS = 
SS PAx/odS = 
where a positive moment is defined as nose-up for either 
Finally, then, 


Ms/c\ 
ms > (42) 
WU, ¢ { 


direct or reverse flow. 


~ wc Me rs . we Me 
In +4 =In+1— (+5) 
l 0 C l 0 ¢C 

If, for instance, the axis of x is along the quarter-chord 
position for a rectangular wing in direct flow, it is along 
the three-quarter-chord position in the reverse flow, and 
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Eq. (42) provides a relation between the lifts ang 
moments for oscillating rotations around these positions. 


(3) Cross Relation Between Lifts and Pitching Moment for 
Pitching Rotation and Vertical Translation 


Let the two cases considered be w; and P, of | and 
@, and P,» of 2. 
h/c = 1, 


/ / PI - ‘) dS = / / Pyi - dS 
0 € ( 


or using the previous definition of lift and moment, 


Applying the general theorem for 6, = 


we M, 06 = 
LL; +1— =%4— le 14 
l 0 € Uo 

where L, is the lift and 1/7, the moment about the axis 
x = 0 for vertical translation of the wing in direct flow 
and where J» is the lift of the wing rotating about the 
axis x = Qin reverse flow. For a two-dimensional wing 
or one having a spanwise axis of symmetry (such as a 
rectangular wing), this provides a relation between the 
lifts and moments of the wing for two different motions 
with the same flow direction. 


4) Relation Between Rolling Moments of Wings 
Oscillating in Roll 


Let 
z; = aye) 15 
- { ” 
z; = aye 
" tw 
W3 = waye 
7 “ twl 
@2. = tuaye 
Applying the general theorem, 
JS JS Psy dS = JS Psy dS 1) 


But J JS Psy dSis the rolling moment in direct flow and 
SS Psy dS is the rolling moment in reverse flow, so that 
the rolling moments due to oscillation in roll are equal. 


(5) Influence Function for Lift 
The total lift for a plan form executing any oscillatory 
motion whatever may easily be obtained in terms of the 
pressure distribution of the plan form in reverse flow 
executing vertical translational oscillations only. 
Let 


(47 


where again c is some representative dimension and 


z/c = (20o/c)(x, ye | 


w = U (020/0x) + iwe(zo/c) 4 


(48 


To ®, corresponds the presumably known pressure dis- 
tribution P,, while to w corresponds the unknown dis- 
tribution P. 
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Applying the general theorem, 


Sf Gr oes)nas- f frw 

— = n 

Ox Uo c ™ Uo , 

, f fr 1S l ly IS (= ‘ 1w0e =) P 1S 

s = “ = — —— e ( ) 
J joe dx ' Ue) ' ae 


gives the total lift for the wing in direct flow with arbitrary 20. 
It is of interest to compare this last result with a formula given by Schwartz’? for the total lift on a two-dimen- 
In the present notation, Schwartz'* gives for the pressure 


then, 


sional oscillating airfoil oscillating with arbitrary w(x). 
distribution on a two-dimensional wing oscillating in vertical translation in incompressible flow 


- = = we\? , 
P, = pUi’e™ | (1 + 7 : 2 ( ) \ ie =| 5()) 
1 plore | a U, \; — + U, l x (50 


where 7 is a function of we/ Uy only. 
For a reversal of flow direction, this becomes 


= oe al i+x ma. z 
P, = pUvte™| (1+ 7 2 ( ) Vv — | (51) 
4 . | Uso ¥; _ = + Uo ’ 


Substituting this result into Eq. (49), there is obtained, for the lift on an airfoil with arbitrary w(x, t) = wo(x ei" 


L= Uy? | (1 1 2 V1 — x?| dx 
u. ? | + oper U, | v? | ax 


F aQ+T) | . Wwe, ‘ 
L = Qal), / w(x, | - 1. + x 4 - Vi — “| dex (52) 
2 l--<s 0 


This is just the result obtained by Schwartz as a result of rather lengthy calculation. 


or 


6) Influence Function for Pitching Moment M i. a a 
, , eta ee =- Pr, — a5 (O4) 
An influence function for the pitching moment on a r 10 . Uo 


wing performing oscillations with arbitrary w(x, y) : ; : ee 
: : which gives the desired result for pitching moment 


may also be obtained in terms of a combination of sim- : ai a ; : 
: A ys about the axis x = 0. The formula given by Schwartz 
ple solutions for the wing in reverse flow. Let ; ape ; : ike 7 
for pitching moment in the two-dimensional incom- 
= —xe™ + (Uy/iwe™ (53) pressible case has been checked by Eq. (54). 
Then 7) Influence Function for Rolling Moment 
@ = —Ue™ + ioxe™ + Ue rhe influence function for the total rolling moment 
on a wing oscillating with arbitrary w(x, y) is also 
sal readily obtained in terms of the pressure distribution on 
. ns Sle the wing in reverse flow oscillating in roll. Let 
@ = —wxe 
= lwl : 
; a at Zz = ye } 
It may be seen from Eq. (53) that the pressure distribu- os (55) 
tion P corresponds to the superposition of the reverse- a a \ 
ie _ Ww 
flow solutions for the wing oscillating in rotation about 
the axis x = 0 and for the wing translating vertically with pressure distribution P, presumed known, and 
. ° r . . ” ‘. wl 4 ae ow . : . 
with amplitude U/)/iw. Applying the general theorem, W = Wo(%, ye with pressure distribution P presumed 


unknown. 


i If P. (- $ aS « [ fr 4 dS Applying the general theorem, 
bi , ™ Iw ; ft 
: Sf PydS = / / P—ds 
Uo : ae Uo 


or 
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or leading to this result, an’ adjoint variational principle, 


, IJ PydS = - v7 PwdS (96) 
C « TW 


The integrai on the left side is the rolling moment that 
gives the required result. 


CONCLUSIONS 


A general relation between linearized solutions of 
lifting-surface problems in direct and reverse flow has 


been established for compressible nonsteady flows. This 
relation is a direct extension of that already known for 


steady-flow solutions. On the basis of the analysis 
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Differential Shear-Flow Analysis for 
Nonprismatic Semimonocoque Beams 


WILLIAM B. KING* ann WALTER R. GARRISON? 


Boeing Airplane Company 


SUMMARY 


This paper gives the derivation of a method for obtaining the 
shear-flow distribution in a nonprismatic semimonocoque beam 
with nonlinear axial-load variations due to externally applied 
loads. This method is based on the actual rates of change of 
segment axial loads at a particular station rather than the aver- 
age rate of change of segment axial loads between two reason- 
ably separated stations. The particular application presented 
is for a wing; however, the method is applicable to any semi- 
monocoque structure. It also has special adaptations for many 
loading conditions at each station. 

The shear-flow distribution may be determined for each com 
ponent of the loading condition such as torsion and beam and 
chord shear and The of the 
may be determined if the shear flow due to pure torsion is de- 


moment. location shear center 


sired. 

This method of analysis is adaptable to International Business 
Machine operations. The IBM calculations may be performed 
on machines such as the relatively slow 602 A Multiplier or the 
faster Card Program Sequence Calculator. The time required 
by IBM to calculate and check the shear center and unit shear 
flows and determine, for ten complete loading conditions, the re- 
sultant shears for a 20-segment station is approximately 30 min. 

The results obtained by using the method presented in this 
paper on two beams, one with straight surfaces and one with 
curved surfaces, are compared with those obtained by VQ/J and 
the ‘Unit Method of Beam Analysis.”’! 


SYMBOLS 


a = segment area 

1 = total section area 

c.g = center of gravity 

C = | (Ixlz Iyz*) 

C; = Clyz 

t. ad Clx 

C; = Cz 

E = modulus of elasticity 

f = axial stress due to bending 

G = shear modulus 

hx = chordwise coordinate with respect to Lr.a. 

hz = beamwise coordinate with respect to L.r.a. 

Ty, Iz moments of inertia about Y and Z axes, respec- 
tively 

Ixz = product of inertia 

Lr.a. = load reference axis 

Al = periphery of the skin or spar web between two 
adjacent stiffeners 

Vx, Mz = bending moments about the XY and Z axes, respec- 
tively 

n = segment number (counting clockwise ) 

N = total number of segments in a particular section 


Received March 26, 1952. Revised and received September 
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a = segment axial load 

AP = difference in segment axial loads at two adjacent 
sections 

q = corrected shear flow 

qo = correcting shear flow 

qu = uncorrected shear flow 

Q = first moment of area 

s = twice the sweep area enclosed by the skin or spar 
web € between two adjacent stiffeners and the 
section c.g. 

S.c. = shear center 

Ss = twice the section enclosed area = Dx 

t = skin or spar web thickness 

Ti.r.2 = torsion moment about Y axis through the |.r.a 

Ts.c. = torsion moment about Y axis through the s.c 

Vx, Vz = shear loads parallel to Y and Z axes, respectively 

w.b.p. = wing beam plane 

w.c.p = wing chord plane 

x = segment c.g. chordwise coordinate with respect 
to the section c.g 

Xs.c = s.c. chordwise coordinate with respect to the sec- 
tion c.g. 

¥ = section c.g. chordwise coordinate with respect to 
w.b.p. 

X = horizontal axis (chordwise ) 

y = section spanwise coordinate with respect to the 
wing tip 

Ay = spanwise distance bet ween two stations 

y = horizontal axis (spanwise ) 

Z = segment c.g. beamwise coordinate with respect 
to the section c.g 

o¢ = skin € beamwise coordinate at segment m with 

R respect to the section c.g. 

Sa.¢ = s.c. beamwise coordinate with respect to the sec 
tion c.g. 

z = section c.g. beamwise coordinate with respect to 
W.c.p. 

Z = vertical axis (beamwise ) 

6 = shearing strain 

‘Os = d/dy 

2 = should equal 


Subscripts Vx, Vz, Mx, Mz, and T refer to segment axial loads 
or shear flow, which are results of the indicated loading condi 


tion. 


INTRODUCTION 


— CONFIGURATIONS OF AIRPLANES designed today 
are vastly different from the configurations of those 
designed a few yearsago. The presently designed semi- 
monocoque types of airplanes have developed from 
those with wings of constant thickness ratios with un- 
machined stringers and skin to those with wings of 
rarying thickness ratios, machine-tapered skins and 
stringers, and tapered stringer spacing. Combinations 


97 








128 JOURNAL OF THE 


























RIGHT HAND 
hy RULE 











Fic. 1. Sign convention and coordinate system (loads and 
directions positive as shown). 


of these variables result in nonlinear segment axial- 
load variations and consequently ever changing shear- 
flow distributions. It has been realized for some time 
that the assumptions upon which the classical pris- 
matic beam theory (¢ = VQ/J) is based are seldom 
fulfilled. A second method, ‘“‘Unit Method of Beam 
Analysis,’’ was developed to account for taper in the 
beam, thus allowing a more accurate determination of 
the shear flow. This method is based on the difference 
in axial loads at two adjacent stations, and the accu- 
racy of the resulting shear flow depends upon how 
closely the segment axial loads between these two 
stations approach a linear variation. Should a wing 
with a relatively long span be designed in such a man- 
ner that the segmeut axial-load variations are non- 
linear, an excessively large number of sections must 
be analyzed in order to obtain sufficient accuracy for a 
shear-flow analysis. 

The method of shear-flow analysis developed in this 
paper is based on the rate of change of segment axial 
load at a particular station and accounts for the actual 
rates of change of segment areas, segment coordinates, 
and section properties. The accuracy of this method is 
dependent upon the accuracy of the means by which the 
axial stress distribution is determined. For a beam 
with no abrupt changes in continuity or unusual end- 
fixity conditions, the normal procedure of determining 
elastic bending stresses will be sufficiently accurate. 
Should a beam have abrupt changes (such as sweep- 
back, cutouts, or unusual end-fixity conditions), this 
method or any method of analysis will be accurate only 
to the extent of the capabilities of the engineer to de- 
termine the effectiveness of the structural material. 
Abrupt changes may be accounted for in the analysis 
by using effective segment areas in place of the actual 
segment areas. 

The assumptions upon which this method is based 
are similar to those used in most practical methods of 
shear-flow analysis. The material is assumed to obey 
Hookes’ Law, the moments of inertia of the individual 
segments about their own c.g.’s are assumed negligible, 
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and the segment areas are assumed concentrated at 
their centers of gravity. The load derivative is taken 
with respect to the Y component of the segment load. 

An interesting comparison is made in an example 
problem worked by two presently accepted methods 
of shear-flow analysis and the method presented by 
this paper. The two methods for comparison are the 
VQ/I method and the ‘‘Unit Method of Beam An. 
alysis.”’ 

The basic principle of the following development js 
simple. Primarily, it consists of setting up the equa- 
tion for the segment axial load and obtaining its rate 
of change with respect to the span by taking the first 
At a first glance it may 
Actually, 
the man-hours expended by engineering personnel in 


derivative of that equation. 
appear to be a long and laborious method. 


using this method may be short in view of the fact that 
it is adaptable to IBM operations. If the basic section 
properties used in the bending stress analysis are calcu- 
lated by IBM, the only additional items required to 
obtain the shear-flow analysis are A//t, (a)’, (hy)’ and 
(hz)’. The values of A//t are required only if the shear 
center is desired. In all widely used methods of shear- 
flow analysis which determine the shear center, the 
values of A//t are required. The remainder of the 
additional items will be discussed in later paragraphs. 


EQUATIONS FOR THE RATE OF CHANGE OF SEGMENT 
AXIAL LOAD 


Let axes X and Z of Fig. | be any two perpendicular 
axes through the section c.g. 

With the sign convention shown in Fig. 1, the elas 
tic bending stress at any point whose coordinates with 


respect to the c.g. are x and 21s 
IyMz\x + UnxzMz — IzMx|z 
IxIz — Ixz? 


[IxzMx 


This equation may be rewritten as 


f = Mx[Cox — Caz] + Mz[Cz — Cox ? 
where 
C, = 1/[Ixlz — Ixz? 3 
Cy on Clyz } 
C3 = Cx 3) 
Ci ” Cl, (6 


The segment axial load then is 


P = fa Mya[Cox — Cyz|] + Mzal| Coz C3x] (4 


The rate of change of P with respect to y is obtained 


by taking the derivative of Eq. (7) with respect to J, 


(P)’ = Mya(Cox — Cys)’ + [Cox — Cy2] (Alxa)’ + 
Mza(Coz — C3x)’ + [Cos — C3x] (Meza 
= Mya[Co(x)’ + x(Cy)’ — Cy(2z)’ — 2(Cy)’] + 
[Cox — Cyz] la( Wx)’ + Mx(a)’] + 
Mza[Cx(z)’ + 2(C2)’ — C3(x)’ x(C3)"] + 
[Coz — Cex] [(a(Mz)’ + Mz(a)’| (8) 
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SHEAR-FLOW ANALYSIS 


By combining terms, Eq. (8) results in 


(P)' = [a(x)’ + x(a)’]* [CQ:Mx — CsMz]? + 
[a(z)’ + 2(a)’]* (C.Mz — CiMx]® + 
[ax]* [C2( Mx)’ — C3(Mz)’ + 
Mx(C2)’ — M;,(C;)']® + 
[az]* [C2(Mz)’ — Ci(Mx)’ + 
Mz(C2)’ — Mx(C,)']® (9) 


It should be noted that in Eq. (9) the first term (*) of 
each of the four multiplications is a function of the seg- 
ment, whereas the second term (@) is a function of the 
section. 

The (P)’ derived is the change in shear flow between 
adjacent segments. The value obtained for a segment 
from Eq. (9) indicates an instantaneous change in load 
at segment c.g. The change is actually a gradual 
change between the c.g.’s of adjacent segments if the 
skin is effective in bending. 

Often it is desirable to know the breakdown of a par- 
ticular shear value—i.e., the values of the component 
parts due to each of the beam and chord shears and 
moments, and torsion. Therefore, the development of 
this method of analysis will deal with shear flows due 
to each of the above components. 

Eq. (9) may be broken down into four equations, each 
dealing with a (P)’ due to a single loading condition 


of shear or moment. Since 
(Mx)’ = V; 
x) al (10) 
(Wz) = | 2 


the following may be written: 


Py,)’ = aVx|Coz _— C3x | (333 
(Pyz)’ = aVz[Cw — Cis] (12) 


Puy )’ a(x)’ + x(a)’] CM, — 
la zs)’ + 3 a)'] CiMx +- ax My(C.)’ _ az.Mx(C,)’ 


(15) 


Pus) = la(z)’ 4. z(a)’ | C.M;z _ 
la(x)’ oo x(a)’] C,M, +. azM,(C2)' _ axM,(C3)’ 
(14) 


Eqs. (11), (12), (13), and (14) deal with the (P)’ 
It may be noted that three 
(a) area, (b) coordinates, and 


of an individual segment. 
basic items are involved : 
The section properties, of course, are func- 
The area is 


¢) loads. 
tions of segment areas and coordinates. 
assumed to be concentrated at the c.g. of the combined 
stringer area and the area of the skin acting with that 
stringer in bending. The most convenient way of ob- 
taining (a)’ is by plotting area vs. y for each segment 
and by graphically determining the slopes at any sec- 
tion. Before obtaining the rates of change of the 
coordinates of a segment and the section properties of a 
section, the coordinates must be defined. In the de- 
rivations to follow, a wing will be used as an example. 
The equations may be modified to suit any beam, such 


as a fuselage. 
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Fic. 3. Rear view, left-hand wing 


To obtain the correct relationship between section 
properties, such as (x)’, (z)’, and (a)’, at different sec- 
tions all coordinates should be referred to three basic 
perpendicular planes. Figs. 2 and 3 illustrate the co- 
ordinate system assumed. The horizontal plane will 
be taken as the w.c.p. The spanwise vertical plane 
will be taken as the w.b.p., which is perpendicular to 
the w.c.p. and through the Lr.a. 

The L.r.a. shown in Fig. 2 is the axis about which the 
The bending sections are taken 


torsions are computed. 
Normally, the load axis is 


perpendicular to this axis. 
chosen as one that is midway between the spars and in 
the w.c.p. The chordwise vertical plane will be taken 
at the tip perpendicular to the L.r.a. 

To correspond with the coordinate planes set up in 
Figs. 2 and 3, properties that increase positively going 
inboard will have a positive rate of change with respect 
to y. 

From Fig. 2 the horizontal coordinate of a segment 
may be expressed as 

x= hy — Tr (15) 


Thus, it can be written that 
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Fic. 4 Example Problem No. 1: sketch, areas, coordinates, slopes, and I6ads. 
(x)’ = (hx)’ — (%)’ (16) ([z)’ = 230 [ex(x)’] + DS) [x? (a)’] (23) 
iia Ix)" = 22, [a2(z)’] + >o[z? (a)"] . (24) 
(xz) = do [ax(z)’] + do [az(x)’] + do [xz(a)’] (25) 
tr = >o[ahy]/A (17) ~ ; ; oF ald 
rhe rates of change of the section coefficients C,, C, 
Therefore, C3, and C; are then determined. From Eqs. (3), (4), 
(x)’ = (hx)’ — R (3), and (6), 
4 [Scathy)’ + Shx(a)’] — [Sahx] (A) 18) Cc) UxIz — Ixz?)' 
= (I (Ci) = - 
A? fel ~ Ice?’ 


Similarly, from Fig. 3, ()’ may be written as 


(2) = (hz)’ i 
A [Svathz)’ + Sohz(a)’] — [Sahz] (A)’ 


A* 


(19) 


Spanwise plots of hy and hz for each segment may be 
used to determine (hx)’ and (hz)’ at any section. 


For most airplane wings, the moments of inertia of 


the individual segments about their own c.g.’s are com- 
paratively negligible. In this paper they are neglected. 


On that basis, the moments of inertia and the product of 


inertia are 


Iz = (ax?) (20) 
Ix = (az?) (21) 
lez = >> (axz) (22) 


and their rates of change are 


Iy(Iz)’ 4. I,(Ix)’ Bi 2Iyz (Ixz)’ 


[Ixlz — Ixz*]? 
(C2)’ = C\(Ixz)’ + Ix2(Ci)’ (97) 
(C3)’ = Ci(Ix)’ + Ix(Ci)’ (28) 
(Cs)’ = Cz)! + I2(G)’ (29) 


The above derivatives may now be used in Eqs. (11), 
(12), (13), and (14) to determine (P)’ 


EQUATIONS FOR SECTION SHEAR-FLOW DISTRIBUTION 


In view of the fact that the shear flows are to be 
computed for each of the loading conditions (shear, 
moment, and torsion) it will be necessary to determine 
The s.c. is immaterial as far 
however, it is the only 


the location of the s.c. 
as the final net g is concerned; 
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Fic. 5. Example Problem No. 2: sketch, areas, coordinates, slopes, and loads. 


point about which the true torsion’ and, consequently, 
dr may be computed. Should the engineer desire to 
have a net shear-flow distribution rather than a break- 
down as indicated above, the location of the s.c. need 
not be determined; any reference point (such as the 
cg.) may be used, and the equations of equilibrium 
should be written about this point. 

Having obtained (P)’ for each segment of a section 
for each of the shear and moment loads, the section 
Shear-flow distributions for each of the loading condi- 


tions may then be obtained. 


Shear-Flow Distribution Due to Vz 


Assume the shear flow at any point in the periphery 
of the section to be zero. For example, assume the 
rear spar of the structure to be cut so that it can trans- 
fer no shear between the rear spar upper and lower 
chords. With this assumption, start at the rear spar 
just below the cut and sum, in a clockwise direction, 
the (P)’ The resulting uncorrected 


shear flow is 


qu = >, (Py:z)’ [Ref. 
It will be clockwise at any point at which q, is positive 
and counterclockwise at any point at which q, is neg- 
ative. Asa numerical check, after summing clockwise, 
q. Should be equal to zero at the cut. The sum of the 


for each segment. 


Eq. (12)] (30) 


forces in the vertical direction should equal the applied 
Vz. The sum of the horizontal forces should equal 
zero. Thatis, 


> [q. As] £ Vz (31) 

DY Iq. Av] £ 0 (32 
where 

Az Sn41 — Zn (33) 

AS @ 2454 — Xe (34) 


The total shearing strain around the section perimeter 
due to q, is 
5 ~ X[q,Al/t] (35) 
If the section is to deflect without twisting this strain 
must be balanced by an applied correcting shear flow 
go, where 


go = — X[quAl/t}/Z[Al/t] (36) 
The resulting corrected shear flow is 
Qvz = Qu + GW (37) 


Since gyz deflects the section without twisting, the 
resultant must pass through the s.c. For tke hori- 
zontal location of the s.c., take moments about the Y 


axis through the c.g. This may be written as 
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TABLE | 
Example Problem No. 1: Station ‘‘A’”’ Using Differential Shear Flow Analysis 


4 
/IN. | #/IN. 


* REF. FIG. 4 
4 SEE DETAILS BELOW 


© *-M, 9/1, *-£54/16219=-@/3_ €Q.(1) 


WHERE My=54 IN.@ ; Iy=Z0a¥*7=162 IN* EQ(21) 


®@ =Ca(y! +%4(0)'IE-My/1,1=-CO+@® 1/3 


@ *LayIe(My)/Iy+My(Iy)/IxI 
=(@) 1£38/10533 
WHERE (My)'=V,= 2# EQ.(10) 
(Iy)'=2E00.9(9)'I+ELY7(a)'3 
=220@@3+F0GI © 


=306/13 IN? EQ. (24) 
VeXs0. = —GS — Z[qQus] (38) 
Therefore 
Xsc. = —[qoS + =[qus]]/ Vz (39) 
where 
S = 2enXnt1 — Xn2%n+1 (40) 


Shear-Flow Distribution Due to Vx 


The derivation of the shear-flow distribution due to 
\’x follows the line of reasoning employed for the shear- 


flow distribution due to Vz. The resulting equations 


are 
Q, = =(Prx)’ [Ref. Eq. (11)] (41) 
[qulv 2 0 (42) 
Z[quAx] = Vx (43) 
VIquAz] £0 (44) 
do = —Z[q,Al/t]/Z[Al/t] (45) 
Qvx = qu + Go (46) 


Zc. = [goS + Z[qus]]/ Vx (47) 


CA(P)'3 (Pp) Hn. In Furi 


16 17 


A @+@ 6 @G 2 


W/IN. | #/IN. IN? 


4 2i 


(5) =LINE £(@) 


© = Gat Yor @ +5/78 
WHERE Go" EL-q,.01/S 
= 5C-(19)(@ 1/2) 
=5/78#/IN. EQ.(59) 


VERTICAL SHEAR CHECK 


XC-P(V)'I+EC 


-2@©@ +=C 


ayaz vy, 
= -@,1L®) zie 
36/13 +£-135/1053IC 3 +3722 
36/13 -10/132 2 
26A3=2 





Shear-Flow Distribution Due to Mz 

The uncorrected shear flow due to A/ is similar to 
thatdueto l. Thatis, 

du = 2(Paz)’ [Ref. Eq. (14)] (48) 
[du lx : 0] (49) 
The equilibrium numerical checks at this point for .\/, 
Should the 
segments not be perpendicular to the plane of the sec- 


however, are different from those for V. 


tion, the segment axial loads due to bending moments 
will produce X and Z components. Therefore, as a 


check, one may use the following equations: 


= [duAr] : Y[Puz(hx)'] (50 
Z[quAz] = =[Puz‘hz)’] (51) 


The correcting shear flow is determined by equating 
the moment about the Y axis of qg,, the X and Z com- 


ponents of Pz, and go to zero. That is, 


> [Puz2(hx)’] —_ > (Paz x(hz)"| “ Z[qus] = goS = () (52 
Therefore, 

Y[Pwz2(hx)'] _ Y[Puzx(hz)’] a wl dus ne 

go = 5 (ed) 
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TABLE 2 
Example Problem No. 1: Station ‘“B’’ Using VQ//J Shear Flow Analysis 
n|asly|a3 ¢ P (3)'} @ Fu X |Fn%ne:|A%nFnei| + a 
CoL.|] | J2]{ 3 4 5 6 | 7 » 4 10 i i2 13 
rmepi + |? OMO@! a |O®@ | # [e® A $ @n@Ond®n@nu|O-O A 
UNITS! IN2| IN.| IN? | #/IN2 + IN3 #/IN. IN. in? Inz@| Ine #/IN. 
| 13721-4|-6 | 67/92 | 201/184 |-2/13| -6 | -1035/711I0032| il | -22 -44 | 22 |5865/220064 
2 15/21-4| -10 | 67/92 | 335/184 |- 2/13| -16 | -2760/110032| 11/2 0 -22 | 22 | 24157220064 
3, |7/2|-4| -14 | 67/92 | 469/184 |-2713| -30| -51 757110032 O 22 0 22 |-24157220064 
4 1572|-4| -10 | 67/92 | 335/184 |-2/713| —40| -69007110032 |-l172 | 44 22 | 22 |+5865/220064 
5 |3/2|-4| -6 | 67/92 | 201/184 |-2/13| -46| -7935/110032 |-I! 44 -44 | 88 |+-79357220064) 
6 13/2] 4] 6 |-67/92 -201/184| 2/13|-40/-6900/110032 |-II -22 -44 | 22 |+5865/220064 
7 15/72] 4] 10 |-67/92 +335/184| 2/13|-30/-51 75/110032 |-11/2 i} -22 | 22 |-2415/220064 
@ |7/2| 4| 14 |-67/92 |-469/184| 2/13|-16 |-27607110032| O 22 (e} 22 | 24157220064 
9 15/2] 4| 10 |-67/92 -335/184| 2713/- 6 [- 10357110032 | 1172 44 22 | 22 |5865/220064 
10 |3/2| 4] 6 |-67/92 -201/184| 2/i3| O 0 TT 44 | -44 | 88 | 79357220064 
= [23 ene 352 
¢ REF. FIG. 4 ( = aut qo*@ +7935/220064 
WHERE qo = 2 A/S 
4 SEE DETAILS BELOW = 3£-(@) @ 3/2@ 
*698280/(550163(352] 
@ *-My Y/ Ly *-C 67/3681 F=-67/368@) E£Q.(1) * 7935/220064 #/IN. 
WHERE M,=54+2CI3/21]= 67 IN. # 
I,*30a%7*1=368 IN* EQ(21) VERTICAL SHEAR CHECK 
*-VQ/I, =345Q/[5981[ 368] = 345 (® /220064#/IN EC-P(BI+ZCqmAYI= vz, €Q.(3!) & (58) 
WHERE V=V, - 3C-P(9)'] 1541/598- (7935I(83/11I003222 
=2- 1541/598 =-345/5988 2=2 
The resulting corrected shear flow is Ts. = Tire + V2zl% + X%c.]| — VxlZ + 2...) (61) 


+ do (54) 


QMz = Qu 


Shear-Flow Distribution Due to Mx 


The equations for shear-flow distribution due to /x 


are sinilar to those for \/7. They are as follows: 


The torsional shear flow is a constant value for the en- 


tire periphery of the section and is 


dr = 7. c S 


(62) 


Net Shear-Flow Distribution 


qs = > (Pux)’ [Ref. Eq. (13) ] (55) The net shear-flow distribution is the algebraic sum 
Iqulx + O (56) of the individual shear-flow distributions due to each 
DV lquAxv] £ SO[Pux(hx)’] (57) of the beam and chord shear and moment, and torsion. 
[g,Az] 2 > (Pw (hz)’] (58) 
dig = d= Qvz2 + Qvx + uz + Ix + 97 (63) 
qo = 
> [Pury2(hx) 4 ey > [Pury hz) ‘| oo > (qus) 50) Unit Shear-Flow Distributions 
(5! ne , 
a Since the design conditions on a wing are often nu- 
Qux = Qu + Qo (60) merous, it is desirable to compute unit shear-flow dis- 


Shear-Flow Distribution Due to T..- 


Normally, torsions are computed about Lr.a. as a 
matter of convenience. To obtain the torsion about 
the s.c., the. beam and chord shears must be trans- 
ferred to the s.c. and the resulting torsion added to 
T\,. For this transfer the following equation may 
be written: 


tributions for each section analyzed rather than to 
carry the complete analysis through for each of the de- 
sign conditions. The magnitude of the unit values, 
Vz, Vx, Mz, Mx, and 7,,.. will depend upon the size 
of the structure analyzed. The shear-flow distribu- 
tions for any particular design condition may then be 
obtained by multiplying the unit shear flows by their 
appropriate shears, moments, and torsion. 
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TABLE 3 


Example Problem No. 1: Station “gt Using Unit Method of Beam Analysis 


CH 4|QA\7X| F 
6 
+ 


5 


¢ REF. FIG.4 


4 SEE DETAILS BELOW 


@=-Mx9/Ix *- @/3 Gi) =LINE 29 
WHERE My=54 IN# 
Iy*ZCaP1=i62IN° 3 = gat Go 
= B +7/182 
WHERE 


©-O© 


9 [lojiifi2| 13 
BO} #/ti#] o 


@ =-My3/Ix *-4@ /35 
WHERE M,=80 IN# 
Iy=20L04"I=7001N* 


> 


(WF nXnui|Xn a | ap 
i4/is|] 6 


STATION 


VERTICAL SHEAR CHECK: 


IC-R(D)1+zCG.AVIEZ Ve 
£609 (9 343{E) CO, 3} 2 
54/9142 VERT. SHEAR DOES 
NOT GHECK 


Yo* -ILGus/ La 


@-@-© 


EXAMPLE PROBLEMS 


Discussion 


Two example shear-flow problems are shown. For 
each problem the shear-flow values are computed at 
three adjacent stations by the) Differential method, by 
the Unit Beam method, and by VQ/J. It should be 
noted that for reasons of simplicity the problems 
chosen are tapered beams with rectangular cross sec- 
tions. The example problems are presented merely to 
indicate that a difference can exist in the calculated 
values for shear flows obtained by the several methods 
noted above. The problems are not intended to be 
representative of any actual design. To illustrate the 
effects of surface curvature, the only difference in the 
two problems is in a side view. Problem No. | has 
straight upper and lower surfaces, whereas Problem 
No. 2 has curved upper and lower surfaces. The 
shear-flow values have been computed using fractions 
in order to show the checks on the summation of ver- 
tical forces. The shear-flow values shown in 
Table 4 are given in decimal form for comparative pur- 
poses. 


= 7/182#/IN 





Normally, for a wing, one chooses stations to check 
bending stresses and shear flows that are spaced span- 
wise approximately the same distance as exists be 
tween spars. In the example problems, however, sta- 
tions “‘A’’ and “‘C”’ are spaced slightly less than two- 
thirds of the interspar distance. An additional inter- 
mediate station is also included to show that the Unit 
Beam method, when calculated on the basis of the dif- 
ference in end loads between stations ‘‘A’’ and ‘‘C”’ but 
applied at a station midway between the two, will give 
results that closer approximate the more accurate 
values obtained at that mid-station by the Differential 
method. 
additional information needed to obtain good results 


It may appear at first glance that the only 


with the Unit Beam method are the sweep areas, 5, 
at a station midway between the two for which the 
segment axial loads are computed. This is true in the 
example problems because the sections are symmetrical 
about the X and Z axes. In actual design this is seldom 
true; therefore, the section properties, bending loads, 
etc., must be calculated in order to determine what 
portions of the external shear and torsion are carried 
by the slopes of the segments. 
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TABLE 4 


Comparison of Shear Flo ws for the Example Probl ems 


EXAMPLE PROBLEM NO. | 


STATION A 


DIFFEREN- 
TIAL SHEAR 


DIF FEREN- 
TIAL 


vQ/I 
SHEAR 
F 
FIG. 4 

/ 
.0495 


.04 


2.5 # 2.00 # 2.00 # 


EXAMPLE 
STATION 


STATION A 


TIAL SHR 
FLOW 


SHEAR SHEAR 
FLOW aF 


* 
.04 
0275 


- .0385 
-.04 
- O2T 


5 |-.1667 
-.46% | 2.00# 


CONCLUSIONS 


The selection of a method for determining shear 
flows of sufficient accuracy for practical design of a 
semimonocoque beam depends upon a number of vari- 
ables. No attempt is made in this paper to describe 
the limits of the parameters which determine whether 
or not the Unit Beam method, or any other similar 
method, gives sufficiently accurate results. The Unit 
Beam method, considering the classical bending stress 
assumptions noted in the introduction, gives accurate 
results if the segment axial load variations are linear. 
Under the same assumptions, the Differential method 
gives accurate shear-flow results for beams with linear 
or nonlinear variations in segment axial load. The 
QI method gives accurate results only for prismatic 
beams. 

In the design of the wing for a large intercontinental 
bomber now in production, a comparison of the shear- 
flow values obtained by the Unit Beam method and 
the Differential method was made. For the same num- 


STATION 





STATION © 
VQ/TI 
SHEAR SHEAR 
FL FLOW 
FIGS. 4&7 .4 


VQ/I 


SHEAR SHEAR 


67 ‘ 0495 .0165 


67 


1.96 # 2.00 # 
PROBLEM, NO.2 


2.00 # 1.69 # 


STATION 


SHEAR 
FLOW 


SHEAR IAL 
FLOW FLOW 


SHE AK 
FLOW 


IN. 
5 


2.00 # 2.92 # 2.00% 


2.00% 


2.41% 


ber of stations analyzed, spar shear differences of as 
much as 30 per cent were noted. The wing configur- 
ation was such that, had the Unit Beam method alone 
been considered, errors of at least one gage in shear 
material would have resulted. Accuracy obtained by 
the Unit Beam method comparable to that obtained 
by the Differential method would have required that 
approximately twice as many stations be analyzed. 
With the use of IBM calculating machines, the engi- 
neering man-hours required would have been greater 
had the Unit Beam method been used. 


The Differential method provides an expedient tool 
in the preparation of a unit-load analysis, which is 
desirable if a large number of loading conditions re- 


quire investigation. 
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Some Remarks on ‘‘A Solution of the Laminar 
Boundary-Layer Equations for a Compressible 
Fluid with Variable Properties, Including 
Dissociation’’! 


M. Z. Krzywoblocki 
Professor of Aerodynamics, University of Illinois, Urb&na, III. 
October 6, 1952 


— UNDERSTAND THOROUGHLY the interesting and highly 
elaborated results presented by Moore,! it seems that some 
remarks might be in order. These are: (a) As the fundamental 
equations, the Navier-Stokes equations are assumed. This 
means that the analysis refers to that domain of flow in which 
those equations are valid. As is well known, there is today a 
highly controversial question up to which altitude those equations 
are valid. Although the calculations were performed up to = 
20, they may be not necessarily valid in the régime of rarefied 
and highly rarefied gases for the simple reason that no one knows, 
as yet, which fundamental equations should be used correctly in 
those régimes. Certainly, the resuits may be correct in the hyper- 
sonic flow in a gas of a higher density. (b) The author ignored 
the second equation of motion (momentum). This means that 
the Prandtl hypothesis was applied—i.e., the boundary layer is 
extremely thin, and the vertical velocity component and the 
pressure gradient in the vertical direction can be neglected. As is 
well known, in the hypersonic flow in a rarefied gas the boundary 
layer is thick, the vertical velocity component and the pressure 
gradient in the vertical direction cannot be neglected, and the 
second equation of motion must be preserved. This is another 
hint that the calculation refers rather to a hypersonic flow in a 
gas of higher density. (c) The third remark refers to purely for- 
mal transformation of a system of partial differential equations 
into a system of ordinary differential equations. In general, such 
a transformation is impossible in a one-to-one correspondence. 
It was shown in the past by the author of the present note that 
in the transformation from x, y, to x, z, the Jacobian of the trans- 
formation vanishes, which clearly means that the transformation 
cannot be one-to-one. Thus, automatically the uniqueness of 
the solution cannot be proved. Although the applied transforma- 
tion may give a solution, it cannot give the unique solution. 
The next remark refers to the neglecting of terms containing x 
Again, in general, a priori, one is not allowed to neglect the terms 
containing x, since they may be of large magnitude. One should 
not diminish the degree of determination of the class of functions, 
which are supposed to give a solution of the differential system 
in question, by neglecting a priori some terms in the original non- 
linear differential systems. 

If one neglects them a priori, one should show a posteriori that 
his step was justified. In his work Mr. Moore assumes a priori 
that the velocity and temperature pattern depend only upon g. 
This may be a sufficiently strong point that the terms containing 
x may be neglected, although no proof is presented. But the de- 
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pendence on z only seems to indicate strongly that the calcula 
tion refers to a particular case—namely, to a parabolic distriby 
tion of velocity and all the other functions. The question arises 
by itself: ‘‘Which are those cases (altitudes, density relations, 
etc.) in which one is allowed to assume such a pattern?” 


REFERENCE 
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Some Remarks on the Effect of Flapping Hinge 
Offset on Rotor Blade Stall 


John R. Meyer, Jr., and G. Falabella, Jr. 
Massachusetts Institute of Technology, Cambridge, Mass 
October 7, 1952 


N A READERS’ FORUM ITEM in the JOURNAL OF THE AERONAU- 
I TICAL SCIENCES, July, 1952, Chawla presented a preliminary 
investigation of the effect of introducing a flapping hinge offset 
on rotor blade stall. It was concluded that, in the helicopter 
régime of flight and with normal values of blade mass constant 
(y = 12), the stall condition at the retreating blade tip is not 
materially improved by hinge offsets up to 20 per cent of the 
blade radius. 

The authors have recently completed a study! concerning the 
above effect. It should be pointed out that Chawla’s result, al 
though correct for the conditions that he imposed, is misleading, 
since an alternate treatment of the theory yields entirely different 
results. 

Whenever a first harmonic aerodynamic loading is produced on 
a rotor element, a hub moment proportional to this harmonic 
loading times its distance to the center of rotation must result 
Therefore, if the hub moment is forced to zero either by employ 
ing a zero offset hinge or by cyclic feathering of an offset or rigid 
rotor, then the first harmonic aerodynamic loading must also go 
to zero or disappear. The result is to sacrifice the possibility of 
obtaining a variation in angle of attack which would postpone or 
eliminate stall in the reverse flow region. Chawla has chosen 
6, so as to produce a zero hub moment and, therefore, eliminates 
at the outset any chance of revealing the overall effect of flapping 
hinge offset on the stall of the retreating blade tip. In addition, 
he has not included a sufficient variation in blade mass constant 
If the cyclic pitch condition is not imposed and a hub moment is 
accepted, it will be found that the flapping hinge offset is aero- 
dynamically significant for proper combinations of offset and 
blade mass constant. 

The local blade angle-of-attack variation with azimuth is na 
turally determined by the loading distribution and the flapping 
motion required for equilibrium about the flapping hinge. The 
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Fic. 1. Effect of blade mass constant and flapping hinge offset 
on the stall of retreating blade tip, no cyclic pitch, blade loading 
(C/o) =.10 


experimental results given in reference 2 indicate that a pro 
nounced decrease in the maximum angle of attack is experienced 
by the blade tip at all advance ratios for the rotor tested. The 
theoretical approach given in reference 1 verifies these experi- 
mental results to a close approximation. 

As an example of the effect of offset hinge taken from reference 
1, the combination of an offset of 20 per cent and a blade mass 
constant of 4 results in the following characteristics at an advance 
ratio of 0.60: 

(1) A concentration of aerodynamic loading at ¥ 140 

(2) First harmonic flapping coefficients of a, = 5.7° and ) = 
-—5.0°. 

(3) A maximum blade tip angle of attack of 14° 

(4) Tip path plane approximately horizontal with shaft axis 
inclined forward 5.5°. 

(5) No sacrifice in lift. 

A theoretical comparison of maximum angles of attack over a 
large range of advanced ratios between rotors of zero and 20 per 
cent offset is given in Fig. 1 as further indication of offset effect. 

The question of possible full-scale rotor application of the offset 
rotor of course can only be answered after the investigation and 
evaluation of such unknowns as propulsive force characteristics 


weight penalties, vibration, controllability, and stability 
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mental Aerodynamic Loading on a 
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Approximate Solutions for the Detached Shock 
Problem* 


S.F. Borg 
Stevens Institute of Technology, Hoboken, N.J 
September 4, 1952 


& DERLEY! POINTED OUT the essential simplification obtained 
by considering the hodograph plane for the transonic de 
tached shock flow about a wedge Vincenti and Wagoner” 
solved the problem using relaxation methods. 
describes a rapid procedure for obtaining approximate values of 
These values may then be used either as 


The present note 


the stream function. 


* This work was done while the writer was a Research Engineer at Grum 


man Aircraft Engineering Corporation, Bethpage, L.I., N.Y 
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the first values of a relaxation solution or, depending upon the 
accuracy desired, directly for transforming back to the physical 
plane 

Using the notation of reference 3 except that unbarred letters 
/ 


are used throughout, the problem requires finding a function y¥ 


which satisfies the equation 
(O%y/dn?) — 2n(d2y/d02) = 0 (1) 


in the region OBEAO, Fig. 1, subject to the boundary conditions 
shown. 

For constant values of », Eq. (1) becomes, under suitably 
transformed coordinates, Laplace’s equation. Hence, in narrow 
strips adjacent to » = constant, the equation of the field will be 
approximately satisfied by potential functions. In particular, 
the boundary conditions may also be satisfied or approximated by 
means of a potential doublet pair reflected about axis OB. We 


must now distinguish between two regions. 


(I1)—Recion CBEA 


In this region, by using doublets, we may satisfy (a) the differ 
ential equation in the strip 7 = constant; (b) the boundary con 
dition along AE; and (c) the boundary condition along CB. 

It is possible, based upon physical reasoning, to determine be 
forehand approximate doublet locations (1, 60), (m0, 20. — @) 
However, in view of the similarity of the stream pattern for various 
values of 6, and because of the form of the dependence of the 
stream pattern upon 7, it would be preferable to proceed as fol 
lows: Determine 7 for the case @, = 1.6 (which is solved in 
reference 3) so that good agreement is obtained for the stream-line 
numbers. For other values of 0, use the same values of no. % 
will vary corresponding to the change in @, 


(I1)—ReGion OCAO 


In this region, the differential equation and the boundary con 
dition along OC are satisfied However, the dual requirement 


along OA 


y = 0) 
oy /dn = Of 
cannot be satisfied by the simple doublet pair. It is possible, 
however, to approximate the required boundary conditions in 
any one of several different ways. For example: 

(A) We can satisfy 0¥/O0n = 0 and can locate the doublets so 
that (a) ~Y = O occurs for small values of 6; (b) on physical 
grounds Yar. will occur at @,/2 for values of —y somewhat 
greater than —1. We may locate the doublets in such a way as 


to obtain y,,er. at the estimated value of @ 
7 
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The products F(y)a(y) and R(y)a(y) cannot be used for comput- 
ing load or moment distributions for the forward wing. 

In this note, approximate solutions are developed for the span 
loadings, F(y) and R(y), at supersonic speeds for a flat reverse 
delta wing with subsonic trailing edges. These solutions, to- 
gether wit h the reciprocity theorems, provide a convenient method 
for computing the total lift or rolling moment on a narrow for 
ward delta wing having arbitrary spanwise geometric or apparent 
The latter results from wing-wing or wing-body inter- 


twist 
For such problems, the downwash can usually be con 


ference 
sidered to vary only along the span 

The following solutions are based on nonviscous linearized 
supersonic thin-wing theory. The wings are assumed to be un 
yawed. The rolling moment is computed with respect to the 
wing centerline. 


A) Span LOADING FOR REVERSE DELTA WING AT CONSTANT 
ANGLE OF ATTACK, @ 


The solution to this problem is known for two cases. The 
first occurs when the tip Mach cones are coincident with the trail 


ing edge of the wing. In this case the span loading is triangular 


and is given by 
l AL ta ly 
er (3) 
qs A(y/s) 8 5 


(Sis the wing area, g is the dynamic pressure, s is the semispan, 
y is the spanwise coordinate measured from the wing ceiterline, 
8 is the Mach Number parameter Y J/? — 1, and d is the edge 
parameter 8 tan 6, where 6 is the semivertex angle of the trailing 
edge. In Eq. (3), d is equal to unity and need not be written ex 
plicitly; it is included for later convenience.) The second case is 
the slender wing theory result (d = 0). Here, the spanwise load 
ing is elliptic and independent of Mach Number.* 


] AL ta ] , (? 2 F 
= ¢ (4) 
gS A(y/s) aN 5 


For the intermediate cases, 0 < d < 1, partial load distributions 
may be found by applying the reciprocity concepts in the reverse 
of their customary usage. That is, the loading at a station near 
the tip of the reverse wing is the same as the exact “‘strip’’ solu 
tion for the lift due to the corresponding deflected chordwise 
element of the forward wing. This solution was derived in refer 


ence 4 for a flat delta wing with subsonic leading edges as 


] AL la ( | y » 
= i : (5) 
gS My/s) BY 5 ' 


Eq. (5) is valid at any spanwise station outboard of that at which 
the tip Mach cone first intersects the opposite edge—i.e., |y/s| > 
(l1—d)/(1 +d). At the limit |y/s|) = (1 — d)/(1 + d), Eqs 
(4) and (5) give identical results. In addition, Morikawa® has 
shown that the lift due to the central strip of a delta wing is inde 


pendent of Mach Number. Therefore, at the wing centerline, 


l AL la 
; = d (6 
gS A(y/s)] 4 a 8 


(B) SpAN LOADING FOR REVERSE DELTA WING ROLLING WITH 
CONSTANT ANGULAR VELOCITY, p 


Sunilar solutions can be found for the rolling reverse delta wing 


The slender wing theory result (d = 0) is® 


I AL 1 p ae (? 2 _ 
= U (4) 
(2s )Sg A(y/s) BV § \ s 


When the tip Mach cones lie along the trailing edge of the wing 


(d = 1), the loading is’ 


l AL 2p .9 ( y 
= = d ] (S) 
(2s) Sq Al y/s B I ) .) 


For intermediate values (0 < d < 1), the correct loading on the 


-10 -8 -6 -4 «= 0 =z J 6 8 Lo 
(4) 
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Fic. 1. Lift influence function for delta wings. 
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Fic. 2. Roll influence function for delta wings 


outer part of the wing is the same as the exact solution for the 
rolling moment due to the corresponding deflected (at angle p/V) 
chordwise strip of the forward delta wing. This result is found 


by the method of Evvard* to be 


l AL ip 3 ( 2) 
(2s)Sq ACy/s) a 2p V Vd , ‘ 


y 
E + v(?)s (1 | (9) 


when y/s| > (1 —d)/(1 +d). [The sign of the last term of Eq 
(9) must be opposite to the sign of y/s.] At |y/s| = (1 d) 
(1 + d), the loadings computed by Eq. (9) and Eq. (7) are equal 
Furthermore, it is known that the lift at the centerline of the re 
verse wing (or the rolling moment due to deflecting the central 


strip of the forward wing) is zero, 


l AL 
= () (10 
. (2s) Sq A(y s) y/s 0 


The span-loading parameters, 


r(’) 1 Bp AL 
Ss 7 ad qS Al y/s 


from Eqs. (3)-(6) for the constant angle of attack wing, and 


y l B AL 
R = . 
s (p/V)d (2s)Sq A(y/s) 


from Eqs. (7)-(10) for the rolling wing, are plotted in Figs. 1 and 2, 
respectively. These curves suggested that a reasonable approxi 
mation for the loading of the inner sections of the reverse delta 
wing would be obtained from slender wing theory. That is, in the 
region 0 < |y/s|; <(1 d)/(1 + d), the loading is assumed to be 
given by Eq. (4) for the wing at constant angle of attack and by 
Kq. (7) for the wing that is rolling with constant rate of roll. In 
the region.(1 d)/(1 + d) < y/s| < 1, the loading is given 


correctly by Eq. (5) or (9) 
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Since the reverse wing span loadings, if correct, are also influ- 
ence functions for computing total lift or moment on forward 
delta wings, the reliability of the assumed loadings was checked by 
using the analytical formulation of the curves in Figs. 1 and 2 to 
compute the lift and rolling moment coefficients for some for- 
ward delta wings for which exact solutions are known. In the 
lift case, comparisons were made for delta wings having constant 
angle of attack? and linear symmetric twist.? The assumed lift 
influence function led to total lifts that were lower than the true 
ones by less than 2 per cent. In the rolling moment case, com 
parisons were made for a delta aileron and for the damping in 
roll of a plane delta wing. These calculations employed the roll 
influence function from Fig. 2, where, for any d, the influence 
function inboard of the station |y/s_/ = (1 — d)/(1 + d)is approxi 
mated by the slender wing result. The solution obtained by this 
method for the delta aileron problem was the exact result reported 
in reference 10, BC; = a(2/3)d. This means that there is the same 
total area under any one of the curves in Fig. 2, and the amount of 
area is correct. For the second example, damping in roll of a 
plane delta wing, the present method yields an analytical expres- 
sion different from the exact result ;!! yet the numerical difference 
between the two solutions is less than can be determined with 
four-decimal-place accuracy. This agreement is probably due to 
the fact that the local angle of attack of this wing is relatively 
small over the inner part of the span where the solution is ap- 
proximate and is larger over the outer span where the influence 
function is exact. 

From these examples, it is concluded that, for any Mach Num 


ber, slender wing theory gives a gcod approximation for the span 


Energy Removal from a Supersonic Turbulent 
Boundary Layer 


John Joseph Martin 
Research Engineer, Aerophysics 
Aviation, Inc., Downey, Calif 


October 6, 1952 


Laboratory, North American 


_— SUPERSONIC BOUNDARY LAYER Is apparently a convenient 
source of thermal energy for use within aircraft propelled at 
high velocities. Before considerable time is spent determining 
a suitable means for utilizing this energy, an analysis to deter- 
mine the desirability of such a utilization should be carried out 
Let the system being analyzed be represented by a section of a 
flat plate of 1 sq.ft., the restrictions being that the section’s dis 
tance from the leading edge be sufficient so that the dimension 
of the 1 sq.ft. in the direction of flow is small and that one side 
These condi 


” 


of the plate is “‘inside’’ and the other ‘‘outside 
tions are imposed so that the skin-friction and heat-transfer co 
efficients (cy, and A) will be essentially constant over the 1 sq-ft 
ind so that the high temperature on one side of the plate may be 
used on the other. Assume, also, that the following constants 


obtain for air 


Molecular weight: m = 29.0 Ibs.m/mole 

Gas constant: R = 1,544 ft.lbs.h7/mole °R 
Isobaric specific heat: c, = 0.24 B.t.u./Ib.m °R 
Specific heat ratio: y = 1.40 


Mechanical equivalent of heat: J = 778.2 ft.lbs.¢/B.t.u 
Acceleration due to gravity: g = 32.2 ft./sec.* 


In addition, the following relations are used: 


Speed of sound: a,? = ygRT,,/m 
Drag force: Dy = C7 .02Uq*/2g 
Drag power: P = u,,D;/J 


Heat transferred: g = A( Tins 1 
Temperature ratio: 6; = 7,/T 
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loading on a reverse delta wing inboard of the spanwise station 
at which the tip Mach cones first intersect the trailing edg, 
Outboard of this station, exact loadings are known. These soly 
tions are useful in reducing the complexity of many problems 
involving the calculation of lift or rolling moment on forward 


delta wings in nonuniform flow fields 
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Insulated wall temperature ratio: @Oins. = + (1/2 


(4 1).W.? 


The primary requirement in this analysis is a set of values fot 
cy, and h which is consistent with variation of @,, Such a con 
sistent system is available in reference 1 and is, in fact, the basis 


of this note. If the assumption is made that the only net work 


that must be done in propelling the plate is the difference (? 
ng), then there is a basis for comparison of various pairs of power 
requirement and heat flux. Making this assumption is the same 
as saying that a percentage (», which may correspond to a heat 
engine overall efficiency ) of the heat extracted from the boundary 
layer is useful in some way in driving the plate and that this 
energy 1s available at the temperature 6 
The power required is 
P=Uu,Dy/J = (Cf Pollo! w/2) (y¥RM ,,2/Jin 


The heat flux is 
q = h( Tins Tw) = (C6oPaotol o/2) (Gms 6 
since, according to refe-eice 2, cy = h/p,uoG rae 
Hence, if 7 is a heat e-gitie efficiency, the net work is 


(P nq) = Poe = (CfPotlal~/2)((yRM,,?/Jm) 


n¢,,( Bins 6 
or, Sizce Ojins. = 1 + [1/2(y L) MV 
Pow = (C¢aP otal o/2)(1/Ay 1)M.°A2 n) + (Oy 1) ¢ 
in which c,, = c,.(@u, R 
For \/ 5.0, R.. = 107, and various values of 6, and 7, Fig. | 
shows the variation of $6,/®éjns.. For values 0 n 0.4 (ap- 


proximately), the net work apparently increases with decreasing 
§,, so that heat extraction from the boundary layer is undesir- 


able. For 0.4 S n 1.0, the net work is decreasing with de- 

* A more accurate value of c~ is given by Van Driest as cs ./2P 
In Fig. 1, Pr = 1 was chosen for convenience; if Pr # 1, a modified efficiency 
is given by n’ Pr*/in 





POWER REQUIRED RATIO, %p,, | Pains 
° 


Fic. | 
wall | 


creasil 
the fa 
clencic 
from <¢ 


ess 


Vai 
Journa 
2 Val 
ne 
Americ 
er3, | 


Rela 
fc 


Profes 
Bou 
Octo! 


It is 
funda 
Ma is 
funda 
$s vali 


discar 


Ex 


wher 
the r 


body 


tion 
rdge 
olu 





lems 


Vari 


iry 


ne 








READERS’ 














14 — ae T 
r- 4 | 
| 
| 
12 

o 

c 

o 

9 | 
ss +6 
z | 
> 

ry | 
e | 
| 
= 08 

°o 

ve) 

3 

5 | 
o 

ra 

© 06 | 
& | 
= | } 
£ | 
2 04} = (POWER REQUIRED) - (NET WORK FROM HEAT TRANSFER ) 


| 
| 


| 
| | 
02 { | + + t 








| 


| 
1 
2 3 





| 
| 
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Fic. 1. Variation of dimensionless net power required versus 
wall temperature ratio with parametric heat engine efficiency 


creasing 0, so that heat extraction seems desirable. However, 


the fact that most stationary heat engines have overall effi 
extraction 


ciencies less than 0.4 means that, generally, heat 


from a supersonic (urbulent boundary layer is an inefficient proc 


ess 
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Relation of Force, Acceleration, and Momentum 
for a Body Whose Mass Is Varying 


H. W. Sibert 

Professor of Aeronautical Engineering, University of Colorado 
Boulder, Colo 

October 2, 1959 


SUMMARY 


It is the belief of many engineers that the relation force 1(MV)/dt is 
fundamental for both constant and variable mass, while the relation force 
Wa is true only for constant mass It is shown that the latter equation is 
fundamental for both constant and variable mass while the former equation 
s valid for variable mass only when the added mass is initially at rest or the 


discarded mass is finally at rest 


EQUATIONS 


Joxrexnn NTS HAVE DEMONSTRATED that 


Fres. = Ma V (1 


(constant 


where .V/ is the mass of a body, a is its acceleration, and Fres. is 


the resultant force acting on the body. If V is the speed of the 


body, 


d( MV )/dt (constant 17) (2) 


Fre W(dV/dt) 
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Since a jet airplane or missile has% high rate of fuel consumption, 
its mass cannot be considered to be constant Hence, the de- 
signer tnust know whether formulas that have been derived for 
constant mass can also be applied to a body of varying mass. In 
particular, can either or both of Eqs. (1) and (2) be applied to 
a body whose mass is varying rapidly? 

To avoid the use of vectors, all forces and velocities will be as- 
sumed to be parallel. Let the body be picking up mass, let all 


velocities be measured relative to the ground, and let 


V = mass of the body at any time / 
V = velocity of the body at time ¢ 
= initial velocity at time ¢ of the mass that is being 


picked up 
t, exclusive 


F..,. = external force acting on the body at time 
of any force due to jet action 
During an infinitesimal length of time from ¢ to ¢ + A¢, the mass 
I’ ») S 


picked up is (dM /dt) At, and its increase in momentum is ( | 7 


(dM /dt)At During time Af, the increase in momentum of mass 
Wis WAV. Since impulse equals change in momentum, 
Fext.At = MAV + (V dM /dt) At 3 
Dividing by At and then letting At ~ 0, 
Fext. = M(dV/dt) + (V vid WM /dt } 


Eq. (4) is given in reference 1. 
Note that ./(dV/dt) is the rate of change of momentum of the 
instantaneous mass while ( V v)d.\J/dt is the rate of change of 


momentum of the added mass. Hence, Eq. (4) becomes 


Fext. = rate of change of momentum of instantaneous mass 4+ 


rate of change of momentum of added mass (5 
At first it might be thought that Eqs. (2) and (5) are identical 


rhis is not true because 


d( MV WidV/dt + Vidi dt 


dt = 
and, hence, Eq. (4) becomes 


Fext. = d( MV)/dt o(d Md 6 


Thus Eq. (2) is valid only when the mass picked up its initially at 
rest and Fres. = Fext.. 

The mass being picked up by the body at time 
body a jet force Fiee = (V — v)d.VJ/dt that acts in a direction 


f exerts on the 


opposite to V Hence, from Eq (4 


Fext. = M(dV/dt) + F ‘ 


Eq. (7) can be written 


W(dV/dt 8 


Fres. = 


with | an = Fox: Fier Hence, Eq 1) is also valid for a body 
of variable mass, provided Fres. = Fex: F 

By a similar analysis it can be shown that Eqs. (4) through (8 
also hold when the body is discarding mass. Moreover, these 
equations will also be valid when the forces and velocities are 
not parallel, provided each velocity and force in them is written 


as a vector 


CONCLUDING REMARKS 


2) is fundamental 


It is the belief of many engineers that Eq. (2 
for both constant and variable mass, while Eq. (1 
They have justified the use of Eq. (1) for vari 

r 
[ 


is true only for 


constant mass. 
able mass on the basis that for most aircraft the rate of change « 


‘ass i j 
mass is small compared with the mass at any instant and, henes 


the effect of the change in mass is negligible It has been shown 


above that Eq. (1) is fundamental for both constant and variable 


mass provided the resultant force includes any jet force on the 
body while Eq. (2) is valid for variable mass only when the added 
at rest or the discarded mass is finally at rest 


mass is initially 
1 jet airplanes or rockets, because 
] 


rhis last condition is not met it 
the jet speed relative to the aircraft is always larger than the speed 
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of the aircraft and, hence, the get is not at rest relative to the 
ground. It should be noted, however, that the resultant force is 
equal to the rate of change of the momentum for both constant 
and variable mass provided that the resultant force does not in- 
clude any jet force and the rate of change of momentum is calcu- 
lated separately for the instantaneous mass and the added (or dis- 


carded ) mass [see Eq. (5) ]. 


REFERENCE 


1 Timoshenko, S., and Young, D. H., Advanced Dynamics, p. 114; Me 
Graw-Hill Book Company, Inc., New York, 1948 


Further Comments on ‘‘Turbulent Boundary- 
Layer Characteristics at Supersonic Speeds’’ 


R. J. Monaghan* 
Royal Aircraft Establishment, Farnborough, Hants., England 
September 29, 1952 


T° A CONTRIBUTION! TO THE READERS’ Forum in the May, 1952, 
issue of the JOURNAL, Luther presented a simple interpolation 
formula based on Wilson’s equation? for turbulent skin friction in 
compressible flow, under zero heat-transfer conditions. 

An equally simple formula can be obtained if we write Wilson's 


equation in the form 


0.242 
aay = logio( R.-Cr,,) (1) 
Cp/* 
where 
@ = sin"! V0/Vo 
with 0 = J/?/(.\/* + 5) if y = 1.4; subscript w denotes that 


density and viscosity are to be evaluated at wall (plate surface) 
temperature, and it is assumed that p « 7-76, 
By comparison of Eq. (1) with the corresponding incom- 


pressible flow formula 


0.242 CH = logi( R;Cr,) (2) 
we see that 
Cr; = Cr, a 
when , (3) 


R = Reet } 

Eqs. (3) can then be applied to the more convenient Blasius or 
Prandtl-Schlichting formulas for mean skin friction in incompres- 
sible flow to obtain corresponding formulas for compressible flow 
Thus, if we take 

Cr; = 0.46(logiwR;)~2-* (4) 


and apply Eqs. (3), we obtain 

Cr, = 046 &*%(log wR?) ~2-6 (5) 
which gives results in close agreement with those from Eq. (1) and 
is more simple in application. 

Wilson’s analysis? starts by assuming the shearing stress rela- 
tion to be of incompressible flow form and by allowing for vari- 
ations in the mean density obtains expressions for velocity profile 
and then for skin friction in compressible flow. If, however, we 
start by assuming that the log law velocity distribution 


. . ’ 
u/ur = Cy + C, log yu,z/v,, (6) 


where u;? = 79/p and C),C2 are the constants found valid in in- 
compressible flow, will be valid in compressible flow, with or with 
out heat transfer, provided density and viscosity are evaluated at 
wall temperature, and if we assume a temperature distribution of 


* The writer is indébted to the Chief Scientist, Ministry of Supply, for 


permission to publish this note 
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the form 
T/Tw = 1+ B(u/m) — A%*(u/m,)? (7 
then by substitution in the momentum equation and by the com. 
parison already described, we obtain® 
Cr; = Cr, | 
when 4 g 
Ri = RATi/Tw)) 
where 7) is free-stream static temperature. 
Applying Eqs. (8) to Eq. (4), we obtain 
Cr, = (0.46 [logioR wl T; ie )] 2.6 (Q 


which gives equally as good agreement with experimental results 
as Wilson’s equation [Eq. (1) or (5)] gives under zero heat-trans 
fer conditions. It also gives results within 10 per cent of the 
values given by an extension of Wilson’s analysis to include heat- 
transfer conditions for Mach Numbers up to 4 and for values of 
Tm/T. between 0.5 and 2.0, with R = 107 (Ty is free-stream 
stagnation temperature). 

The reason for this agreement is presumably to be found in the 
fact that the velocity profiles as given by Wilson’s analysis? are 
forced by boundary conditions to agree with the log law profile 
[Eq. (6)] in the inner portion of the layer, which region is of major 
importance when evaluating the momentum integral. 

Note that the ‘“‘compressibility parameter’’ in Eq. (9) is the 
ratio 7;/T,. Thus, if the formula were valid, the skin friction in 
a high-speed flow under zero heat-transfer conditions should be 
the same as that in a low-speed flow with heat transfer from plate 
to stream for the same value of 7;/7,. A partial check of this 
was obtained by an analysis‘ of the skin-friction and heat-transfer 
coefficients obtained from subsonic flow tests in pipes with high 
heat-transfer rates> when a good correlation was obtained on the 
basis of Eqs. (8). 

On the other hand, recent measurements at the R.A.E. of the 
boundary layers on a flat plate at JJ = 2.48 and 2.82 under heat- 
transfer conditions (7, = 373°K., Ty; = 237°K., 276.5°K.) have 
shown no variation of skin-friction coefficient from its zero heat 
transfer value at the same Mach Number. These results (un 
published M.O.S. reports) were obtained from pitot traverses, but 
heat-transfer rates were measured directly and showed a variation 
somewhat as given by Eqs. (8). 

Finally, the analysis of reference 3 also gave 

H/H; = (Tw/Ti) + ((y — 1)/2] Mi? (10) 
for the variation of the ratio (/7) of displacement thickness to 
momentum thickness with Mach Number and heat transfer, and 
this formula has given good agreement with our experimental re- 
sults to date (with and without heat transfer) if 7; is taken to be 
1.3 and if the second term on the right-hand side is multiplied by 


the cube root of Prandtl Number—i.e., 
H = 1.34(Twe/T:) + Pr'/*[(4 1)/2)A1,?} (10a 


(The original analysis* was for Pr = 1, and it can be argued that 
a logical extension to other values would include the factor P 
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On the Boundary-Layer Equations in 
Hypersonic Flow and Their Approximate 


Solutions} 


Lester Lees 

Department of Aeronautical Engineering, Princeton University, 
Princeton, N.J. 

October 2, 1952 


SUMMARY 


Analytical solutions of the Prandtl boundary-layer equations are ob 
tained for the problem of the “‘strong’’ interaction between the leading 
edge shock and the viscous layer over a flat plate at hypersonic velocities 
As the Mach Number increases and the interaction region spreads down 
stream over the plate, the local skin-friction coefficient increases rapidly 
over its conventional value. The local heat-transfer coefficient at first re 
mains practically unaffected but then also begins to increase with Mach 


Number 


List OF SYMBOLS 
Free-stream quantities are denoted by the subscript ©, while the sub 
denotes quantities evaluated at the surface and the subscript 6 


script 
denotes quantities evaluated at the outer edge of the viscous layer 
r*, y* coordinates parallel and normal to surface, respectively, with 
origin at leading edge of plate 
u*, v* components of velocity parallel and normal to surface, respec 
tively 
p* pressure 
p* density 
T* absolute temperature 
u* viscosity 
y* kinematic viscosity, u*/ p* 
ratio of specific heats, Cp/C, 
6 thickness of viscous layer 
1 characteristic length 
VM free-stream Mach Number 
Re,* Reynolds Number based on free-stream quantities 
p p*/Pao* 
p p*/pa*, etc 
\ x¥*/1., y*/6 . 


DISCUSSION 


I’ A RECENT NOTE, Shen! shows that the Prandtl boundary 
layer equations are valid even for the interaction between the 
leading-edge shock and the viscous layer over an insulated 
flat plate or wedge at hypersonic velocities, provided that 
(6/x*)? <<< 1 However, the argument leading to the statement 
that the static pressure variation across the layer is of order 
1/M.2, and is therefore negligible, seems inconclusive. This 
argument is based on the postulate that p ~ 1/T ~ 1/2, in 
spite of the fact that 7 is certainly not of order ./,,? near the 
outer edge of the viscous layer or near the surface for the case of 
heat transfer. Thus a few important statements must be added 
to Shen’s discussion 

From the y-momentum and continuity equations, py ~ y* 
near the surface, and the contribution of this region to the pres 
sure variation across the layer is negligible, as expected Near 
the outer edge, 

py = yM .*( 6: puvrx + pvv,) + viscous terms 

-yM ,*(6/x*)? 


where the subscript denotes differentiation with respect to the 
indicated variable. But the energy equation shows that 07°/Oy ~ 
M2, so that the thickness of the outer zone in which p = 0(1) is 
of order 1/M 2. Therefore the contribution of this zone to the 
pressure variation across the layer is of order (6/x*)? and is neg 
ligible when (6/x*)? << 1. In the interior of the layer p ~ 
1/\M.2; thus Shen’s argument is completed. 

Once it is established that the pressure variation across the vis- 
cous layer is negligible when (6/x*)? << 1, the problem of the 
interaction of shock and viscous layer over a flat plate for large 
values of M.3/+~/ Rest = x treated by Shen? and by Li and 


Tt This work was supported by the U.S.A.F., Air Research and Develop 
3(038)-250 
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Nagamatsu’ can be considerably simplified. By solving the von 
Karman momentum integral equation approximately for an in 
sulated surface, these investigators find that 6 ~ x/‘ and p ~ 
x~'/? for x >> 1. But this result can be obtained directly by 
means of an order of magnitude consideration. Flows with large 
x are characterized by the fact that the pressure gradient 0p*/Ox 
and the viscous stress term (0/0y*) [u*(Ou*/Oy*)] are of the same 
At the outer edge of the viscous layer the oblique shock 


-* 


order. 
relation 
op ds d% 


= ¥(7¥ + 1)M,? 
Ox = ax* dx** 


holds approximately; therefore, 


> 6° x * u 
a > va : 
© al o- 
or 
. VM 
6 = dy (s* (1) 
se. * . 
and 
Pb = pox (2 
where x = VC x and C is the constant in the linear viscosity- 
temperature relation ~ = CT; 6) is a constant, and po = 


(°/s2)y¥(y + 1)60?.t 

In reference 3 it is also shown that the Pohihausen parameter 
\ determining the shape of the velocity profile u/us = (y, A) is 
approximately constant for xy >> 1. This result, together with 
the relations (1) and (2), suggests that “‘exact’’ solutions of the 
boundary-layer equations can be obtained directly for both in 
sulated and noninsulated surfaces, without resorting to the ap 
proximate Pohlhausen method. The asymptotic series solution 
for u(x, y), for example, is of the form 


u l 
(X,¥) = Uo) + ule 
us x 
where o = y*/(x*)’t and “(o), u:(@), etc.; each satisfies an ordi 


nary second-order differential. equation. In the ‘‘von Mises 
coordinates” x,y along and normal to the stream lines, the x 


momentum equation is 


1 p 
uu, = U(puny)y 3 1 (3 
7 V/ 2 p 
and the energy equation is 

: Y Ss 

T; = T+ (4 1) *pu(uy)? + pluTy)y (4) 
7 p 

Here yy = y*/V v,,*u,*LC, while pu* = dy*/oy* and pv* = 
Oy*/ox*. By employing Eq. (2) for the pressure and intro 


ducing the nondimensional stream function ¢ ~ (Y yp) (1/V x), 


or 


the momentum equation |[Eq. (3)] is reduced to the following 
equation for uo: 


t In references 2 and 3, the slope of the outer edge of the viscous layer 


dé/dx* is regarded as identical with the oblique shock angle @,, so that 


27 dé 
p Vv 
y+ 1 dx* 


But this statement is inconsistent with the fact that 7 
+ 1)/2]e when VM >> 1, according to 


16/dx* by the con 
tinuity equation, whereas 0; [(4 
the oblique shock relations In the present note the approximate boundary 


condition 


9 > 


is employed instead 
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duo d duy or x . M.@/2) 1 
tuo + 2uo io + T) = 0 (5) Cfy ™ 7k 
dt d¢ d¢ yM? (Rez+)/* V Rez 
where the term xuu, is dropped because it is of order (6/x*)? . Min! he)" 


compared to the other terms (us* ~ uo* for Mao >> 1). The 
energy equation is 


OT> y-1. duy \? 
2x + To = (¥ — 1) Maw*M 
Ox a dé 


d ( dT» af dT - 
0 ») 
dé i dé 3 & 


Similar equations are obtained for m, 7), etc., but they will not 
be given here. 
(a) Insulated Surface 


For a Prandtl Number of unity and zero surface heat transfer, 
Crocco found that 


T=1 + [(y7 — 1)/2] Mo°A1 — 4?) 
Therefore, when 1/.. >> 1, Eq. (5) becomes 
du d duo y¥—-—1 _ 
Uo + 2u Uy + ° (1 — us?) = O (7) 
dé dé dé ¥ 


and this equation is exactly the Falkner-Skan differential equa- 
tion’ for the low-speed laminar boundary layer with ug ~ x”. 
If the coordinate n; defined by the relation d{/dn; = up is intro- 
duced, then Eq. (7) takes the more familiar form 


d% ae 


Sorel) 
~ dnt : dn; : dni 


B = 2m/(m + 1) = (vy — 1)/¥ 


where 


In other words, u(¢) is the same known function of ¢ as in the 
Falkner-Skan case with B = (y — 1)/y. (For air, (y — 1)/y = 
0.286.) The corresponding actual distance from the surface is 
found by means of the simple quadrature 


F ™ dé dé dy ; 
y= dn, ; \S 
0 dn / \dy dy* 
or 
y* ! as ”) (y 1) 1 n , 
(Rew) /* = — CAM. : (1 — uo?)dn; 
x* 3 V y(v¥ + :1) do Jo 
(Sa) 
where 
1 
2 V4 l ta” Uco* 


and the quantity in brackets is the sum of the displacement and 
momentum thickness for the Falkner-Skan profile with 8 = 
(y — 1)/y. (The necessary values are given by Pretsch.°) 

The expression for the local skin-friction coefficient for large 


Jin” (“*) (2 
meal po*tto* \ dni /w\ dy" Js 
3 du 25) ae I 
ep = |— Voy + 1b cys (9) 
} dn; (Rezx) /4] V Rez 


Ifu ~ 7”, then 


x is given by 


or 


Thus 


for an insulated surface. 


These results agree with those obtained in reference 3, but for 
various reasons that will not be discussed here the numerical 
constants are different. (See footnote, for example.) These re 
sults can also be obtained as a special application of the Stewart. 


son transformation.® 


(b) Heat Transfer 


In free flight the surface temperature is much lower than the 
equilibrium temperature, and 7) is approximately equal to 
((y — 1)/2]M,,2 u(1 — «u) in the interior of the viscous layer fora 
Prandtl! Number of unity. The maximum value of the term 
(2/yMo?)To in Eq. (5) is (y — 1)/4y, which is small. (For 
air, (y — 1)/4y = 0.071.) Therefore, uo(¢) approximately satis 


fies the differential equation 


¢ duo d ( =) 
+ Uy = () 
2 at dé dt 


and is therefore approximately equivalent to the Blasius velocity 


distribution in the ¢-coordinate. More precisely, 


uf) = up(f) + Au(¢) 


~% 
fat) «= ¥. +(1 + 1— Ma? - Te) w 


- J 
? M w2uo? + f(o) + T5(x)eli 


where Au(¢), f(¢), and g(¢) represent small corrections satisfying 
ordinary linear, second-order inhomogeneous differential equa 
tions; these equations are now in process of being integrated 
With good accuracy the actual distance normal to the surface is 
J ¥ 

given by 

¥" | : eh 

; (Re,« ‘t= 4 ‘V/ 
¥ 3 Vy(y¥y4+ 1 


I , 
Nol 1 u an 
6 0 


where 

. ~ Vy 1 V 0.664 
V3 Ly(y + 1) 

and w(7;) is the Blasius function. The approximate expression 

for local skin-friction coefficient is 

0.664 ¢ 


i % 0.664 ¢ X wh 10 
Ve Y 


. Vs. “a7 
ey Ny tI V4 
and the local heat transfer coefficient is given by 
ho = qo/Po Ua *¢,(T,* To") = ¢y/2 
where 7.*/7.* = 1 + [(7 1)/21M ~? 


C m= {T,*/T 


x = VC(M V Re,x ) 


These expressions (9) and (10) for local skin-friction and heat- 
transfer coefficients are expected to be good approximations for 
x > 15 (approximately); presumably these solutions can be im- 
proved by calculating the 1/x term in the asymptotic series If 
the upper limit of applicability is taken roughly as (6/x*)? = 0.10, 
then the minimum value of local Reynolds Number at which these 

9x 


= 40 


solutions are applicable would be about 1.6 X 10‘ for Wa = = 
for example, for the case y = 1.40, Ty*/Tao* = 5, and m = 0 76 
However, the slip-flow corrections are no longer negligible when 


Ma/V Re,« > 0.10 (see reference 7), or for Rese < 10% at Wo = 
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FIGURE-!. LOCAL LAMINAR SKIN-FRICTION 
AND HEAT TRANSFER COEFFICIENTS FOR 


FLAT PLATE IN HYPERSONIC FLOW. 
RATIO OF SURFACE TO FREE-STREAM 
TEMPERATURE =5.0 ¥=1.4, “~T™ m=076 


10, for Re,+ $ xX 10‘ at Wo = 20, and for Rez« < 9 X 10% at 
Meo = 30. 

All of the above results can be extended to the case of arbitrary 
Prandtl Number, to the case of the flow over a wedge, and prob- 
ably also to the flow over slender bodies of revolution. 

Solutions of the boundary-layer equations for small values of 
\(X < 6, approximately) have already been obtained as asymp- 
totic series in ascending powers of X.5 Except for a region in the 
immediate vicinity of the leading edge where slip-flow exists, a 
large portion of the range of practical interest is now covered. 
In Fig. 1 the local skin-friction and heat-transfer coefficients are 
plotted as functions of Re,« for séveral Mach Numbers, for the 
casey = 1.40, T,*/T.* = The first approxi- 
mation to the skin-friction and heat-transfer coefficients for X >> 
lis given by Eq. (10), while the second approximation includes a 
temin 1/X. As the Mach Number increases and the interaction 
region spreads downstream over the plate, cy increases rapidly 

The local heat-transfer coefficient 
unaffected’ but then also begins to 


5, and m = 0.76. 


over its conventional value. 
at first remains practically 
increase with VV. for large x. 
cussion will be presented in a forthcoming report 


More extensive results and dis 
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On the Use of the Neutral Point as a Stability 
Parameter 


Judson R. Baron 

Naval Supersonic Laboratory, Massachusetts Institute of Technology, 
Cambridge, Mass 

September 29, 19592 


oe LONGITUDINAL STATIC STABILITY of an aircraft is given 
by the rate of change of the pitching moment coefficient with 
respect to the lift coefficient. Considering conventional subsonic 
aircraft, the well-known result follows: The above derivative is 
equal to the (nondimensional ) distance (the ‘‘static margin”’) be- 
tween the center of gravity and neutral point locations, the latter 
being aft of the former for a stable configuration. The neutral 
point is defined as usual as that point at which both the moment 
and the rate of change of moment with angle of attack vanish. 
The simplicity of the result has made it extremely useful. How 
ever, it should be noted that the above is based upon an assump 
tion of constant elevator effectiveness (rate of change of moment, 
or force, with control deflection). When considering super- 
sonic missiles, the presence of nonlinear downwash effects pre- 
clude this assumption. The consequences are outlined below. 
Assume only that the elevator effectiveness is independent of 
the control deflection for a fixed angle of attack. The moment 


about a point 6 in terms of information applicable at point a is 


‘ — oC u OCyN : 
(Ca CyX* + ( ) ( )x*]s (1) 
06 Ja 06 


ire the moment and normal force coefficients, 6 is 


(Fig. 1) 
(Cub = 


where Cy, Cy : 
the control deflection, and X * is the nondimensional transfer dis 
tance from a to }, positive forward. Differentiating Eq. (1) with 


respect to angle of attack a with 6 constant yields 
(<<) - (<<) (1°) xs ' 
da } da a da 
d (> ” - (2% ole én 
da 06 /a 06 


point 6 as the neutral point, it is re- 
The 


Momentarily considering 
quired that the left-hand sides of Eqs. (1) and (2) vanish. 
resulting relations may be used to rewrite Eq. (2) in the form 


. IC 
(< *) = & Xy*)— 
da b da 
x3 _— d CS ( a.) d (26s (3 
* 4 * 6 - 0) 
‘- ‘ ait’ * 06 , da 05 Ja 


where 6 is again a general point. The subscript N refers to the 
Of interest is dCy/dCy evaluated at the center 


Letting ’a and 6 coincide at the center of gravity, 


neutral point 
of gravity. 
Eq. (3) gives the in-trim stability as 


dC j . . d (2s 
= Hy 2* ~ % * 1 y : = 
(“~) - (Xe . )| + On toy \ 8 


(3 5 d ‘Ge \ 
= r) - ~ 
. “ dCyn 06 : aN 


(4) 
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where the subscript 7 refers to the trim condition. In addition, 


since nonlinearities are due to downwash effects, it may be shown 


that 
d (26s d Ge }* (5 
= — »”) 
d( y 06 dCy 06 


where /* is the nondimensional tail length from the moment 
center to the aft surface center of pressure. Thus, from Eqs. (4) 


and (5), 


(<<) ies JX = 
dCn Jo.g..1 
ee ae oe 
( . )s rT (dr — dy) ( ; ! “ 
- dCy 05 Jen | 


from which the classic result is obtained when the elevater effec- 
tiveness is constant. It is apparent that in the case considered 
here the stability is no longer a simple function of the center of 
gravity and neutral point locations. 

The criterion for a stable attitude is that the quantity within 
the braces of Eq. (6) be positive. When 67 = 0, the bracket 
vanishes in Eq. (6), since then a = é6y = 0 for a symmetric missile. 
The familiar form results and the c.g. must be forward of the 
neutral point for a stable configuration. Remembering that 
(Yo.g.* — Xy*) 2 0 when (67 — 6y) 2 Oif a positive deflection in- 
crement produces a positive moment, it follows that for 67 + 0 
a stable attitude results when 


(d/dCy) (0Cy 05)e.g. S A (7) 
where 
Rau’ = Ae 
A= : us : (8) 
|. Fas Vy*)/l* lin + (6r — by) 


The “less than’’ sign applies to the c.g. forward of the neutral 


point, and vice versa. Writing 
(d/dCy) (0Cy OS)eg. = kA 
where k = k( Cy), simplifies Eq. (6) to 


(dCu/dCn e.g. t = —(Xe-g.* Xy*) (1 — k) (9) 


from which it ts more easily seen that & S 1, in addition to the rela 
tive positions of the c.g. and neutral points, determines a stable 


attitude. 
A graphical interpretation follows by noting that in a practical 


case 

((Xeg.* — Xn*)/l* bn = 0 (10) 
From Eqs. (7) and (8) then 
d (2Ex) ohs Xe.n.* — Xn* _ 
dCy O86 sex 67 by 

1(Cur)v/(Cw of (OC 1/05). 

(Car) v/(OC 4 /05)e.. | + Cy, ti 

at a specific Cy = Cy,. Therefore, on a plot of moment effective- 


ness as a function of normal force (Fig. 2), one may determine the 
required fore and aft position of the c.g. with respect to the neu- 
tral point by comparing the slope of the curve at Cy, with the 


slope of the line passing through both (0C/065)..¢., cy, and the 


origin. 


AT Cy=Cy,: 







STABLE, C.G.AFT OF NP > 


oCu 
35 cg 


—. 


STABLE, C.G. FORWARDoF 
NP, k<I 


Cn so Cn, 


ELEVATOR MOMENT 
EFFECTIVENESS, 


NORMAL FORCE COEFFICIENT, Cy 


Fic. 2. 


In the vicinity of Cy = 0, in all cases k < 1 and the c.g. must 
be forward of the neutral point for a stable attitude. If the ef- 
fectiveness variation is such that k 1 for some Cy = Cy,, then 


at an intermediate value, 0 < Cy < Cy,, neutral stability (R = | 
is implied irrespective of the c.g. location [within the validity of 


Eq. (10)]. 


On Asymptotic Solutions for the Heat Transfer 
at Varying Wall Temperatures in a Laminar 
Boundary Layer with Hartree’s Velocity Profiles 


H. Schuh 

Royal Institute of Technology, Division of Aeronautics, Stockholm 
Sweden 

October 23, 1952 


_.. AUTHORS! ® HAVE MADE ALLOWANCES for a variation 
of the wall temperature in thermal boundary layers in lami 
nar flow. The problem was simplified by assuming constant 
fluid properties and neglecting heat produced by friction and com- 
pression. Despite the variation of the wall temperature, the 
partial differential equations for the temperature field can be re 
duced to an ordinary differential equation, provided the tem 
perature difference between wall and free stream is proportional 
toa power of x.! > 

te t; = const. x” (1 
where ¢,, = wall temperature (= variable), f; = temperature out- 
side boundary layer (= constant), m = anexponent. The ordinary 
differential equation so obtained reads 


6 _ dé df = 
+ of — no(2 — B) 6 
dn? dyn dn 


ll 
iS 


where the following quantities have been introduced: a dimen 


sionless wall distance n = (y/VW2 — 8B) Vuj/vx (y = distance 
normal to wall, x = distance parallel to wall, u, = x" = speed 
outside boundary layer, v = kinematic viscosity, and 8 = 2m- 
(m + 1), parameter introduced by Hartree®); the function f de 
fined by ¥ = V2 — B Vvxuf(n) (¥ = stream function ) and cal- 
culated by Hartree;® a dimensionless temperature 6 = (f — 4) 
(t~ — 4); and o = c,p/k = Prandtl Number (c, = specific heat, 
um = dynamic viscosity, and k = thermal conductivity). The 


boundary conditions for Eq. (2) are 


oe ae = 2; n> o,0=(0 (3 
With three param 


eters 8, n, and o, the amount of labor for fairly general solutions 


Eq. (2) can only be solved numerically. 
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It is interesting to consider whether asymptotic 


is considerable. 
glutions of Eq. (2) are not easier to obtain for large er small 


yalues of the parameters or a combination of parameters, es- 
necially since these solutions are often good approximations over 
= 


q wide range of the parameter. 

On inspection of some solutions of Eq. (2), the thickness of the 
thermal boundary layer is found to become small compared to the 
velocity boundary layer if either » or o becomes large. Then it 
is sufficient to replace the velocity profile by its tangent at the 
wall, since, for calculating the temperature field, only that part 
of the velocity profile is of influence which lies within the thermal 
boundary layer. With the exception of the separation point and 
its immediate surroundings, where the velocity gradient at the 
wall vanishes, we put 


df/dn = bn (4) 


where b = constant = df/dn at » = 0. Now Eq. (2) becomes 


d*6 b . dé 
2 nab(2 


: o-% (5) 
dy? 2° dn 


B)n? = 0 


The first asymptotic solution is obtained by introducing a new 


independent variable 


c= n' ’n (6) 
in Eq. (5), 
d*6 ab de " 
c? aobh(2 B)c@ = VO (7) 
de? n 2 d¢ 


For values of ¢/n — 0, the middle term in Eq. (7) can be omitted; 
then we obtain 
ab(2 (8) 


(d°6 d¢?) B)c@ = VO 


this differential equation can be solved in a closed form; a solu 
tion, which fulfills the boundary conditions (3), is 


: 2 
“r( ) eu 1 (r”) Ii 
3 


, 9! ‘ 
rf = (279/3"") [ob(2 


@=2 (r*)| (9) 


where 

B)]'*¢" 

and J, the modified Bessel function (Bessel function with purely 
imaginary argument ) of the first kind and order v as defined for 
instance in reference 7 and I’ the Gamma function.* 

If q is the heat flow per unit surface area, the heat-transfer co- 


efficient h becomes 


9/5 1/3 
q ['(2/3) (abn) b My (10) 


ty t, 3791(4/3) (2 — B)'’* V., 
For convenience we note that 


*/31(4/3) = 0.72895 

For a flat plate and o near 1, the asymptotic solution approxi 
mates the accurate solution (see Table 1) by about 1 to 2 per cent 
> 4, if allowance is made for a systematic 


atm = 4 and less for n 


small error of 2 per cent in reference 3; for the same case, at o = 
20, the error at m = 4 is slightly bigger 


our theory, since the accuracy of the asymptotic solution depends 


This is consistent with 








onthe quantity ”/o. For the stagnation flow, the approximation 


deviates by about 3 to 5 per cent at m = 4. The good agreement 
for higher o-values even at lower values of n if fortuitous since the 
two influences cancel each other: (a) replacing the velocity pro 
file by its tangent at the wall increases 4 for 6 > 0; (b) neglecting 
the middle term in Eq. (7) reduces h 

A second asymptotic solution can be found which is valid for 
We introduce a new variable 


c—- « 


with no restriction on 
(11) 


. 1/7)! 
¢=d0 b/*n 


in Eq. (5) and obtain 
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TABLE 1 
h/k V/vx/iti, 


Flat Plate (8 = 0, 6 = 0.4696) 


o = 0.72 : o = 0.7 
Ref.3 Ref. 5 Asympt. Ref. 5 Ref.9 Asympt 
n 
l 0.489 0.48138 0.4524 0.4770 0.480 0.448 
2 0.597 0. 5863 0.5700 0.5811 0.588 0.565 
3. 60.684 0.6633 0.6525 0.6573 0.646 
4 0.744 0.7255 0.7182 0.7193 0.723 0.711 
5 0.799 0.7737 
10 «1.006 0.9747 
o=z=2 o = 10 . o = 20 
Ref.5 Asympt Ref. 5 Asympt. Ref.5 Asympt 
n 
] 0.6783 0.6359 1.149 1.087 1.431 1.370 
2 0.8238 0.8012 1.389 1.370 1.729 1.726 
3. 0.93806 0.9172 1.563 1.568 1.938 1.971 
4. 1.0161 1.0096 1.701 1.726 2.1038 2.175 
Stagnation Flow (8 = 1.0, 6 = 1.2326) 
o = 0.7- o =5 o = 10 
Ref.5 Asympt Ref. 5 Asympt. Ref.5 Asympt 
n 
l 0.708 0.694 1.432 1.3386 1.818 1 684 
2 0.846 0.874 1.683 2.159 2.121 
3 0.957 1.001 1.914 1.927 2.418 2.428 
4 1.048 1.101 2. O86 2.12) 2 630 2 672 
d’@—s =? db 
} «tO = 0) (12) 
de 2 dt 
with e = n(2 8). This is a differential equation with only one 


parameter instead of three as in Eq. (2). The author found no 
solution of Eq. (12) in a closed form and made no attempt to solve 
it numerically. It should, however, be useful for supplementing 
the known solutions of reference 5 for values of o above 20 and 


any value of n. 
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Graphical Solution for ‘‘Turbulent Boundary 
Layer in Compressible Fluids’’ 


John Joseph Martin . 

Research Engineer, Aerophysics Laboratory, North 
Aviation, Inc., Downey, Calif 

October 24, 1952 


American 


\ J AN DRIEST, IN A RECENT PAPER,' derived equations for local 
and mean coefficients of friction and local coefficients of heat 


transfer. These equations are applicable to a flat plate with the 
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restrictions noted by their author and generally involve trial- 
and-error solutions. It is the purpose of this note to show a 
means by which they may be solved directly by graphical means; 


these solutions were made at Dr. Van Driest’s suggestion. 


The equation derived for the local coefficient of friction C;.,. is 


0.242 wer aor 
———— * (sin x + sin 3) = 
A (eq)? (Tw/T @)'”? ‘ 
;. 1+ 2w Tu me 
0.41 + logw Rk, o. — me login / (66) 
and for the mean coefficient of friction C; , 
0.242 rere 
* (sin + sin 3) = 
ACE; y" i om Fe: i/e a h 
; 2w Tw . 
login Ra Cr logio ~~ (71) 
4 Pa 
Using the definition, A = [('/2) (y — 1 \M 37)'?# Tu ’ ae )' 


one may write, for c;.,, 


0.242 


I('/2) (y — 1M, 
2.57 R 


1 
(G4) ee 
| (T/T )8 


0.242 
— 1)]'“M, 


(sin”' @ + sin 


» + logw cy | (66a ) 


(sin~'a + sin~! 8B) = 


and for C; 


z. ; 
(Gia) ‘tos (T/T) + 20/2 + logio Cy | (7la) 
u x 


Now the left-hand sides of Eqs. (66a) and (71a) are identical 
and may be called by y, which is a function of 17, and Ty/7’,. 
Consequently, a family of curves as in Fig. | may be computed 
for y. 


subject to graphical representation as in Fig. 2 


The right-hand sides of Eqs. (66a) and (71a) as written are 


Hence, to find a value of c;,, or C;,,, the value of y, correspond 
ing to(M,), and (T/T): is found in Fig. 1. In Fig. 2 a line is 
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Graph for local and mean skin-friction coefficients, 





Fic. 2. 


and C;_. 


drawn through (7/7); and (R,,); (on the scale corresponding 
to cy, or Cy,,, as marked) and intersecting the reference line. The 
point,on the reference line is then joined with the value of y; and 
the value of (c;,,); or (C;,,): read on the curve. 


x 


Values for Cy. = h/p..U.¢Cp may be found from the tentative 
equation Cy,, = (1/2)c;, given in the reference or from more 
accurate relations of the form Cy ., = (1/2)c;.,-F( Px) 
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+ 
Additional Values of C(k) 
W. B. Brower and R. H. Lassen 
Aerophysics Laboratory, North American Aviation, Inc., Downey 


Calif. 
October 11, 1952 


T° A RECENT WORK, Luke and Dengler! published, among others, 
a table of the Theodorsen circulation function, C(%) for real 
It was found that this table, while the most com 
> 1.0 with 


arguments. 
plete known to the author, does not cover the range k 
sufficiently small increments. 

Additional values of C(k) were computed from the following 
formulas: 


C(k) = F+4G 


F = JJ: + Yo) + ViCNs Ji 
(Ji + Yo)? + ("1 Jy)? 
G 2 NV Yo + Sid 
(Ji + Yo)? + (¥ Jo)? 


The com 


machines 


using the tables of Bessel functions given by Watson.” 
IBM 
Eight decimals were carried, seven decimals being shown in the 
Table 1 is thought to have the 


putations were performed on punched-card 


result; no ‘‘round-off’’ was made. 
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k 
000 
020 
040 
060 
O80 
100 
120 
140 
160 
180 
200 
220 
240 
260 
280 
300 


.320 
.B40 
_360 
380 


400 
420 
440 


.460 


480 
500 
520 
540 


.560 


580 
600 
620 
640 
660 
680 
700 
720 
740 
760 
780 
800 
820 
840 
860 
880 
900 
920 
940 
960 
980 
000 
020 
040 
060 
O80 
100 
120 
140 
160 
180 
200 
220 
240 
260 
280 
300 
320 
340 
360 
380 
100 
420 
440 
460 
480 
500 
520 
540 
560 
580 





F(k) 
0.539 4348 
0.538 3005 
0.5387 2139 
0.5386 1723 
0.585 1733 
0.5384 2148 
0.533 29438 
0.5382 4099 
0.531 5599 
0.580 7425 
0.529 9560 
0.529 1989 
0.528 4697 
0.527 7672 
0.527 O899 
0.526 4367 
0.525 8066 
0.525 1984 
0.525 6112 
0.524 0439 

52% 4957 

“a 9658 
0.522 4533 
0.521 9574 
0.521 4777 
0.521 0132 
0.520 5634 
0.520 1277 
0.519 7054 
0.519 2961 
0.518 8992 
0.518 5143 
0.518 1407 
0.517 7782 
0.517 42638 
0.517 0845 
0.516 7525 
0.516 43800 
0.516 1165 
0.515 8118 
0.515 5157 
0.515 2274 
0.514 9471 
0.514 6743 
0.514 4087 
0.514 1502 
0.513 8984 
0.513 6533 
0.513 4143 
0.513 1818 
0.512 9547 
0.512 7337 
0.512 5181 
0.512 3077 
0.512 1027 
0.511 9026 
0.511 7073 
0.511 5168 
0.511 3308 
0.511 1490 
0.510 9717 
0.510 7985 
0.510 6293 
0.510 4639 
0.510 3023 
0.510 14438 
0.509 9900 
0.509 8391 
0.509 6914 
0.509 5471 
0.509 4059 
0.509 2678 
0.509 13826 
0.509 0002 
0.508 8707 
0.508 7441 
0.508 6200 
0.508 4985 
0.508 3794 
0.508 2628 





—G(k) 
0.100 2729 
0.098 8716 
0.097 5055 
0.096 1733 
0.094 8739 
0.093 6062 
0.092 3693 
0.091 1620 
0.089 9835 
0.088 8328 
0.087 7089 
0.086 6111 
0.085 5386 
0.084 4904 
0.083 4660 
0.082 4643 
0.081 4849 
0.080 5270 
0.079 5900 
0.078 6732 
0.077 7759 
0.076 8977 
0.076 0380 
0.075 1961 
0.074 3717 
0.073 5641 
0.072 7730 
0.071 9976 
0.071 23878 
0.070 4936 
0.069 7629 
0.069 0469 
0.068 3447 
0.067 6560 
0.066 9803 
0.066 3174 
0.065 6667 
0.065 0282 
0.064 4013 
0.063 7860 
0.063 1817 
0.062 5882 
0.062 0053 
0.061 4326 
0.060 8701 
0.060 3172 
0.059 7739 
0.059 2399 
0.058 7149 
0.058 1988 
0.057 6913 
0.057 1922 
0.056 7013 
0.056 2185 
0.055 74385 
0.055 2762 
0.054 8164 
0.054 3638 
0.053 9184 
0.053 4799 
0.053 04838 
0.053 6234 
0.052 2048 
0.051 7927 
0.051 3868 
0.050 9871 
0.050 5932 
0.050 2058 
0.049 8229 
0.049 4463 
0.049 0750 
0.048 7091 
0.048 3485 
0.047 9930 
0.476 6425 
0.047 2969 
0.046 9561 
0.046 6201 
0.046 2888 
0.045 9619 
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k 
600 
620 
640 
660 
680 
700 
720 
740 
760 
780 
800 
820 
840 
860 
SSO 
900 
920 
940 
960 


980 


000 
020 
040 
060 
O8O 
100 
120 
140 
160 
180 
200 
220 
240 
260 
280 
300 
320 
340 
360 
380 
400 
420 
440 
460 
480 
500 
520 
540 
560 
580 
600 
620 
640 
660 
680 
700 
720 
740 
760 
780 
800 
820 
840 
860 
880 
900 
920 
940 
960 
980 
000 
020 
040 
060 
O80 
100 
120 
140 
160 
180 


F(k) 
0.508 1487 
0.508 0369 
0.507 9271 
0.507 8199 
0.507 7145 
0.507 6114 
0.507 5101 
0.507 4110 
0.507 3136 
0.507 2182 
0.507 1247 
0.507 0828 
0.506 9429 
0.506 8544 
0.506 7678 
0.506 6826 
0.506 5992 
0.506 5170 
0.506 4366 
0.506 3576 
0.506 2799 
0.506 2037 
0.506 1288 
0.506 05538 
0.505 9832 
0.505 9122 
0.505 8424 
0.505 7738 
0.505 7065 
0.505 6401 
0.505 5753 
0.505 5114 
0.505 4485 
0.505 3867 
0.505 3258 
0.505 2660 
0.505 2072 
0.505 1492 
0.505 0925 
0.505 0366 
0.504 9815 
0.504 9273 
0.504 8742 
0.504 8216 
0.504 7701 
0.504 7192 
0.504 6693 
0.504 6200 
0.504 5716 
0.504 5237 
0.504 4768 
0.504 4306 
0.504 3851 
0.504 3402 
0.504 2960 
0.504 2525 
0.504 2096 
0.504 1672 
0.504 1257 
0.504 O845 
0.504 0445 
0.504 0044 
0.503 9652 
0.503 9264 
0.503 8886 
0.503 8509 
0.503 8140 
0.503 7774 
0.503 7413 
0.503 7059 
0.503 6708 
0.503 6364 
0.503 6023 
0.503 5688 
0.508 5359 
0.503 5029 
0.503 4708 
0.503 4391 
0.503 4078 
0.503 3770 


—G(k) 
0.045 6393 
0.045 32138 
0.045 0075 
(.044 6980 
0.044 3925 
0.044 0909 
0.043 7985 
0.043 4999 
0.043 2100 
0.042 9239 
0.042 6415 
0.042 3626 
0.042 0873 
0.041 8154 
0.041 5468 
0.041 2818 
0.041 O199 
0.040 7612 
0.040 5057 
0.040 25383 
0.040 00388 
0.039 7575 
0.039 5142 
0.039 2736 
0.039 0359 
0.038 8011 
0.038 5689 
0.038 3395 
0.038 1126 
0.037 8884 
0.037 6668 
0.087 4477 
0.037 2310 
0.037 0168 
0.036 8049 
0.086 5954 
0.036 3883 
0.036 1834 
0.035 9809 
0.035 7803 
0.035 5820 
0.035 3860 
0.0385 1919 
0.0385 0001 
0.034 8101 
0.034 6223 
0.034 4365 
0.034 2524 
0.034 07038 
0.033 8901 
0.083 7119 
0.033 5353 
0.033 3608 
0.033 1879 
0.033 0167 
0.032 8470 
0.082 6793 
0.0382 5133 
0.032 3490 
0.032 1860 
0.032 0248 
0.081 8653 
0.031 7072 
0.031 5507 
0.031 3957 
0.031 2421 
0.031 0901 
0.030 9394 
0.030 7901 
0.030 6424 
0.0380 4961 
0.030 3509 
0.030 2074 
0.080 0650 
0.029 9239 
0.029 78438 
0.029 6459 
0.029 5084 
().029 3725 
0.029 2377 
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TABLE 1 (Continued) 











k Fik) —G(k) k F(k) -G(k) 
4.200 0.503 3465 0.029 1044 6.000 0.501 6845 0.020 5954 | Sub: 
4.220 0.503 3164 0.028 9719 6.100 0.501 6312 0.020 2651 | foil th 
4.240 0.503 2867 0.028 8408 6.200 0.501 5808 0.019 9448 | tostre 
4.260 0.503 2574 0.028 7108 6.300 0.501 5320 0.019 6349 
4.280 0.503 2284 0.028 5819 6.400 0.501 4854 0.019 3339 
4.300 0.503 2000 0.028 45438 6.500 0.501 4415 0.019 0420 Cy. * 
4.320 0.508 1718 0.028 3276 6.600 0.501 3994 0.018 7589 
4.340 0.508 1488 0.028 2022 6.700 0.501 3588 0.018 4840 
4.360 0.5038 1166 0.028 0778 6.800 0.501 3199 0.018 2169 
4.380 0.503 O895 0.027 9545 6.900 0.501 2824 0.017 9572 
4.400 0.503 0627 0.027 8322 7.000 0.501 2468 0.017 7050 
4.420 0.503 O0362 0.027 7110 7.100 0.501 2130 0.017 4595 ~ 
4.440 0.508 0104 0.027 5908 7.200 0.501 1800 0.017 2209 
4.460 0.502 9845 0.027 4717 7.300 0.501 1486 0.016 9887 8 
4+. 480 0.502 9592 0.027 3534 7.400 0.501 1186 0.016 7625 saith 
4.500 0.502 9340 0.027 2363 7.500 0.501 O894 0.016 5423 
4.520 0.502 9093 0.027 1200 7.600 0.501 0613 0.016 3278 C2 
4.540 0.502 8848 0.027 0048 7.700 0.501 0348 0.016 1184 
4.560 0.502 8609 0.026 8906 7.800 0.501 O089 0.015 9147 
4.580 0.502 8869 0.027 7775 7.900 0.500 9837 0.15 7160 
4.600 0.502 8131 0.026 6651 8.000 0.500 9598 0.015 5220 
4.620 0.502 7901 0.026 5537 8.100 0.500 9364 0.015 3329 
4.640 0.502 7670 0.026 4431 8.200 0.500 9142 0.015 1484 “ 
4.660 0.502 7443 0.026 3336 8.300 0.500 8923 0.014 9677 
+. 680 0.502 7222 0.026 2250 8.400 0.500 8718 0.014 7918 
a _ 7 - os - ° va for 4 
+. 700 0.502 7000 0.026 1170 8.500 0.500 8518 0.014 6197 ae 
$.720 0.502 6781 0.026 0100 8.600 0.500 8323 0.014 4518 a 
4.740 0.502 6564 0.025 9040 8.700 0.500 8136 0.014 2873 the e 
+. 760 0.502 6350 0.025 7987 8.800 0.500 7954 0.014 1269 condi 
+. 780 0.502 6137 0.025 6948 8.900 0.500 7779 0.013 9695 eed 
AU 
+. 800 0.502 5930 0.025 5907 9.000 0.500 7607 0.018 8159 
+. 820 0.502 5722 0.025 4877 9.100 0.500 7443 0.013 6658 value 
4.840 0.502 5519 0.025 3861 9.200 0.500 7288 0.013 5187 Th 
+. 860 0.502 5318 0.025 2848 9.300 0.500 7133 0.013 3746 mann 
4+. 880 0.502 5119 0.025 1846 9.400 0.500 6981 0.013 2339 
4.900 0.502 4923 0.025 0849 9.500 0.500 6842 0.013 0960 
4.920 0.502 4729 0.024 9862 9.600 0.500 6698 0.012 9605 
4 940 0.502 4536 0.024 8880 9.700 0.500 6563 0.012 8285 
4.960 0.502 4348 0.024 7909 9.800 0 500 6430 0.012 6986 wher 
4.980 0.502 4158 0.024 6941 9.900 0.500 6801 0.012 5711 
& . 9 ART . 
5.000 0.502 3974 0.024 5985 10.000 0.500 6180 0.012 4467 C 
Ps Sle tb ; riety 10.500 0.500 5611 0.011 8586 
5.100 0.502 3076 0.024 1303 = “ean Pet , 
ae =i tat dg ae. aes 11.000 0.500 5121 0.011 3234 (2 
5.200 0.502 2228 0.023 6794 el = “ee = tt 
eas ptt ae tlle Be 11.500 0.500 4684 0.010 8343 
5.300 0.502 1426 0.023 2451 9 ~ > F : 
v2 nts se =p dics 12.000 0.500 4304 0.010 ( 
5.400 0.502 0663 0.022 8259 a oe Popa 
ey ¥ d 12.500 0.500 3972 0.009 
5.500 0.501 9946 0.022 4216 13.000 0.500 3673 0.009 5905 
5.600 0.501 9260 0.022 0311 13.500 0.500 3410 0.009 2373 
5.700 0.501 8609 0.021 6541 14.000 0.500 3170 0.008 9091 
5. 800 0.501 7992 0.021 2895 14.500 0.500 2952 0.008 6030 
5.900 0.501 7404 0.020 9369 15.000 0.500 2766 0.008 3174 D 
15.500 0.500 2589 0.008 0502 
16.000 0.500 2433 0.007 7993 
| 
: 4 ’ - For | 
same estimated accuracy as those of Luke and Dengler—i.e., the REFERENCES 
error is £3 X 107%. 1 Luke, Y. L., and Dengler, M. A., Tables of the Theodorsen Circulation VM 
With this table supplementing those of Luke and Dengler and — F'"¢lion for Generalized Motion, Journal of the Aeronautical Sciences, Vol 
‘ ‘ - F P = ae 7 18, No. 7, p. 481, July, 1951 ‘ 
using the three-point Lagrangian interpolation table,’ it is possible 2 Watson, G. N.. A Treatise on the Theory of Bessel Functions, 2nd Ed., P for 4 
. : . : . : 666; Cambridge University Press, London, 1945 Ther 
“ ‘ 4 r k yiven s1X -Ccimie w I < - é - » 4 S : 2 
to interpolate fo given to six decimals with an estimated error : Tables of Legrangien Iuterpolation Cocficients; Columbia Universit 
te 5 xX 10-6 Press, 1944 ( 
C; 
+ ' 
( 
Comments on ‘‘Two-Dimensional Airfoils at ¢,, 2 + +1 ¥+1 C 
Moderate Hypersonic Velocities’’ =" sft = 7 * —. 
. The 
L.E. Flanagan, Jr. (y+) Y) K2fs 4 . (la as 
Research Engineer, Aerophysics Laboratory, North American 18 : 
Aviation, Inc., Downey, Calif F : C,, 
‘or a compression, 
October 27, 1952 é, } 
c 2 vy+1 (y + 1)? 
“= —f 4 i4 = Kf 4 
6? K 2 16 
yy HIS RECENT PAPER ON HYPERSONIC AIRFOILS, Dorrance! has : 
developed power series for pressure coefficients for an expan (y¥ +) K3f6 4 Ib 
1,024 


sion and a nonisentropic compression. For an expansion, 
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Substituting the relations K = Méand f = 6/6, where 6 = air- 


foil thickness ratio and 6 = angle of surface 


tostream, in Eqs. (la ) and (1b) gives 


9 wy +4 3 + 1 
Cc. =—0+- e+ Me 4 
pe M 2 6 
( + 1) (8 ) 
Z —— 
48 
- . 7 1.2 @2 + 0.40 Me O}107 1A? |a4 
M 
for y = 1.4, and 
2 + (y + 1) 
C= c= e+ — Ve* 4 
; WV 2 16 
(y 4+ F 
: M96 
1,024 
2 
= 6+ 1.26? + 0.36.1/88 + O[10~2.173 |@® 
M 
for y = 1.4 


+ 


+ 


inclination relative 


(2b) 


It is the purpose of this note to compare these equations with 


the expressions developed by Busemann? and Miles* under the 


condition that VM? — 1/M = 1 


Number for which Dorrance’s equations apply is 3.19. 


value of 17, V M? 1/M = 0.95 


The lower limit 


of Mach 
At this 


The expressions for pressure coefficients as developed by Buse 


mann and corrected by others are 


Cre = C10 + Co? + C08 4 
“oe Ce + GC + (C D)e* 4 
where 
C, = 2/V M? -— 1 
C, = [yA + (MM? 2)?2|/2( M2 1)? 
l 1. 3 
C; = : us 4 
(M 1)? ’ I if) 
27? ; (y+1 
- - oS i 
6 3 
l 7 + l 
D = M x 
(M2 — 1)7/Q 5 12 ) 


For hypersonic flow, 


VM*-—-12M: (M 2)? = M*: 
(Af? 1)¥/' 
for y = 1.4. 
Therefore, 
C, ~2/M 
C; “(y+ 1)/2 
y+] 2y2—7y —5 
( ne Ee i eed ag 
6 6 
; (y 4 5y2 6 3 
( D) Y 1) v4 rY l6y a 
16 12 
The pressure coefficients then become 
9 
2 ¥ + 1 Y 
l,, = 6 + ge? + \7e* 4 
I rT, r ; M 
ay” = *9 
6 
2 
= 6+ 1.267 0.40170 O[ AM ~' a8 


O[M 3] 


+ 


O|M 
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for y = 1.4, and 
: 2 y¥+1 (7 1)? 
Cc, = 6 + a 4 Me + 

M 2 16 

YY 164 13 
12 Ww '0 (4b) 
2 
= 6 1.202 + 0.36 1768 O|M~'\6 

M 

for y = 14 


These equations agree with those of Dorrance for the first three 


terms. 
Miles has developed the following expressions for the pressure 


coefficients, for y = 1.4: 
: 2 1O(3.474 5M? + 5 
C,, = T - - “4 4 
VV 1 25( M? 1)? 
10(15.M* 68.\7° + 150174 75M? + 50) 
+ 
375( M? 1) 
. 2 2(3.M‘i — 5M? + 5) 
Cre = 64 - e +4 
V Me 1 5( MM? 1)? 
27 M8 112M* + 270M! 1504/2 + 100 
75( M? 1) 
Applying the conditions of Eqs. (3) gives 
Coe = (2/M)O + 1.262 — O| M~? 6? + 0.401108 — O[M~']08 § (5a) 
Cy, & (2/M)0 + 1.262 ~ O[.M-2]e? + 0.36.Me* — O[.M-1J6* (5b) 


The effect of the terms of the order of 1/~*6? may be determined 
For hypersonic flow, 1/6 is of the order of 1, and 6 is of 
Therefore, @ is of the order of 1J~', sothat \/~°6? is 


readily. 
the order of 6 


of the order of 17~'6°. Thus, Eqs. (5a) and (5b) become 
Cp. = (2/M)0 + 1.267 + 0.40.Me O[M~')e8 
Cp, = (2/M)0 + 1.26? + 0.36.M05 O[M~—' a 


It is apparent that the first three terms of Miles’ equations agree 
with the corresponding terms of Busemann and of Dorrance and 
that the remaining terms are of the same orders in 6 and M as 
are those of Busemann. 

It should be pointed out that the fourth terms of Eqs. (2a) and 
(2b), being of the order of 10~'A7°6* and 1072.65, are, under the 
condition that M@ is of the order of 1, of the second order in 6 
All subsequent terms of these series developed by Dorrance are 
also of the order of 6?, but, since the series converge rapidly, it is 
practical, as Dorrance points out, to neglect all terms beyond the 


third in the expressions for C,, and C,, 
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On the Transition Through a Weak Shock 
Front* 


Paul Mandl 
Institute of Aerophysics, University of Toronto, Toronto, Canada 
September 29, 19592 


TRANSITION OF THE one-dimensional steady flow of real 


r | SHE 
fluids through shock fronts has received much attention re 


* The complete report is available at the Institute of Aerophysics, Uni 


versity of Toronto, Toronto (U.T.I.A. Report No. 14 
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cently. Becker’s exact solution for a Prandtl Number of */;,' has 
been extended by several authors;” * approximate solutions have 
been derived for other Prandtl! Numbers.‘ 

This investigation deals rigorously with the shock transition of 
a perfect gas passing through a weak shock front for arbitrary 
values of the Prandtl Number. The mathematical analysis is 
based on the differential equation describing the transition of ve 
locity which can be deduced readily from the general equations of 
motion of a continuous fluid. If the coefficients of viscosity and 
heat conduction are assumed constant, then the basic equation 


for the velocity can be expressed as 


att ‘ (<2 2 vr -_ dU ’ 
- Ca = ) 
dkR? dR dR 


C(U 2 


1)(U — Uz =0 (1) 


where U is a velocity ratio equal to unity at the initial state and 
Uz, at the final state; R is a shock-wave Reynolds Number pro 
portional to the distance normal to the shock front; and A, B, C 
are constants depending on the free-stream Mach Number J/,, 
the ratio of specific heats of the gas y and the Prandtl Number a. 
The condition of uniform flow both ahead of, and behind, the 
shock front requires that the solution of Eq. (1) satisfy the fol- 
lowing boundary conditions: 
Ui! 


U— U: 


as R—~ — o-, 
asR~+0, 
The additional condition that d?U/dR? = 0 when R = 0 de 
termines the position of the centre of the shock front. 
If, and only if, the constants of the gas satisfy the relation 


(o — 444) (vy — 1) = 0 


Eq. (1) can be factorized into the differential operators 


d 2A 
— B x 
ee”) 


aU A : 
[ - nam 1)(U — U2) 
dR 1+ U, 


Eq. (3) is satisfied if the second factor equals zero. 
tion can be integrated at once, yielding a solution in a closed form 


This equa 


that satisfies the above boundary conditions—viz., 


Lit Bl (-7) (4-4) 
- | log - UU» log 
434-8 ti ae ft. — 0, 


(4) 


R= 


Hence, the basic equation for the velocity degenerates to a 
first order equation if o = */, [cf. Eq. (2); the value y = 1 must 
be rejected for a real gas]. This result coincides with the solu- 
tion Becker obtained by entirely different reasoning.! 

To solve Eq. (1) for arbitrary values of o, the substitution dU 
dR = n transforms it into an equation of the hodograph type 


_ dn 
"20 


The corresponding boundary conditions become 


l +? +(A — BUn + C(U 1)(U — Us) = 0 (5) 


n(1) = n( U2) = 0 
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The initial state is a nodal point and the final state a saddle 
point of Eq. (5). An infinite number of integral curves Passes 
through the node but only two through the saddle. The problem 
is, therefore, to find a curve connecting the neighborhoods of these 
two singular points. 

If calculations are started from the initial state, an integral 
of Eq. (5) can be expressed in terms of the following series: 


m mn 7 yb 
>. Cmy & + > S. tnt (CEP (6) 
m ln=1 


m= 1 


3 = 


where — = (1 — U)/(1 U2); bisa constant that depends on the 
physical quantities y, 7, M, and is generally irrational; and Cjg 
an arbitrary constant characterizing the single infinity of solutions 
through the node. 

The coefficients Cm, are determined from the differential equa. 
tion but are of such a form that the solution of Eq. (6) breaks 
down for a discrete set of Mach Numbers; the series must then 
be modified by including logarithmic terms. 

Because of the presence of the arbitrary constant C, the series 
solution (6) is not, in general, suitable for numerical calculations, 
However, for infinitely weak shock fronts ( Cc| + 0), Eq. (6) re- 
duces to an ordinary power series which converges slowly. (If 
only the first two terms are retained, the solution is identical with 
Taylor’s approximate formula. )® 

Since the final state is a saddle point of the differential equa- 
tion, a ‘‘backwards”’ solution should yield a more satisfactory re- 
In fact the solution can always be expressed in terms of an 
( The rate of con- 


sult. 
ordinary power series, which converges rapidly 
vergence increases as 1, > 1.) 

This power series was used to calculate the variation of ve- 
locity through the shock front for various values of M, and o. 
The following results were obtained: 

(1) Increasing the Prandtl Number at a constant Mach Number 
steepens the local velocity gradients and, hence, shortens the 
transition zone. 

(2) In the neighborhood of « = */;, the results are relatively in- 
sensitive to small deviations in o, particularly at Mach Numbers 
near 1. Therefore, Becker’s results can be applied in this range, 
though caution should be used at higher Mach Numbers. 

(3) Taylor’s solution® predicts slightly high values for the di- 
mensions of the transition zone. 

(4) The present results agree with Mises’ approximate values 
for the thickness of the shock front.‘ 
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